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Foreword 


Once  again,  it  is  a  great  pleasure  to 
report  on  a  year  of  accomplish- 
ments in  all  areas  of  the  Center.  Although  the 
shutdown  was  prolonged  by  unforeseen  problems 
in  an  auxiliary  system,  the  results  when  we  re- 
started were  everything  that  we  had  hoped  for.  The 
new  cooling  tower  is  functioning  well,  reducing  the 
visible  plume  during  cold  weather  dramatically.  This 
system  is  now  well  prepared  for  another  twenty-plus 
years  of  operation.  The  work  on  relicensing  is  also  on 
track  —  the  new  analyses  are  almost  complete,  and 
confirm  the  existing  Safety  Analysis  Report.  A  new  liquid 
hydrogen  cold  source  of  improved  design  was  installed 
during  the  same  shutdown.  It  is  working  well,  producing 
nearly  twice  as  many  cold  neutrons  as  its  predecessor, 
exactly  in  agreement  with  the  detailed  MCNP  predictions. 
In  every  respect,  the  future  for  neutron  production  at  the 
NIST  source  looks  excellent. 

As  always,  we  continue  to  improve  existing  instru- 
ments and  work  on  new  ideas.  An  optical  filter  was 
installed  in  the  NG-3  30  m  SANS  instrument  (a  key  part  of 
the  NIST/NSF  Center  for  High  Resolution  Neutron 
Scattering),  providing  further  gains  of  up  to  a  factor  of 
three  at  long  wavelength.  A  modification  to  the  guide  that 
runs  through  the  choppers  of  the  Disk  Chopper  time-of- 
flight  Spectrometer  has  given  intensity  increases  of  2  to  3 
over  the  whole  wavelength  range.  So  these  two  instru- 
ments have  increased  in  performance  by  factors  of 
between  3  and  5  over  the  past  year,  a  direct  result  of  our 
continuing  investment  in  instrumentation.  A  consortium  of 
researchers  from  five  universities  has  been  awarded  a 
grant  from  the  National  Institutes  of  Health  to  begin 
development  of  a  new  Research  Resource  called  Cold 
Neutrons  for  Biology  and  Technology,  which  includes 
construction  of  a  new  diffractometer/reflectometer  for 
membrane  and  related  structures.  The  National  Science 
Foundation  has  increased  its  partnership  with  the  NCNR, 
and  now  helps  to  support  the  high-resolution  inelastic 
scattering  instruments  in  the  guide  hall.  These  and  other 


partnerships,  such  as  the  ones  with  Johns  Hopkins  Univer- 
sity, the  University  of  Maryland,  and  Brookhaven  National 
Laboratory,  are  a  key  part  of  our  plans  to  move  forward, 
and  have  been  an  essential  part  of  our  past  success. 

And  again,  as  always,  the  science  being  produced 
fully  justifies  the  dedication  and  work  of  the  people  who 
provide  the  neutrons  and  instruments.  While  walking 
through  the  facility  during  operating  cycles,  I  am  always 
amazed  at  the  breadth  and  depth  of  research  now  being 
done  with  neutrons,  and  delighted  by  the  energy  and 
intensity  of  the  young  people  doing  it.  I  derive  great 
pleasure  helping  to  choose  scientific  highlights,  and  then 
from  reading  the  results  of  those  choices.  They  form  the 
greatest  portion  of  this  report,  and  I  hope  that  you  enjoy 
them  as  much  as  I  do. 


iv 


TheNIST  Center  for 
Neutron  Research  (NCNR) 


Neutrons  are  powerful  probes  of  the  structure  and  dy- 
namics of  materials  ranging  from  molecules  inserted 
into  membranes  mimicking  cell  walls  to  protons  migrating 
through  fuel  cells.  The  unique  properties  of  neutrons  (see 
inset)  can  be  exploited  by  a  variety  of  measurement  tech- 
niques to  provide  information  not  available  by 
other  means.  They  are  particularly  well  suited 
to  investigate  all  forms  of  magnetic  materials 
such  as  those  used  in  computer  memory  stor- 
age and  retrieval.  Atomic  motion,  especially 
that  of  hydrogen,  can  be  measured  and  moni- 
tored, like  that  of  water  during  the  setting  of 
cement.  Residual  sfresses  such  as  those  in- 
side stamped  steel  automobile  parts  can  be 
mapped.  Neutron-based  research  covers  a 
broad  spectrum  of  disciplines,  including  en- 
gineering, biology,  materials  science,  poly- 
mers, chemistry,  and  physics. 

The  NCNR's  neutron  source  provides 
the  intense  beams  of  neutrons  required  for 
these  types  of  measurements.  In  addition  to 
the  thermal  energy  neutron  beams  from  the 
heavy  water  or  graphite  moderators,  the 
NCNR  has  a  large  area  liquid  hydrogen  mod- 
erator, or  cold  source,  that  provides  intense 
neutron  beams  for  the  only  cold  neutron  fa- 
cility presently  operating  in  the  United  States. 

There  are  currently  28  experiment  stations:  six  pro- 
vide high  neutron  flux  positions  for  irradiation,  and  22  are 
beam  facilities  most  of  which  are  used  for  neutron  scattering 


Instruments  in  the  cold  neutron  guide  hall. 


research.  The  following  pages  show  a  schematic  layout  of 
the  beam  facilities.  More  complete  descriptions  of  instru- 
ments can  be  found  at  www.ncnr.nist.gov. 

The  NCNR  supports  important  NIST  research  needs, 
but  is  also  operated  as  a  major  national  user  facility  with 

Why  Neutrons? 

Neutrons  at  the  NCNR  reveal  properties  not  available  to  other  probes.  They  can  behave  like  micro- 
scopic magnets,  can  diffract  like  waves,  or  set  particles  into  motion  as  they  recoil  from  them. 

Wavelengths  -  range  from  0. 1 A  to  1 00  A,  allowing  them  to  form  observable  ripple  patterns  from 
structures  as  small  as  atoms  to  as  large  as  proteins. 

Energies  -  of  millielectronvolts,  the  same  as  that  of  motions  of  atoms  in  solids  or  liquids,  waves  in 
magnetic  materials,  or  vibrations  in  molecules.  Exchanges  of  energy  between  neutrons  and  matter  as 
small  as  nanoelectronvolts  and  as  large  as  tenths  of  electronvolts  can  be  detected. 

Selectivity  -  in  scattering  power  varies  from  nucleus  to  nucleus  almost  randomly.  Specific  isotopes 
can  stand  out  from  other  isotopes,  even  of  the  same  kind  of  atom.  Specific  light  atoms,  difficult  to  observe 
with  x-rays,  are  revealed  by  neutrons.  Hydrogen,  especially,  can  be  distinguished  from  chemically 
equivalent  deuterium. 

Magnetism  -  makes  the  neutron  sensitive  to  the  magnetic  spins  of  both  nuclei  and  electrons,  allowing 
the  behavior  of  ordinary  and  exotic  magnets  to  be  detailed  precisely. 

Neutrality  -  of  the  uncharged  neutrons  allows  them  to  penetrate  deeply  without  destroying  samples,  and 
pass  through  walls  controlling  a  sample's  environment  allowing  measurements  under  extreme  conditions. 
Properties  ranging  from  the  residual  stresses  in  steel  girders  to  the  unfolding  motions  of  proteins  are 
amenable  to  measurement  by  neutrons. 

Capture -characteristic  radiation  emanating  from  specific  nuclei  capturing  incident  neutrons  can  be  used 
to  identify  and  quantify  minute  amounts  of  material  in  pollutants  or  ancient  pottery  shards. 


merit-based  access  made  available  to  the  entire  U.S.  techno- 
logical community.  Each  year,  over  1700  research  partici- 
pants from  all  areas  of  the  country,  from  industry,  academia, 
and  government  use  the  facility  for  measurements.  Beam 
time  for  research  to  be  published  in  the  open  literature  is 
without  cost  to  the  user,  but  frail  operating  costs  are  recov- 
ered for  proprietary  research.  Access  is  gained  mainly  through 
a  peer-reviewed,  web-based  proposal  system  with  beam  time 
allocated  by  a  Program  Advisory  Committee  twice  a  year. 
For  details  see  www.ncnr.nist.gov/beamtime.html.  The  Na- 
tional Science  Foundation  and  NIST  co-fund  the  Center  for 
High  Resolution  Neutron  Scattering  (CHRNS)  that  operates 
six  of  the  world's  most  advanced  instruments.  Time  on 
CHRNS  instruments  is  made  available  through  the  proposal 
system.  Some  access  to  beam  time  for  collaborative  mea- 
surements with  the  NIST  science  staff  can  also  be  arranged 
on  other  instruments. 
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A  Cold  Neutron  Depth  3 
Profiling  instrument  (not 
shown)  for  quantitative  profiling 
of  subsurface  impurities 
currently  at  this  site  will  be  ^ 
moved  to  another  position. 
Shown  is  MACS,  a  triple  axis 
cold  neutron  crystal  spectrom- 
eter under  construction  with 
double  focusing  monochromator  ^ 
and  multiple  crystal  analyzer/ 
detectors  that  can  be  flexibly 
configured  for  several  energies 
simultaneously  or  for  high 
throughput  at  one  energy.  6 

BT-7  Triple  Axis 
Spectrometer  with  fixed 
incident  energy  for  measure- 
ments of  excitations  and 
structure. 


BT-8  Residual  Stress 
Diffractometer  optimized  for 
depth  profiling  of  residual  stress 
in  large  components. 

BT-9  Triple  Axis  Crystal 
Spectrometer  for  measure- 
ments of  excitations  and 
structure. 

Thermal  Column  a  very  well- 
thermalized  beam  of  neutrons 
used  for  radiography,  tomogra- 
phy, dosimetry  and  other 
experiments. 

BT-1  Powder  Diffractometer 

with  32  detectors;  incident 
wavelengths  of  0.208  nm,  0.154 
nm,  and  0.159  nm,  with  highest 
resolution  of  8d/d  =  8x10-4. 


7  BT-2  Triple  Axis  Crystal 
Spectrometer  with  polarized 
beam  capability  for  measure- 
ment of  magnetic  dynamics  and 
stmcture. 

8  BT-4  Filter  Analyzer  Neutron 
Spectrometer  with  cooled 
Be/Graphite  filter  analyzer  for 
chemical  spectroscopy. 

9  BT-5  Perfect  Crystal 
Diffractometer  SANS  small 
angle  neutron  scattering 
instrument  for  microstructure  on 
the  W  nm  length  scale, 
sponsored  by  the  National 
Science  Foundation  and  NISI, 
part  of  the  Center  for  High 
Resolution  Neutron  Scattering 
(CHRNS). 


10  NG-7  Horizontal  Sample 
Reflectometer  allows 
reflectivity  measurements  of  free 
surfaces,  liquid  vapor  interfaces, 
as  well  as  polymer  coatings. 

11  Neutron  Interferometry  and 
Optics  Station  with  perfect 
silicon  interferometer;  vibration 
isolation  system  provides 
exceptional  phase  stability  and 
fringe  visibility. 

12  Spin  Polarized  Triple  Axis 
Spectrometer  (SPINS)  using 
cold  neutrons  with  position 
sensitive  detector  capability  for 
high  resolution  studies  —  part  of 
CHRNS. 
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NIST  Center  for  Neutron  Research  Layout 


16 


15 


13  Spin  Echo  Spectrometer 

offering  neV  energy  resolution, 
based  upon  Julicli  design, 
sponsored  by  NIST,  Julich  and 
ExxonMobil  —  part  of  CHRNS. 

14  Prompt  Gamma  Activation 
Analysis  cold  neutron  fluxes 
allow  detection  limit  for  H  of 

1  |ig  to  10  ng.  Focused  beams 
are  available  for  profiling. 

15  NG-7  30  m  SANS  for 

microstojcture  measurements 
sponsored  by  NIST,  ExxonMobil, 
and  the  University  of  Minnesota. 

16  Neutron  Physics  Station 

offering  three  cold  neutron  beams 
having  wavelengths  of  0.5  nm, 
0.9  nm,  and  "white"  that  are 
available  for  fundamental  neutron 
physics  experiments. 


17  Fermi  Chopper  hybrid  time- 
of-flight  (TOF)  Spectrometer 

for  inelastic  scattering  with 
incident  wavelengths  between 
0.23  nm  and  0.61  nm  chosen 
by  focusing  pyrolytic  graphite 
crystals.  A  simple  Fenni 
chopper  pulses  the  beam. 

18  Disk  Chopper  TOF  Spec- 
trometer a  versatile  time-of- 
flight  spectrometer,  with  beam 
pulsing  and  monochromatization 
effected  by  7  disk  choppers. 
Used  for  studies  of  dynamics  in 
condensed  matter,  including 
macromolecular  systems  —  part 
of  CHRNS. 


19  NG-3  30  m  SANS  for 

microstructure  measurements 
sponsored  by  the  National 
Science  Foundation  and  NIST 
-  part  of  CHRNS. 

20  Bacl(scattering 
Spectrometer  high  intensity 
inelastic  scattering  instrument 
with  energy  resolution  <  1  i^eV, 
for  studies  of  motion  in  molecular 
and  biological  systems  —  part  of 
CHRNS. 

21  NG-1  8  m  SANS  for  polymer 
characterization  being  modified 
to  10  m  and  to  be  made 
available  for  CHRNS  use  along 
with  its  current  use  by  the  NIST 
Polymers  Division. 


22  Vertical  Sample 
Reflectometer  instrument  with 
polarization  analysis  capability 
for  measuring  reflectivities  down 
to  10-8  to  determine  subsurface 
structure. 

23  Cold  Neutrons  for  Biology 
and  Technology  vertical 
sample  reflectometer,  to  be  built 
with  polarization  analysis 
capability  for  measuring 
reflectivities  down  to  10-^.  It 
will  have  a  position- 
sensitive  detector  for 
measuring  off- 
specular 
reflections. 


Removal  of  Unit  1 

The  bear  trap  is  placed  on  special 


The  "bear  trap"  transport  device  is  moved  into  place 
to  remove  Unit  1. 


Installation  of  Unit  2 
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Cold  Source  Replacement 


The  Advanced  Liquid 
Hydrogen  Cold  Neutron  Source 


The  14  cold  neutron  instruments  at 
the  NCNR  received  a  substantial 
increase  in  performance  when  the  liquid 
hydrogen  cold  source,  successfully  operated 
since  1995  (hereafter  called  "Unit  1"),  was  replaced 
with  an  improved  design  (called  "Unit  2")  during  last 
year's  maintenance  outage.  The  NCNR  cold  source  is 
designed  to  maximize  the  useftil  intensity  of  neutrons  in 
the  wavelength  band  between  0.2  nm  and  2.0  nm  by 
providing  a  low  temperature  medium  that  further  slows 
down  the  neutrons  exiting  the  heavy  water  moderator/ 
reflector  of  the  NIST  20  MW  reactor.  The  basic  operating 
principles  for  the  cold  source  are  the  same  in  both  designs 
[1].  Using  the  large  diameter  cryogenic  beam  port  at  the 
NCNR  reactor,  a  cryostat  containing  a  thin,  annular  shell  of 
liquid  hydrogen  is  placed  in  a  high  intensity  region  of  the 
reactor's  reflector.  Interaction  with  the  slow  moving 
hydrogen  molecules  modifies  the  neutron  energy  spectrum. 
A  neutron  guide  tube  system  beginning  close  to  the  exit  of 
the  cryostat  transports  the  intense  cold  neutron  beams  to 
the  neutron  scattering  instruments.  The  radiation  heat 
deposited  in  the  moderator  is  removed  by  boiling  of  the 
liquid  hydrogen.  Stable  operation  of  the  cold  source  is 
maintained  by  a  thermosiphon;  cold  liquid  hydrogen  flows 
by  gravity  into  the  cryostat  from  a  helium  cooled  con- 
denser and  the  hydrogen  vapor  from  boiling  returns  to  the 
condenser  in  a  closed,  convective  circulation  loop.  The 
cryostat  also  includes  a  heavy  water  cooling  jacket, 
insulating  vacuum  regions  and  helium  containment  vessel 
that  surrounds  all  portions  of  the  system  containing 
hydrogen. 

Following  the  commissioning  of  Unit  1,  detailed 
simulations  of  the  NCNR  reactor  and  moderators  became 
sufficiently  accurate  that  it  was  possible  to  explore  further 
optimization  of  the  cold  source  geometry  for  improved 
performance  using  a  sophisticated  MCNP  model  of  the 
NCNR  source  [2].  The  Unit  2  design  differs  from  Unit  1 
in  many  key  respects,  as  can  be  seen  in  Fig.  1 .  Most  of 
the  volume  of  the  Unit  1  cryostat  assembly  was  empty 
space,  leaving  a  large  void  in  the  reactor  reflector  and 
lowering  the  neufron  flux  near  the  cryostat.  The  heavy 
water  cooling  jacket  of  Unit  2  was  expanded  to  fill  nearly 
half  of  the  insulating  vacuum  void  of  Unit  1.  This  change 


reduced  the  void  in  the  reflector  and  increased  the  neutron 
flux  near  the  cryostat  by  about  40  %  thus  providing  a  cold 
neutron  gain  of  the  same  magnitude.  Because  of  these 
changes  and  to  fully  illuminate  the  rectangular  aspect  ratio 
of  the  NCNR  guide  system,  the  Unit  2  moderator  chamber 
is  smaller  than  the  original  and  is  an  ellipsoidal  annulus, 
rather  than  a  spherical  one.  With  a  smaller  vessel  size,  the 
thickness  of  the  liquid  hydrogen  shell  could  be  increased 
from  20  mm  to  between  25  mm  and  30  mm  without 
increasing  the  overall  hydrogen  inventory.  Finally,  the 
inner  ellipsoid  of  the  newer  design  is  evacuated  through  a 
small  vacuum  port,  removing  most  of  the  hydrogen  vapor. 
In  Unit  2,  the  cold  neutron  beam  passes  through  only 
20  mm  of  vapor,  rather  than  through  the  300  mm  path 

Liquid  Hydrogen 


Fig.  1.  Comparison  between  tlie  original  (top)  and  the  advanced  liquid 
hydrogen  cold  sources  for  the  NIST  research  reactor  The  heavy  water 
volume  is  indicated  in  blue,  liquid  hydrogen  in  orange,  and  hydrogen 
vapor  in  yellow. 
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length  in  Unit  1 .  This  resuhs  in  less  neutron  adsorption  by 
hydrogen.  Collectively,  these  changes  in  the  moderator 
chamber  add  another  15  %  to  20  %  to  the  predicted  gain. 

The  complex  geometry  of  Unit  2  posed  many 
challenges  for  the  design  and  fabrication  of  the  cryostat. 
Complicated  elliptical  vessels  were  designed  with  the  aid  of 
finite  element  analysis  to  model  their  mechanical  strength. 
Extensive  pressure  testing  supported  this  detailed  modeling. 
Because  of  the  non-standard  shapes,  components  were 
precisely  cut  from  solid  blocks  of  Aluminum-6061  using  a 
high-speed  mill  at  the  NIST  instrumentation  shop.  The 
geometry  also  placed  high  demands  on  the  welding  of  the 
various  components  into  a  leak-tight  cryostat  satisfying 
stringent  quality  control  tests.  The  entire  cryostat  system 
was  completely  fabricated,  assembled,  and  tested  at  NIST. 
The  replacement  of  Unit  1  with  the  new  design  was  also  a 
significant  technical  challenge.  Detailed  work  planning 
was  required  to  ensure  that  the  highly  radioactive  Unit  1 
could  be  safely  removed  with  the  least  amount  of  disrup- 
tion to  the  sensitive  alignment  of  the  neutron  guide  net- 
work. Radiation  damaged  neutron  guides  near  the  source 
were  also  systematically  replaced  during  the  installation. 
The  installation  was  completed  on  schedule  with  a  minimal 
dose  exposure  to  the  installation  team. 

When  the  reactor  resumed  normal  operation  on 
March  6,  2002,  the  thermal  and  neutronic  performance  of 
the  new  cold  source  was  tested.  Startup  tests  measured 
the  helium  refiigerator  heat  load  as  reactor  power  was 
increased.  At  20  MW,  the  heat  deposited  in  the  source 
was  measured  at  1300  watts  in  excellent  agreement  with 
calculations.   The  cold  neutron  flux  was  monitored  to 
ensure  that  the  thermosiphon  functioned  smoothly  and  that 
the  moderator  volume  remained  full  of  liquid  hydrogen  as 
the  evaporation  rate  increased  to  over  3  g/s  at  full  power. 
Further  tests  established  the  normal  operating  parameters 
for  the  cryostat  and  verified  the  protection  systems  under 
various  abnormal  conditions. 

The  performance  of  the  new  source  was 
benchmarked  against  Unit  1  by  duplicating  flux  measure- 
ments of  the  intensity  as  a  function  of  neutron  wavelength 
at  several  Guide  Hall  instruments.  Figure  2  is  a  plot  of  the 
cold  source  gain,  defined  as  the  ratio  of  intensities  between 
Unit  2  and  Unit  1 ,  as  a  function  of  neutron  wavelength. 
The  figure  includes  both  the  gain  factors  measured  at  three 
spectrometers  and  as  obtained  from  the  MCNP  Monte 
Carlo  simulafions  [2].  While  the  agreement  is  excellent,  the 
measured  gains  were  actually  somewhat  greater  than 
predicted,  especially  at  long  wavelengths.  This  additional 
gain  is  likely  due  to  the  new  guide  sections  replaced  during 
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Fig.  2.  Gains  in  cold  neutron  flux  measured  at  the  NG-7  and  NG-1  small 
angle  neutron  scattering  (SANS)  instruments,  and  at  the  disk  chopper 
spectrometer  (DCS)  on  NG-4,  compared  to  the  predicted  gains,  based  on 
Monte  Carlo  simulations  using  MCNP. 

the  cold  source  installation,  an  effect  that  is  not  included  in 
the  MCNP  models.  The  new  cryostat  design  also  substan- 
tially reduces  the  number  of  fast  neutrons  from  the  cold 
source  resulting  in  improved  instrument  signal-to-noise 
levels. 

The  advanced  liquid  hydrogen  cold  source  has 
successfiilly  met  all  of  its  design  goals.  Its  success 
demonstrates  how  the  effective  use  of  detailed  computer 
models  to  optimize  a  design  can  lead  to  substantial  perfor- 
mance improvements.  However,  it  also  demonstrates  that 
the  complex  geometry  that  results  from  such  an  optimiza- 
tion requires  more  complex  stress  modeling,  new  fabrica- 
tion technologies  and  even  more  stringent  quality  assurance 
testing.   Most  importantly,  the  success  of  the  Unit  2  cold 
source  design  significantly  enhances  the  measurement 
capability  of  the  cold  neufron  scattering  instrumentation  at 
NIST  with  intensity  gains  of  nearly  a  factor  of  two  in  the 
important  long  wavelength  region  of  the  spectrum. 
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Based  on  Polarized  ^He 
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'  pin-polarized  neutrons  are  useful 
'probes  of  magnetic  matter  since 
both  the  magnetic  scattering  of  neutrons  by 
unpaired  electrons  and  the  scattering  of  neutrons 
by  nuclei  of  non-zero  spin  can  be  strong  functions 
of  the  neutron  spin  state  [1].  Unfortunately,  the 
limitations  of  current  polarizing  and  analyzing  devices 
have  substantially  restricted  the  application  of  this 
powerful  technique.  ^He  spin-filters  have  the  potential  to 
yield  broadband  neutron  polarizers  and  analyzers  that  can 
be  used  for  cold,  thermal,  and  epithermal  neutrons.  Such 
spin  filters  could  make  new  classes  of  neutron  scattering 
experiments  possible. 

For  certain  applications  there  are  presently  no  suitable 
devices.  For  example,  polarization  analysis  of  scattered 
beams  having  wide  angular  divergence  is  impractical  with 
reflection-based  devices  such  as  supermirrors  because  of 
their  inadequate  angular  acceptance  and  small-angle 
scattering.  At  present,  we  are  developing  ^He-based 
polarization  analyzers  for  diffuse  reflectometry  and  small- 
angle  scattering  (SANS)  with  neutrons,  both  of  which 
require  analysis  of  a  large  divergence  beam.  ^He-based 
neutron  spin  filters  would  also  be  advantageous  on  crystal 
spectrometers  since  any  monochromator  or  analyzer 
crystal  can  be  used,  rather  than  just  Heusler  alloy,  greatly 
increasing  the  range  and  flexibility  of  the  instrumentation. 
Furthermore,  polarizing  or  analyzing  thermal  and  hot 
neutrons  at  spallation  sources,  where  the  time-of-flight 
method  is  employed,  often  requires  the  broadband  capabil- 
ity of  ^He  spin  filters.  In  addition  to  their  utility  for  neutron 
scattering,  polarized  ^He  spin  filters  are  also  of  interest  for 
nuclear  and  particle  physics  studies  with  neutrons,  such  as 
determination  of  weak  coupling  constants  and  tests  of  the 
Standard  Model  via  accurate  measurements  of  decay 
correlation  coefficients  in  polarized  neutron  beta-decay  [2]. 

Neutron  spin  filters  based  on  polarized  ^He  rely  on 
the  strong  spin-dependence  of  the  neutron  capture  cross 
section  for  ^He.  For  a  sufficient  thickness  of  100  % 
polarized  ^He  gas,  essentially  all  of  the  neutrons  with 
antiparallel  spin  would  be  absorbed,  while  nearly  all  of  the 
neutrons  with  parallel  spin  would  be  transmitted,  resulting 
in  100  %  neutron  polarization  and  50  %  transmission.  In 
Fig.  1  we  show  calculated  values  of  the  neutron  polariza- 


tion    and  neutron  transmission     for  ^He  polarization 
Pjjg  =  60  %,  an  experimentally  achievable  value.  Since 
there  is  a  tradeoff  between  neutron  polarization  and 
transmission,  we  also  show  P^^T^^,  which  is  a  useful 
figure-of-merit  for  many  experiments.  A  transmission 
analyzer  can  be  characterized  either  by  the  transmission 
asymmetry  A  or  the  flipping  ratio  F.  For  T+  and  T_  defined 
to  be  the  transmissions  for  neutrons  with  polarization 
parallel  and  antiparallel  to  the  ^He  polarization,  respectively, 
the  asymmetry  is  given  by  ^  =  {T+  -  T_)I{T+  +  T_)  and  the 
flipping  ratio  by  F=  TJT_  =  (1  +  A)l{\  -  A).  The  flipping 
ratio  F  is  also  shown  in  Fig.  1 .  The  asymmetry  for  a  spin 
filter  used  as  an  analyzer  is  the  same  as  the  neutron 
polarization  produced  when  it  is  used  as  a  polarizer.  For 
the  specific  case  of  cold  neutrons  (A  =  0.5  nm),  a  spin 
filter  used  as  a  polarizer  with  P^e  ^  60  %  and  a  pressure- 
length  product  of  7  kPa-m  would  yield     =  90  %  and 

=  20  %,  or,  when  used  as  an  analyzer,  would  give 
A  =  90  %  and  F=  19. 

We  produce  polarized  ^He  gas  by  two  optical  pump- 
ing methods:  spin-exchange,  which  is  performed  directly 
at  high  pressure  (0.1  MPa  to  0.3  MPa),  and  metastability- 
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Fig.  1.  Neutron  polarization  P„  (solid  thick  line),  transmission  7„  (solid 
thin  line),  a  figure  of  merit  given  by  P^l^  (dashed  line),  and  the  flipping 
ratio  F  (dotted  line)  as  a  function  of  the  product  of  ^He  pressure  in  Pa 
(assuming  a  room  temperature  cell),  cell  length  in  m  and  wavelength  in 
nm.  The  calculations  are  shovi/n  for  ^He  polarization  Pne"    "Z"-  "^^^ 
scale  on  the  right  hand  y-axis  is  used  for  the  flipping  ratio. 
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exchange,  in  which  the  gas  is  polarized  at  low  pressure 
(s;  100  Pa)  and  then  compressed  [3].  The  spin-exchange 
method  is  convenient  and  well  matched  to  continuous 
operation  on  a  beamline.  whereas  the  metastability- 
exchange  method  has  a  higher  polarizing  rate.  At  present, 
it  has  been  more  com  enient  (for  either  method)  to  polarize 
gas  off-line  and  transport  it  to  the  beam  line.  Maintaining 
the  polarization  in  the  absence  of  optical  pumping  requires 
a  homogeneous  magnetic  field  and  specially  prepared  glass 
cells  with  slow  wall  relaxation.  We  have  produced  cells 
with  relaxation  times  as  long  as  one  month,  dominated  by 
the  intrinsic  dipole-dipole  relaxation  in  the  ^He  gas  itself 
[4].  In  the  future,  continuously  operating  spin  filters  can 
be  installed  directh  on  the  neutron  instrument. 

Polarization  analysis  allows  for  separation  of  mag- 
netic from  nuclear  scattering,  and  also  separation  of 
coherent  from  spin-incoherent  scattering.  In  our  first 
demonstration  experiment  for  SANS,  we  used  a  ^He  spin 
filter  to  extract  a  small  component  of  spin-incoherent 
scattering  in  the  presence  of  strong  coherent  scattering 
[5].  We  ha^  e  continued  with  tests  of  separating  magnetic 
from  nuclear  scattering. 

Recently  we  carried  out  successful  tests  of  a  polar- 
ized ^He  spin  filter  on  the  NCNR  NG- 1  reflectometer.  The 
first  test  experiment  on  a  specular  reflection  was  a  careful 
comparison  of  results  obtained  with  a  ^He  analyzer  to  those 
obtained  with  the  current  technique  that  employs  a 
supermirror  anah  zer.  The  test  sample  was  an  epitaxial 
Muq  52Pdo.48/Fe  (001)  bilayer,  for  which  chemical  ordering 
and  magnetic  exchange-bias  have  recently  been  reported 
[6].  By  var\'ing  the  magnetic  field  on  this  sample,  it  was 
possible  to  test  the  ^He  analyzer  imder  conditions  of  both 
significant  and  negligible  spin-flip  scattering.  We  employed 
a  compact,  magnetically  shielded  solenoid  that  was 
interchanged  with  the  analyzing  supermirror  on  the  NGl 
mstrument.  The  solenoid  adequately  shielded  the  stray 
fields  from  the  0.6  T  magnet  and  guide  field  magnets.  The 
results  obtained  with  the  ^He  analyzer  were  identical  to 
those  obtained  with  the  supermirror  analyzer. 

For  efficient  studies  of  diffuse  scattering  on  a 
reflectometer,  the  ^He  analyzer  will  be  combined  with  a 
position  sensitive  detector  (PSD).  Two  issues  not  tested 
in  the  specular  experiment  are  relevant  in  this  case: 
(1)  possible  depolarization  of  neutrons  that  follow  trajecto- 
ries off  the  axis  of  the  ^He  analyzer,  and  (2)  possible 
background  from  the  ^He  cell.  We  have  established  that 
most  of  the  area  of  apertures  in  the  magnetic  shield  is 
usable,  with  depolarization  only  occurring  for  a  small  range 
of  extreme  trajectories.  We  have  also  determined  that  small 


angle  scattering  from  the  cell  is  negligible  and  therefore 
should  not  pose  an  issue  for  the  low  neutron  fluxes 
expected  in  diffuse  scattering  experiments.  This  series  of 
tests  established  the  suitability  of  polarized  ^He  spin  filters 
for  diffuse  reflectometry  experiments,  which  we  will 
pursue  in  the  near  future. 

We  are  making  continual  improvements  in  the 
polarization  and  relaxation  time  of  our  ^He  cells.  In  the 
recent  reflectometry  experiments,  the  initial  -''He  polariza- 
tion was  57  %,  resulting  in  a  neutron  flipping  ratio  of  3 1 
and  a  transmission  for  the  desired  spin  state  of  24  %.  The 
relaxation  time  of  the  polarized  gas  was  1 5  days  on  the 
beam  line,  which  was  somewhat  reduced  from  the  intrinsic 
value  of  24  days  for  the  cell  due  to  a  small  fractional 
gradient  of  3  x  lO-'^  cm"i  in  the  magnetic  field  of  the 
shielded  solenoid.  Our  aim  is  to  increase  the  ^He  polariza- 
tion to  70  %  and  obtain  relaxation  times  at  the  dipole-dipole 
limit. 
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Development  of  an  Ion-Implanted 
Phosphorus  in  Silicon  SRM 
for  the  Semiconductor  Industry 


'  ilicon  is  the  foundation  material 
of  the  electronic  age,  lying  at  the 
base  of  the  200  billion  dollar  world-wide 
semiconductor  industry.  Silicon's  properties  are 
manipulated  by  precise  implantation  and  distribution 
of  other  atoms  (dopants),  which  makes  the  accurate 
measurement  of  impurity  distributions  imperative. 
SEMATECH  (a  consortium  of  semiconductor  manufac- 
turers) recently  listed  Standard  Reference  Material  (SRM) 
implant  of  phosphorus  in  silicon  as  a  high  priority  industrial 
need. 

The  U.S.  semiconductor  industry  relies  heavily  on 
secondary  ion  mass  spectrometry  (SIMS)  for  characteriza- 
tion of  the  depth  distribution  of  dopants  in  silicon.  To 
achieve  high  accuracy  in  the  concentration  determination 
by  SIMS,  standards  of  known  dopant  concentration, 
conveniently  provided  by  ion  implants  of  certified  dose,  are 
required.  Standard  Reference  Materials  of  boron  and 
arsenic  implants  in  silicon  (SRMs  2137  and  2134  [1])  have 
already  been  developed  by  NIST  as  SIMS  calibration 
standards. 

The  SIMS  community  in  the  United  States  undertook 
a  round-robin  study  to  calibrate  the  implanted  dose  of 
phosphorus  in  silicon  by  consensus.  Dose  determinations 
among  laboratories  varied  by  nearly  a  factor  of  2,  reflect- 
ing primarily  the  errors  of  the  respective  in-house  stan- 
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Fig.  1.  Round-robin  measurements  of  phosphorus  in  silicon  by 
secondary  ion  mass  spectrometry  in  16  U.S.  industrial  laboratories, 
demonstrating  the  need  for  a  common  phosphorus  reference  material  to 
improve  inter-laboratory  reproducibility. 


dards  (Fig.  1 ).  These  results  demonstrate  the  need  for  a 
common  phosphorus  reference  material  to  improve  inter- 
laboratory  reproducibility. 

In  pursuit  of  a  phosphorus  in  silicon  SRM,  a  radio- 
chemical neutron  activation  analysis  (RNAA)  procedure 
has  been  developed,  critically  evaluated,  and  shown  to  have 
the  necessary  sensitivity,  chemical  specificity,  matrix 
independence,  and  precision  to  certify  phosphorus  at  ion 
implantation  levels  in  silicon  [2,  3].  During  sample  irradia- 
tion, 3ip  undergoes  neutron  capture  to  form  ^^P,  which  is 
then  separated  from  the  matrix  and  measured  by  beta 
counting.  The  procedure  is  used  here  to  assign  a  value  to 
the  quantity  of  phosphorus  in  SRM  2133  (phosphorus 
implanted  silicon). 

SRM  2133  was  prepared  by  phosphorus  implantation 
of  a  200  cm  diameter  silicon  wafer.  The  wafer  was  cut 
into  1  cm-  pieces  by  means  of  a  wafer  saw,  and  twelve 
pieces  were  chosen  for  analysis.  The  area  of  each  piece 
was  determined  by  measuring  its  dimensions  using  a  digital 
micrometer.  During  irradiation,  ^^p  may  be  produced  via 
neutron  capture  by  the  implanted  phosphorus,  by  bulk 
phosphorus  impurity  in  the  wafer,  and  by  silicon.  The 
latter  two  sources  lead  to  measurement  error,  which  must 
be  corrected.  These  corrections  were  made  by  measuring 
silicon  blanks,  which  were  prepared  by  cutting  a  non- 
implanted  silicon  wafer  into  ~  1  cm^  pieces.  Standards 
were  prepared  by  deposition  of  standardized  solutions  of 
phosphorus  on  aluminum  foils. 

Samples,  standards,  and  blanks  were  irradiated  for  3 
hours  in  irradiation  tube  RTl  of  the  NIST  neutron  source 
at  a  neutron  fluence  rate  of  1.05  x  10'^  cm-^-s"'.  Targets 
were  then  processed  using  the  method  developed  previ- 
ously [2,  3].  Silicon  samples  and  blanks  were  mixed  with 
a  few  milligrams  of  non-radioactive  phosphorus  (carrier) 
and  dissolved  in  a  mixture  of  nitric  and  hydrofluoric  acids. 
Phosphorus  was  separated  from  the  matrix  first  by 
precipitation  as  ammonium  phosphomolybdate,  then  as 
magnesium  ammonium  phosphate.  The  latter  was  col- 
lected by  filtration,  washed  with  dilute  ammonium  hydrox- 
ide and  ethanol,  air  dried,  weighed,  and  packaged  for 
beta-ray  counting.  The  yield,  (fraction  of  recovered 
phosphorus  carrier),  was  determined  gravimetrically, 
assuming  a  composition  of  MgNH4P04-6H20.  Standards 
were  processed  using  a  similar  procedure. 


10 


Instrumentation  Development  and  Applications 


1 


Counting  of  ^-P  (ti/2  =  14.28  d)  was  performed  using 
a  beta  proportional  counter.  Purity  of  the  ^-F  from  other 
radionuclides  was  ascertained  by  gamma-ray  spectros- 
copy, and  by  monitoring  beta  count  rates  over  a  period  of 
four  half  lives.  No  discernible  contamination  from  other 
radionuclides  was  noted.  ^-P  count  rates  were  adjusted  for 
decay  time,  blank  correction,  and  for  differences  in  beta 
self-absorption. 

NIST  currently  certifies  elemental  concentrations  in 
SRMs  using  one  of  three  modes:  (1)  a  primary  method  at 
NIST  with  confirmation  by  other  method(s);  (2)  two 
independent  critically  evaluated  methods  at  NIST,  and 
(3)  one  method  at  NIST  and  different  methods  by  outside 
collaborating  laboratories.  Certification  using  a  primary 
method  is  only  possible  when  all  potentially  significant 
sources  of  uncertainty  have  been  evaluated  explicitly  for 
the  application  of  the  method  and  the  matrix  under  investi- 
gation. Therefore,  in  order  to  certify  the  phosphorus 
concenfration  in  SRM  2133  using  RNAA  as  a  primary 
method  it  was  necessary  to  evaluate  all  significant  sources 
of  uncertainty  explicitly.  For  this  set  of  measurements,  we 
considered  sources  of  relative  uncertainty  greater  than 
0.01  %  to  be  significant.  The  results  of  a  complete 
evaluation  of  all  sources  of  uncertainty  are  hsted  in  Table  1 . 
This  evaluation  yielded  an  expanded  relative  uncertainty  of 
1.27  %  (as  defined  by  ISO  and  NIST),  and  gives  an 
approximate  level  of  confidence  of  95  %.  The  measured 
phosphorus  concentration  is  (9.55  ±  0.12)  x  lO'^  atoms- 
cm-. 

In  conclusion,  we  have  successfully  applied  RNAA 
as  a  primary  method  for  the  certification  of  this  new  SRM. 
The  observed  relative  expanded  uncertainty  of  1 .27  %  is 
smaller  than  the  3  %  value  needed  by  the  semiconductor 
industry.  This  new  SRM  should  greatly  enhance  the  U.S. 
semiconductor  industry's  ability  to  achieve  accurate  and 
reproducible  analytical  results  for  this  key  dopant  in  silicon. 


Table  1.  Individual  Uncertainty  Components  for  Determination 
of  Phosphorus  in  SRM  2133. 


Source  of  Uncertainty 

Uncertainty 

(la)  % 

Measurement  Replication  (a/V n) 

0.31 

Measurement  Replication  (Standards)  (a/V n) 

0.47 

Measurement  of  Sample  Area 

0.012 

Blank  Correction 

0.022 

Blank  Correction  (Standards) 

0.019 

Volumetric  Calibrations  (Standards) 

0.12 

Concentration  of  Standard  Solution 

0.15 

Mass  Standard  Solution  on  Foil 

0.037 

Radiochemical  Impurity 

0.13 

Beta  Self-Absorption 

0.07 

Beta  Self-Absorption  (Standards) 

0.08 

Carrier  Yield  Determination 

0.10 

Self-Shielding 

0.024 

Irradiation  Geometry 

0.09 

Combined  Uncertainty 

0.63 

Coverage  Factor 

2.0 

Expanded  Uncertainty 

1.27 
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|(^  New  Thermal  Neutron  Prompt 

Gamma-Ray  Activation  Analysis  Facility 


'  cientists  from  the  Analytical 
'Chemistry  Division  and  the  US  Food 
and  Drug  Administration  have  maintained  a 
thermal  neutron  prompt  gamma-ray  activation 
analysis  (PGAA)  facility  at  the  NCNR  for  the  past 
25  years.  The  facility  is  used  routinely  for  analysis  of 
major  and  minor  element  compositions  of  a  variety  of 
biological,  environmental  and  industrial  materials.  Recent 
applications  include  determination  of:  several  elements  in 
Standard  Reference  Material  (SRM)  1575a  Pine  Needles; 
cadmium  in  SRM  2702  Inorganic  Sediment;  nitrogen  and 
phosphorus  in  cattail  samples  in  collaboration  with  Florida 
A&M  University;  boron  in  sapphire  in  collaboration  with 
Southern  University;  and  nitrogen,  sodium,  phosphorus, 
sulftir,  chlorine,  potassium,  calcium,  and  cadmium  in  food, 
dietary  supplements,  and  other  food-related  matrices. 

A  new  thermal  neutron  PGAA  instrument  became 
operational  in  2001,  replacing  the  original  thermal  neutron 
PGAA  instrument  that  was  built  in  the  late  1970s.  The  new 
instrument  was  designed  and  built  to  improve  analytical 
sensitivities  and  limits  of  detection,  and  decrease  the  levels 
of  radiation  associated  with  use  of  the  instalment.  This 
was  achieved  by  filtering  the  neutron  beam  to  decrease  the 
fast  neutron  and  low  energy  gamma-ray  components  of 
the  beam,  designing  the  new  facility  components  and 
support  structure  to  minimize  neutron  capture  gamma- 
rays,  and  installing  an  improved  detection  system.  The 

new  thermal  neutron 
PGAA  facility  is 
shown  in  Fig.  1. 

A  cylindrical 
sapphire  beam  filter 
5.3  cm  long  and 
4.3  cm  in  diameter 
was  incorporated  into 
the  shutter  assembly. 
This  modification 
reduced  the  number 
of  fast  neutrons  by  a 
factor  of  5,  decreased 
the  thermal  neutron 

flux  by  approximately 
Fig  1.  Newtherr^al  neutron  PGAA  y  j 

facility  at  VT-5. 


reduced  low-energy  gamma-ray  background.  With  the 
reduction  in  fast  neutrons,  the  amount  of  hydrogenous 
shielding  required  to  slow  and  capture  the  fast  neutrons 
was  reduced  as  compared  with  that  needed  in  the  old 
facility.  As  a  result,  the  background  count-rate  for  the 
2223.3  keV  capture  gamma-ray  of  hydrogen  has  been 
reduced  from  about  1  count  per  second  (cps)  or  the 
equivalent  of  about  0.5  mg,  of  hydrogen,  to  0.05  cps  or 
the  equivalent  of  25  \ig  of  hydrogen. 

A  summary  of  the  design  changes  for  each  compo- 
nent of  the  facility  and  their  effects  is  shown  in  Table  1 . 
The  original  beam  tube  and  sample  chamber  were  fabri- 
cated from  concentric  Plexiglas  tubes,  filled  with  a  mixture 
of  paraffin  and  boron  carbide.  The  new  sample  chamber 
and  beam  tube  were  fabricated  from  aluminum  and  lined 
with  a  2  mm  thick  flexible  lithiated  polymer  (lithoflex). 
The  beam  tube  is  kept  under  vacuum  to  reduce  back- 
ground gamma-rays  arising  from  neutrons  scattered  by  air 
and  captured  in  the  surrounding  materials.  The  new 
sample  chamber  can  be  evacuated  to  fiarther  reduce 
background  associated  with  neutron  scatter  and  capture  in 
air.  As  a  result  of  these  design  changes,  the  background 
count-rates  from  hydrogen,  nitrogen,  and  boron  have  been 
reduced  by  a  factor  of  1 0  or  more  and  carbon  has  been 
eliminated  from  the  background.  The  count  rate  for 
aluminum  is  about  the  same  in  the  new  system,  despite  the 
large  amounts  of  aluminum  in  the  beam  tube,  beam  stop, 
sample  chamber  and  support  structure. 

The  beam  stop  of  the  old  facility  was  a  separate  unit. 
It  consisted  of  a  steel  box  containing  a  lithiated  polycar- 
bonate cube  (30  cm  on  a  side)  surrounded  by  lead  bricks. 
This  unit  was  lifted  into  place  on  top  of  the  beam  tube  each 
time  the  system  was  assembled.  The  new  beam  stop 
consists  of  an  aluminum  box  welded  onto  the  aluminum 
support  structure,  and  is  filled  with  a  cube  of  borated 
polyethylene  (20  cm  on  a  side)  surrounded  by  lead.  The 
elimination  of  steel  from  the  beam  stop  has  greatly  reduced 
the  background  count-rate  from  iron. 

A  new  germanium  detector  with  improved  resolution 
and  peak-to-Compton  ratio  was  installed.  This  detector  has 
a  relative  efficiency  of  40  %  (relative  to  a  7.6  cm  cube 
of  sodium  iodide)  and  2.0  keV  (peak  fiill-width-half- 
maximum)  resolution  for  the  1332.5  keV  gamma-ray  line 
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Table  1 .  Comparison  of  original  and  new  thermal  neutron  PGAA  systems 


Facility 
Component 

Old  PGAA  System 

New  System 

Effect 

Beam  Stop 

Steel  cube  (61  cm  on  a  side)  Li2C03  plug,  cube  of 
Li2C03  in  resin  surrounded  by  lead  bricks 

Aluminum  cube  (43  cm  on  a  side)  ^LiF  polymer  plug,  berated 
polyethylene  surrounded  by  lead  bricks 

Eliminate  Fe 
background  y-rays 
from  beam  stop 

Beam  tube 

Plexiglas  concentric  tubes  with  a  mixture  of  boron 
carbide  in  paraffin,  and  Li2C03  in  paraffin  between 
them;  Teflon  film  end  windows;  airfilled  beam  tube 

Al  tubing  lined  with  ^LiF  polymer  capped  at  each  end  with  Al 
and  center  windows  of  a  Mg  alloy  Upper  and  lower  sections 
evacuated,  lithiated  polyethylene  shielding  on  outer  two  sides 

Reduce  H  and 
eliminate  N 
background  y-rays 

Sample 
Chamber 

Walls  fabricated  from  B4C  in  resin  with  Teflon  windows 
covered  with  fused  enriched  Li2C03  plugs  in  front  of 
and  across  from  the  detector  window 

Welde(;|  Al  plate  covered  with  ^LiF  polymer  shielding;  Mg 
alloy  thin  windows  in  front  of  and  across  from  the  detector 
covered  with  the  same  Li2C03  plugs;  chamber  may  be 
evacuated 

Reduce  B  and 
eliminate  N 
background  y-rays 

Detection 
system 

Ge  detector  (27  %  relative  efficiency)  in  a  Nal 
Compton  suppressor 

Ge  detector  (40  %  relative  efficiency)  in  a  BGO  Compton 
suppressor 

Eliminate  Na 
background  y-rays 

Support 
Structure 

Al  tripod  frame  with  a  cage  to  hold  the  beam  stop 

Al  tripod  frame.  The  beam  tube  and  beam  stop  are  attached 
to  the  frame  and  the  detection  system  placed  on  rails 

Improve  structure 
stability 

Component 
Set-Up 

Four  separate  pieces  (beam  tube,  Al  tripod,  beam 
stop,  detection  system)  that  are  assembled  individually 

Assembled  as  one  unit  with  a  removable  Ge  detector 

Improve  safety  and 
ease  of  system  set-up 

10^ 

10^ 

tn 

10^ 

Coi 

10^ 

10^ 

10^ 

10° 

Background  Spectrum,  Original  PGAA  System 

T 


4  6 
Energy  (MeV) 

Background  Spectrum,  New  PGAA  System 

T 


Energy  (MeV) 

Fig.  2.  Comparison  of  tfie  background  spectra  collected  with  the  old 
PGAA  system  (top)  and  the  new  system  (bottom). 


from  6OC0.  The  original  sodium  iodide  Compton  suppres- 
sion system  was  replaced  with  a  bismuth  germanate 
suppressor.  This  modification  eliminated  sodium  gamma- 
ray  background. 

The  goal  of  reducing  background  radiation  and 
improving  analytical  sensitivity  and  detection  limits  has 
been  achieved.  A  comparison  of  spectra  from  background 
counts  from  the  old  and  the  new  PGAA  systems  is  shown 
in  Fig.  2.  The  reduction  in  the  background  count-rates  in 
the  low  energy  region  is  evident.  This  reduction  will 
greatly  improve  the  detection  limits  for  elements  with 
neutron  capture  gamma-rays  in  the  low  energy  region. 
With  the  sample  chamber  evacuated,  nitrogen  capture 
gamma-rays  have  been  eliminated  from  the  background. 
Elimination  of  the  sodium  background  has  improved  the 
detection  limit  for  this  element.  Improvements  in  shield- 
ing have  resulted  in  a  decrease  in  the  background  count- 
rates  of  hydrogen,  boron,  carbon,  iron,  and  germanium, 
improving  detection  limits  for  those  elements  as  well. 
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Performance  Enhancements  of 
the  CHRNS  30-m  SANS  Instrument 


T 

■  t 


'wo  changes  made  to  the  optics  of 
the  CHRNS  30-m  SANS  instrument 
during  the  installation  of  the  NCNR's  new  cold 
source  have  increased  the  flux  at  the  sample  and 
the  \ow-Q  resolution  of  the  instrument. 

The  increase  in  flux,  over  and  above  that 
provided  by  the  new  cold  source,  comes  from  replac- 
ing the  instrument's  cooled  bismuth-beryllium  filter  with 
an  "optical  filter"  of  the  type  shown  schematically  in 
Fig.  1.  The  optical  filter  replaces  10  m  of  straight  neutron 
guide  and  its  40  cm  long  Bi-Be  filter  with  an  inclined  guide 
section  with  a  high  critical  angle  supermirror  reflective 
coating.  The  optical  filter  channels  cold  neutrons,  with 
wavelengths  as  short  as  3  A,  out  of  the  line-of-sight  of  fast 
neutrons  and  gamma  rays  fi"om  the  source.  The  cold 
neutrons  delivered  to  the  SANS  instrument  undergo  an 
even  number  of  reflections  in  the  optical  filter  to  emerge  in 
the  horizontal  direction,  but  displaced  14  cm  vertically.  To 
accommodate  the  beam  displacement,  the  entire  SANS 
instrument  had  to  be  raised;  a  non-trivial  task  that  required 
detailed  engineering  analysis,  careful  planning  and  skilled 
execution. 

The  design  of  the  optical  filter  was  refined  with  the 
aid  of  detailed  Monte  Carlo  calculations  of  the  angular  and 
spectral  distribution  of  neutrons  transmitted  by  the  optical 
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Fig.  1 .  Schematic  elevation  of  the  NG-3  optical  filter  (note  vertical  scale 
is  exaggerated  for  clarity).  The  inclined  section  shown  in  red  is  coated 
with  supermirror  with  critical  angle  approximately  3.2  times  that  of 
natural  nickel.  The  removable  guide  sections  are  in  the  pre-sample  flight 
path  of  the  NG-3  SANS  instrument. 


filter.  Figure  2  shows  the  calculated  and  measured  gain  in 
flux  at  the  sample,  due  to  the  optical  filter  alone,  for  the 
supermirror  reflectivity  model  shown  in  the  inset.  These 
results  demonstrate  that  the  optical  filter  transmission  is 
about  the  same  as  the  previous  crystal  filter  for  wave- 
lengths around  5  A,  but  becomes  substantially  better  at 
longer  wavelengths  where  absorption  in  the  crystal  filter 
becomes  significant. 

The  optics  for  the  CHRNS  SANS  instrument  were 
further  improved  by  installing  a  system  of  retracting  lenses 
and  prisms  near  the  sample  position  to  focus  17  A  neu- 
trons onto  the  detector  at  its  maximum  distance,  13  m, 
from  the  sample.  A  similar  lens  system,  consisting  of  28 
biconcave  single  crystals  of  MgF2  for  focusing  8  A 
neutrons,  has  been  in  use  on  the  NCNR's  other  30-m 
SANS  instrument  on  guide  NG-7  for  nearly  two  years  [1]. 
The  refracting  power  of  the  lenses  increases  with  the 
wavelength  squared,  but  so  does  the  distance  the  neutron 
falls  between  the  sample  and  detector  due  to  gravity.  The 
vertical  spreading  of  the  focus  for  a  beam  with  a  wave- 
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Fig.  2.  Transmission  gains  of  the  NG-3  optical  filter.  The  black  circles 
are  measured  with  no  removable  guides  in  the  beam.  The  red  squares 
are  measured  with  8  removable  guides  in  the  beam.  The  violet  circles 
are  the  original  Monte  Carlo  simulated  gain  predictions  (March  1999) 
using  the  supermirror  reflectivity  model  shown  in  the  inset.  The  blue 
circles  are  simulated  gains  (July  2002)  using  the  refined  reflectivity 
model  shown  in  the  inset.  These  "hindsight"  simulations  predict  the 
measured  gains  well  for  /W  =  3.2  supermirror  with  R(Q  =  0)  =  0.930  and 
RMS  surface  roughness  equal  to  10  A  with  R(Q  =  0.069  A-^)  =  0.55. 
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length  spread  typical  of  a  SANS  instrument,  about  1 0  %  to 
1 5  %  F  WHM,  restricts  the  utility  of  the  lenses  alone  to 
wavelengths  less  than  10  A.  The  new  development 
implemented  on  the  CHRNS  SANS  instrument  is  to  follow 
the  lenses  with  prisms  that  refract  in  the  vertical  direction 
to  counteract  the  effect  of  gravity. 

Figure  3  depicts  the  arrangement  of  lenses  and 
prisms  now  installed  on  the  CHRNS  SANS  instrument. 
Seven  MgFT  biconcave  lenses  focus  1 7  A  neutrons  at  the 
detector,  13  m  from  the  sample.  For  this  distance  a  single 
prism  with  an  apex  angle  of  161°  would  cancel  the  beam 
spreading  due  to  gravity  for  all  wavelengths  [2].  Such  a 
prism,  however,  would  have  to  be  250  mm  long  at  its  base 
to  intercept  the  full  beam  height  transmitted  by  the  lenses 

2  cm).  A  more  practical  scheme,  as  shown  in  Fig.  3,  is 
to  stack  two  sets  of  prisms,  with  each  prism  30  mm  long 
and  5  mm  high  with  an  apex  angle  of  143",  to  give  the 
same  anti-gravity  effect.  In  this  scheme,  the  bottoms  of 
each  prism  are  coated  with  Gd203  to  eliminate  surface 
reflections. 

The  effectiveness  of  the  lens/prism  combination  can 
be  seen  in  Fig.  4  that  compares  circularly  averaged  SANS 
data  from  voids  in  an  irradiated  crystal  of  aluminum 
obtained  with  pinhole  collimation  and  with  the  focusing 
optics.  The  data  for  the  lens/prism  system  extend  to  a 
minimum  scattering  vector  Q  =  0.00045  A-',  nearly  a 
factor  of  two  lower  than  the  pinhole  collimation.  Further- 
more, the  scattered  intensity  is  3  times  higher  with  the 
lens/prism  system  because  a  larger  area  of  the  sample  can 
be  illuminated  without  degrading  the  ^-resolution.  With 
pinhole  collimation,  intensity  at  the  detector  is  roughly 
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Fig.  4.  SANS  from  voids  in  an  aluminum  crystal  measured  with  and 
without  the  new  focusing  and  gravity  cancellation  optics  installed  on  the 
CHRNS  SANS  instrument.  The  optics  halve  the  minimum  accessible 
Q-value  while  increasing  the  scattered  intensity  roughly  threefold. 

proportional  to  Q^mm^  hence  the  lens/prism  system  repre- 
sents an  overall  intensity  gain  of  48  compared  with  simply 
reducing  pinhole  aperture  size  to  achieve  the  same  mini- 
mum Q. 

The  new  cold  source  and  improvements  in  optics  are 
enabling  microstructural  studies  that  link  nanoscale  with 
microscale  features  in,  for  example,  polymer-clay 
nanocomposites,  gels,  and  fluxoid  lattices  in  superconduc- 
tors. 
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Correction  of  Gravity  Effect  with  Prisms 
Source  Lens    Prism  Sample  Detector 


Fig.  3.  Schematic  diagram  of  the  configuration  of  lenses  and  prisms 
used  on  the  30-m  CHRNS  SANS  instrument.  The  seven  biconcave 
lenses  focus  17  A  neutrons  at  the  detector,  13  m  from  the  sample,  and 
the  double  stack  of  prisms  refracts  the  beam  vertically  to  cancel  the 
effect  of  gravity. 
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Probing  Nanoscale  Disorder 
with  Neutron  Tunneling  Spectroscopy 


■he  presence  of  structural  disorder 
can  have  a  dramatic  effect  on 
macroscopic  materials  properties.  At  the 
atomic  scale,  structural  disorder  is  typically 
described  in  terms  of  distributions  of  geometrical 
quantities  such  as  bond  lengths  and  bond  angles  that 
are  determined  from  diffraction  measurements.  Alter- 
natively, disorder  can  be  quantified  by  determining  the 
distribution  of  potential  energies  experienced  by  the  atomic 
scale  constituents  of  the  system.  The  distribution  is  sharp 
for  a  highly  ordered  (crystalline)  system,  while  it  is  broad 
for  a  disordered  system.  Rotational  tunneling  spectros- 
copy is  a  useful  tool  for  characterizing  the  local  molecular 
environment  due  to  its  exquisite  sensitivity  to  the  local 
potential  [1].  This  approach  has  been  used  to  determine 
the  distribution  of  rotational  potentials  felt  by  methyl  iodide 
molecules  adsorbed  in  a  series  of  mesoporous  silica 
glasses.  The  degree  of  confinement-induced  disorder  in 
the  methyl  iodide  can  then  be  directly  related  to  the  pore 
sizes  in  the  glass. 

The  rotational  potential  felt  by  the  methyl  (CH3) 
groups  in  methyl  iodide  (CH3I)  arises  solely  from  interac- 
tions with  the  local  environment  including  the  neighboring 
molecules.  The  energies  of  the  transitions  between  rota- 
tional states  can  be  found  by  solving  the  Schrodinger 
equation  for  the  model  potential  by  V3(l-cos(39))/2  in  the 
angular  coordinate  about  the  C-I  axis,  9.  In  the  limit  of  an 
infinitely  high  barrier,  V3,  the  motion  of  the  CH3  group  is 
that  of  a  torsional  oscillator.  However,  if  the  barrier  is  not 
too  large,  the  wavefunctions  will  have  significant  overlap 
and  quantum  mechanical  tunneling  between  the  potential 
minima  can  occur.  This  results  in  a  splitting  of  the 
rotational  ground  state.  The  energy  of  this  splitting,  which 
can  be  directly  measured  using  cold  neutron  spectroscopy, 
has  a  nearly  exponential  dependence  upon  the  barrier 
height.  This  strong  dependence  on  V3  is  the  reason  why 
turmeling  spectroscopy  is  so  sensitive  to  the  local  molecu- 
lar environment  and,  in  particular,  the  quenched  disorder  of 
adsorbed  molecular  solids. 

The  porous  hosts  used  in  this  investigation  have  a 
range  of  pore  size  distributions  with  differing  nominal  pore 
diameters  and  widths.  Nominal  pore  diameters,  <d>, 
determined  via  nitrogen  adsorption/desorption  isotherms 
range  from  2.5  nm  to  14.4  nm.  In  all  cases  but  one 


(xerogel  with  <d>  =  10.0  nm)  the  porous  glasses  are  disks. 
The  xerogel  is  a  fine  powder  with  a  nominal  particle  size  of 
1 1  \xm.  Liquid  methyl  iodide  (CH3I)  was  condensed  into 
the  porous  hosts  until  95  %  of  the  open  pore  volume  was 
filled  and  then  cooled  to  5  K.  Spectra  were  collected  using 
the  NIST  High-Flux  Backscattering  Spectrometer  with  an 
energy  resolution  of  0.79  fj.eV. 

The  neutron  scattering  spectra  display  broad  asym- 
metric lines  indicative  of  a  distribution  of  potential  barrier 
heights.  These  distributions  can  be  extracted  firom  the 
neutron  scattering  data  using  a  simple  model  in  which  a 
sum  of  Gaussian  distributions  of  potential  barriers  (which 
yield  an  asymmetric  tunneling  line-shape)  is  used  to  fit  the 
spectra  [2].  A  sum  of  two  potential  distributions  is  found 
to  be  sufficient  to  describe  the  data.  One  of  these  distribu- 
tions is  quite  broad,  indicating  substantial  disorder  of  the 
methyl  iodide  solid,  while  the  other  is  much  narrower 
suggesting  a  relatively  high  degree  of  order. 

The  relative  weights  of  the  two  components  show  a 
pronounced  dependence  on  <d>.  Figure  1  shows  that  the 
fi-action  of  the  total  integrated  intensity  that  appears  in  the 
highly  disordered  component  of  the  tunneling  spectrum 
decreases  monotonically  as  <d>  increases.  This  behavior 
can  be  described  with  a  simple  model  in  which  the  rela- 
tively well-ordered  component  of  the  scattering  is  attrib- 
uted to  CH3I  molecules  near  the  center  of  the  pore  while 


10  15 
<cl>  (nm) 


Fig.  1.  Fraction  of  the  total  integrated  intensity  that  appears  in  the 
disordered  component  of  the  tunneling  spectrum  as  a  function  of  the 
mean  pore  diameter.  The  solid  line  is  a  fit  to  the  model  discussed  in  the 
text.  The  dashed  line  is  drawn  at  twice  the  disordered  layer  thickness. 
Inset  shows  the  geometrical  picture  used  in  the  model. 
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the  highly  disordered  component  is  associated  with  a  ring 
of  thickness  t  of  disordered  molecules  near  the  pore  walls. 
Treating  the  pores  as  circular  cy  linders,  the  fraction  of 
disordered  molecules  is  just  the  ratio  of  the  ring  to  circular 
areas:  /^ordered  =  '^t(<^  -  t)/<ci>^. 

The  solid  line  in  Fig.  1  compares  this  model  to  the 
data  for  t  =  1.41(4)  nm,  which  corresponds  to  about  three 
molecular  la\  ers.  It  is  clear  that  this  model  provides  an 
excellent  fit  to  the  data. 

While  the  relative  weights  of  the  two  components 
depend  on  the  pore  size,  the  mean  values,  <V3>,  and 
widths.        of  the  two  distributions  of  barrier  heights  are 
found  to  be  roughly  the  same  for  each  of  the  hosts.  This 
indicates  that  the  structure  of  the  confined  methyl  iodide 
both  near  the  surface  and  in  the  center  of  the  pores  does 
not  depend  strongly  on  either  <d>  or  the  width  of  the  pore 
size  distribution. 

CH3I  molecules  possess  a  relatively  large  electric 
dipole  moment.  This  results  in  the  molecules  being  anti- 
aligned  in  the  bulk  cry  stal  structure  [3].  Silica  surfaces 
often  \ia\e  h\  drox\  1  (OH)  groups  which  would  be  ex- 
pected to  interact  rather  strongK  with  the  dipole  moment 
of  the  molecules  and  lead  to  a  substantial  surface-molecu- 
lar interaction.  The  surface-molecule  and  molecule- 
molecule  interactions  can  be  of  similar  strength  leading  to 
structural  frustration  in  the  confined  solid.  Thus,  one 
would  expect  that  the  disordered  fraction  and  even  the 
structure  of  the  tvvo  components  to  be  strongly  influenced 
b>  changing  the  relative  strengths  of  these  interactions.  In 
order  to  test  this  hypothesis,  the  surface  of  one  set  of 
porous  disks  (<d>  =  5,8  nm)  was  chemically  treated  such 
that  hydrophobic  OCH3  groups  replaced  the  OH  groups. 
This  change  significantly  reduces  the  surface-guest 
interaction,  and  should  therefore  lead  to  a  reduction  of  the 
disorder  in  the  confined  solid.  The  measured  tunneling 
spectra  for  methyl  iodide  adsorbed  in  both  the  treated  and 
untreated  hosts  (<d>  =  5.8  nm)  are  shown  in  Fig.  2.  As 
seen  by  the  prominence  of  the  narrower  peaks,  the  more 
ordered  component  of  the  tunneling  spectrum  increases 
substantially  for  the  methyl  iodide  confined  in  the  glass 
with  the  hydrophobic  surface  as  compared  to  the  host  with 
the  hydrophilic  OH  groups.  In  fact,  the  disordered  fraction 
decreases  by  about  25  %,  which  indicates  that  the  thick- 
ness of  the  highly  disordered  layer  is  roughly  two-thirds 
that  of  the  untreated  glass.  Moreover  the  width  of  the 
tunneling  line  in  the  highly  disordered  component  has  also 
decreased  substantially  indicating  that  the  structure  of  the 
methyl  iodide  near  the  hydrophobic  surface  is  more  bulk- 
like than  it  is  near  the  hydrophilic  surface. 
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Fig.  2.  Tunneling  spectrum  of  confined  CH3I  in  5.8  nm  pores  with  and 
without  surface  treatment  to  make  the  surface  hydrophobic.  The 
hydrophobic  surface  results  in  a  diminished  disordered  component. 

In  summary,  the  sensitivity  of  tunneling  spectroscopy 
to  small  changes  in  the  rotational  potential  has  allowed  us 
to  probe  the  importance  of  surface-guest  interactions  in 
determining  the  disorder  imposed  by  the  confinement  of 
methyl  iodide  within  mesoporous  glasses.  It  has  been 
shown  that  the  adsorbed  species  in  the  hydrophilic  pores 
form  a  highly  disordered  layer  with  a  thickness  of  approxi- 
mately three  molecular  layers  and  that  a  more  ordered 
region  exists  in  the  pore  center  Interestingly,  the  molecu- 
lar structures  of  these  two  regions  do  not  depend  strongly 
on  either  the  average  pore  size  or  on  details  of  the  pore  size 
distribution.  When  the  surface  chemistry  is  altered  so  as 
to  reduce  the  importance  of  surface-guest  interactions, 
both  the  thickness  and  molecular  scale  disorder  of  the 
region  near  the  wall  are  greatly  reduced. 
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The  '^Bucket  Brigade" 
Mechanism  of  Proton  Diffusion 
in  Superprotonic  Conductors 


|roton  diffusion  plays  a  central  role 
in  many  processes  ranging  from  the 
regulation  of  biological  functions  to  the  produc- 
tion of  electricity  in  fuel  cells.  The  renewed 
emphasis  on  developing  the  next  generation  of 
portable  sources  of  power  requires  the  development 
of  lightweight  protonic  conducting  materials  for  use  as 
electrolytic  membranes  [1].  Furthermore,  understanding 
the  atomic-scale  mechanism  of  protonic  diffusion  in 
candidate  materials  is  necessary  to  optimize  their  perfor- 
mance. Since  neutrons  are  particularly  sensitive  to 
hydrogen,  neutron  scattering  is  an  extremely  valuable 
probe  of  proton  dynamics  (including  diffusion).  Moreover, 
because  the  wavelength  of  neutrons  is  comparable  to  inter- 
atomic distances,  neutron  scattering  is  capable  of  revealing 
the  atomic-scale  geometry  of  protonic  motions  in  a  way 
that  is  unmatched  by  other  spectroscopic  probes.  Compu- 
tational techniques  also  provide  detailed  atomistic  informa- 
tion about  proton  incorporation  and  dynamics.  Because  of 
the  simplicity  of  the  scattering  mechanism,  computational 
results  can  be  directly  compared  to  neutron  scattering  data. 
Thus  neutron  scattering  and  computational  methods  are  a 
powerful  combination  for  understanding  the  performance 
of  protonic  conducting  materials  [2].  Here  we  briefly 
describe  neutron  vibrational  spectroscopy  (NVS)  and 
quasielastic  neutron  scattering  (QNS)  measurements  on 
superprotonic  conductors  of  the  general  formula  MHXO4 
(M  =  Rb,  Cs,  X  =  S,  Se,  etc.)  that  reveal  detailed  informa- 
tion on  the  proton-conduction  mechanism. 

The  crystal  structure  of  alkali-metal  hydrogen 
sulfates  and  selenates  consists  of  chains  of  XO4  tetrahedra 

94  M>i 


(Fig.  1),  which  are  linked  by  hydrogen  bonds,  0-H-O. 
Due  to  the  proximity  of  the  proton  site  to  the  midpoint  of 
the  0-0  separation,  it  is  believed  that  the  proton  is  dynami- 
cally disordered,  yielding  symmetric  hydrogen  bonds. 
However,  at  a  sufficiently  low  temperature,  one  expects 
that  the  protons  may  order,  giving  rise  to  observable 
changes  in  the  vibrational  spectrum.  Above  room  tempera- 
ture, there  is  an  orientational  order-disorder  transition  to  a 
phase  in  which  the  tetrahedra  are  dynamically  disordered 
and  which  displays  a  very  large  protonic  conductivity  that 
should  be  observable  by  QNS. 

Figure  2  summarizes  NVS  resuhs  on  RbHS04 
measured  using  the  Filter  Analyzer  Neutron  Spectrometer. 
At  low  temperatures,  the  NVS  spectrum  shows  sharp 
peaks  at  83.3  meV  and  99.7  meV  The  temperature 
dependence  of  these  peaks  is  quite  unusual.  The  intensities 
decrease  very  rapidly  with  increasing  temperature  and 
almost  disappear  at  temperatures  as  low  as  200  K.  In 
addition,  while  the  high-energy  mode  softens,  as  expected, 
the  low-energy  mode  hardens  by  about  2  meV  as  the 
temperature  is  increased.  Near  200  K,  these  two  modes 
collapse  to  form  a  broad  feature  around  95  meV.  This 
suggests  that  there  is  significant  proton  motion  even  at  this 
low  temperature.  We  believe  that  this  unusual  behavior  is 
related  to  dynamically  disordered  hydrogen  bonds  arising 
fi-om  breaking  and  reforming  0-H  bonds.  Careful  neutron 
diffraction  studies  in  deuterated  materials  are  underway  to 
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Fig.  1.  Projection  of  tiie  MHXO4  structure,  showing  the  zigzag  pattern 
of  the  hydrogen  bonds  between  XO4  tetrahedra.  Yellow  dashed  lines 

indicate  the  hydrogen  bonds. 


Fig.  2.  Temperature-dependent  NVS  measurements  of  RbHS04,  As 
discussed  above,  these  results  suggest  very  large  amplitude  proton 
motions  are  present  at  temperatures  as  low  as  %  100  K.  The  right  panel 
illustrates  the  dynamically  disordered  hydrogen  bonding,  which  is  most 
likely  responsible  for  this  unusual  temperature  dependence. 
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look  for  significant  changes  in  the  apparent  locations  or 
thermal  factors  of  the  protons  at  these  temperatures.  It 
should  be  noted  that  since  this  type  of  dynamical  disorder 
is  localized  to  two  neighboring  oxygen  ions,  it  cannot  be 
responsible  for  long-range  protonic  conduction. 

To  probe  protonic  diffusion  directly,  the  temperature 
dependence  of  the  QNS  was  measured  for 
CsH(S04)i.^.  (Se04),-  (x  =s  0.24)  above  room  temperature 
using  the  Fermi-Chopper  time-of-flight  Spectrometer.  The 
data  were  fit  as  a  function  of  the  scattering  vector  Q 
(assuming  both  elastic  and  quasielastic  components)  in 
order  to  extract  the  Elastic  Incoherent  Structure  Factor 
(EISF),  which  is  defined  as  the  fraction  of  the  total 
scattering  that  is  elastic.  Figure  3  shows  the  EISF  as  a 
function  of  Q  for  several  different  temperatures.  The  fact 
that  a  non-zero  EISF  is  observed  even  at  475  K  indicates 
that,  within  the  time  window  of  the  measurement 
(<  20  ps),  the  hydrogen  motion  is  localized  rather  than 

diffusive. 
Moreover,  the 
fact  that  the 
EISF  changes 
with  tempera- 
ture indicates 
the  existence  of 
non-energeti- 
cally  equivalent 
proton-sites. 
The  best  fit  to 
the  EISF 
(shown  by  the 
solid  lines  in 
Fig.  3)  was 

obtained  assuming  two  types  of  hopping  motions  for  the 
protons  between  pairs  of  inequivalent  sites.  The  agreement 
between  the  data  and  the  model  is  excellent.  The  two  jump 
distances,  1.4  A  and  4.5  A,  are  essentially  independent  of 
temperature  (see  Fig.  4).  Looking  at  the  structure,  these 
can  easily  be  identified  as  the  reorientation  of  the  SO4-H 
group  around  the  quasi-threefold  (1.4  A)  and  the  quasi- 
twofold  (4.5  A)  axes  (Fig.  4).  The  ratio  of  the  times  that 
the  proton  spends  in  each  of  the  two  sites  depends  expo- 
nentially on  temperature  as  expected  (solid  lines  in  the 
upper  left  panel  of  Fig.  4).  For  the  short-distance  reorien- 
tation, the  time  ratio  goes  rapidly  to  unity  around  405  K, 
just  before  the  super-protonic  phase  transition,  indicating 
that  the  two  sites  are  equally  populated.  For  the  long- 
distance reorientation,  the  ratio  approaches  0.2  at  the 
highest  temperature. 
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Fig.  3.  EISF  (dots)  and  fits  (solid  lines)  at  different 
temperatures  for  CsH(S04)i.^(Se04)^  (x  0.24). 
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Fig.  4.  Top:  Temperature  dependence  of  the  ratios 
of  the  residence  times    to  Tj  for  the  pair  of  two- 
site  hopping  motions  identified  by  QNS. 
Bottom:  The  two  jump  distances  obtained  from 
QNS  as  a  function  of  temperature.  These  distances 
correspond  to  the  two  reorientational  motions 
shown  on  the  right.  The  red  circles  correspond  to 
the  motion  depicted  on  top  while  the  blue  squares 
correspond  to  the  motion  shown  on  the  bottom. 


These  findings 
are  consistent 
with  a  previ- 
ous neutron- 
powder- 
diffraction 
study  on  the 
high-tempera- 
,ture  phase  of 
CSHSO4, 
which  revealed 
four  different 
orientations  of 
each  tetrahe- 
dron [3]. 
Most  probably 
the  four 
orientations 
mimic  the 
dynamical 

orientational  disorder  of  the  tetrahedra  that  we  observe  by 
QNS.  We  also  note  that  the  two  types  of  dynamic 
orientational  disorder  shown  in  Fig.  4  do  not  result  in  long- 
range  protonic  conduction  since  the  motions  are  localized 
on  a  single  tetrahedron.  However,  when  one  combines 
these  two  motions  with  the  dynamically  disordered 
hydrogen  bond  implied  by  NVS,  the  mechanism  of  pro- 
tonic conduction  can  be  understood.  Namely,  the  protons 
move  between  two  neighboring  tetrahedra  via  dynamically 
disordered  hydrogen  bonds.  The  tetrahedra  then  reorient, 
allowing  the  protons  to  move  to  a  third  tetrahedron,  again 
via  the  dynamically  disordered  hydrogen  bonds.  Thus  the 
tetrahedra  play  the  role  of  a  "bucket  brigade"  by  accepting  a 
proton  from  an  adjacent  tefrahedron  then  turning  and  handing 
it  the  tetrahedron  on  the  other  side.  On  the  atomic  scale, 
superprotonic  conductivity  occurs  in  the  MHXO4  family  of 
compounds  due  to  this  efficient  mechanism  of  transporting 
protons  through  the  lattice. 
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Motion  in  Intermediates  Along  the 
Protein-Folding  Pathway:  cytochrome  c 


P roteins  are  long  macromolecular 
chains  consisting  of  amino  acid 
repeat  units  that,  in  the  native  state,  adopt  a 
well  defined  and  intricately  folded  three-dimen- 
sional conformation.  The  ability  to  predict  a  priori 
the  native  three-dimensional  structure  of  a  given 
amino  acid  sequence,  the  so  called  "protein  folding 
problem,"  is  essential  for  the  full  promise  of  genetic 
therapies  to  come  to  fruition  [1].  In  fact,  many  apparently 
unrelated  diseases  (e.g.,  Alzheimer's,  Mad  Cow,  cystic 
fibrosis,  and  even  many  cancers)  result  from  protein 
folding  gone  wrong. 

In  order  to  obtain  detailed  knowledge  of  the  protein 
folding  mechanism,  extensive  characterization  of  the 
relevant  interactions,  structure  and  dynamics  is  required  at 
each  step  of  the  protein-folding  pathway.  To  this  end, 
researchers  have  been  investigating  the  properties  of 
proteins  in  the  native,  partially  unfolded  and  the  completely 
unfolded  states  to  obtain  "snapshots"  of  the  relevant 
behavior  of  intermediates  along  the  folding  pathway, 
yielding  information  intended  to  facilitate  rationalization  of 
folding  mechanisms. 

To  deepen  our  understanding  of  the  molecular 
motions  involved  in  protein  folding  and,  hence,  the  role  of 
dynamic  behavior  on  the  overall  folding  mechanism,  we 
employed  quasielastic  neutron  scattering  (QENS)  as  a 
dynamical  probe  of  several  folded  states  of  the  electron 
transfer  protein  cytochrome  c  (see  Fig.  1).  QENS  provides 


Fig.  1.  The  native  state  of  cytochrome  c. 
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Fig.  2.  Illustration  of  the  general  quasielastic  neutron  scattering  profile 
of  a  protein  in  solution  indicating  the  narrow  scattering  from  the  global 
diffusion  and  the  broad  scattering  due  to  the  internal  side  chains  of  the 
amino  acid  residues. 

a  method  for  directly  probing  the  internal  dynamics  of 
biomolecules  on  the  picosecond  time  scale  and  Angstrom 
length  scale.  In  this  technique  neutrons  exchange  energy 
and  momentum  with  the  nuclei  of  the  sample,  allowing  the 
forms,  geometry,  amplitudes  and  time-dependence  of  the 
motions  involved  to  be  determined.  By  making  use  of  the 
large  incoherent  cross  section  of  hydrogen  nuclei  and  the 
capabilities  of  the  Disk-Chopper  Spectrometer  at  «  32  |aeV 
resolution  it  was  possible  to  collect  data  on  five  different 
cytochrome  c  folded  states  (in  order  of  increasing  unfold- 
ing): native  (N),  salt  induced  molten  globule  (MG),  alkaline 
denamred  (D^u^),  acid  denatured  (D^^^^  and  guanidine 
hydrochloride  unfolded  (U). 

In  a  typical  QENS  experiment  of  this  sort,  a  mono- 
chromatic beam  of  neutrons  is  scattered  by  a  buffered 
solution  containing  the  protein,  the  contribution  of  the 
scattering  from  the  buffer  is  removed,  and  the  resulting 
profile  is  considered  to  be  that  resulting  from  just  the 
protein  [2,  3].  Figure  2  depicts  the  general  neutron 
scattering  profile  from  a  protein  in  solution,  illustrating  the 
distinct  scattering  from  global  (bulk  translational  and 
rotational  diffusion)  and  internal  (protein  side-chain) 
dynamics.   A  profile  of  this  sort  is  analyzed  with  a  two 
Lorentzian  fit  model  —  a  narrow  Lorentzian  to  fit  the 
global  motion  and  a  broad  Lorentzian  to  fit  the  internal 
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motions.  The  two  Lorentzian  fit  model  permits  the  average 
ensemble  time  scale  of  global  and  internal  side-chain 
motions  to  be  ascertained  separately. 

The  half  width  at  half  maximum  of  both  Lorentzian 
fit  functions,  corresponding  to  the  global  and  internal 
motions  as  a  function  of  scattering  vector  squared  {q-)  are 
given  in  Fig.  3.  An  effective  diffusion  coefficient  (Deff)  is 
determined  from  the  slope  of  -Tgiobai     a  function  of  q-. 
As  evident  from  the  tabulated  results,  D^^  generally  de- 
creases as  the  protein  unfolds  reflecting  the  correspond- 
ingly higher  specific  molecular  volume  of  non-native 
folding  states.  The  exception  lies  with  the  acid  denatured 
state  that,  due  to  the  sfrong  electrostatic  interactions 
present  in  the  system,  shows  a  near  native  D^^. 


0.04 


I  '  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


State  D„ff/(10-^cm2/s) 


A  decrease  in  the  relative  time  scale  of  the  internal 
motions  is  observed  as  a  function  of  protein  unfolding.  It 
is  believed  that  the  rather  close-packed  native  structure 
provides  more  sharply  defined,  short-range  interactions 
between  the  proximate  side  chains  of  the  amino  acid 
residues.  Upon  unfolding  of  the  protein,  these  interactions 
are  diminished  by  the  increased  separation  between  side 
chains  and  small,  weakly  interacting  buffer  molecules 
between  them. 

The  results  of  this  study  demonstrate  the  type  of 
dynamic  information  that  can  be  determined  from  QENS 
experiments  performed  on  a  protein  in  several  different 
structural  states.  As  the  nature  of  the  unfolding  medium 
displays  a  critical  role  in  the  internal  dynamics  of  the 
unfolded  states  of  cytochrome  c,  any  attempts  to  correlate 
features  between  different  folding  states  should  take  this 
into  account.  In  combination  with  other  characterization 
techniques,  information  of  this  sort  will  prove  beneficial  to 
the  understanding  and  development  of  accurate  protein 
folding  mechanisms. 
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Fig.  3.  Plot  of  the  half  width  at  half  maximum  of  the  global  (top,  native 
only)  and  internal  (bottom)  Lorentzian  as  a  function  of  scattering  vector 
squared.  The  tabulated  effective  diffusion  coefficients  are  given  for  each 
folded  state  of  cytc. 
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Neutron  Reflectivity  from 
Biomimetic  Membranes:  improving 
Reliability  through  Wavelet  Analysis 


Pj  H  nderstanding  such  key  biological 
processes  as  molecular  recognition, 
protein  insertion,  and  molecular  self-assembly 
in  Irving  biological  membranes  remams  a  chal- 
lenge. Membrane  materials  mimicking  biological 
ones  (biomimetic  membranes)  provide  model  systems 
to  address  this  issue,  and,  because  of  its  special 
sensitivity  to  hydrogen,  neutron  reflectivity  (NR)  offers  a 
unique  way  to  reveal  thin  film  structures  in  biomimetic 
membranes. 

The  feasibilits'  of  using  phase-inversion  techniques  in 
NR  to  reveal  such  structural  details  has  been  demonstrated 
[1].  It  is  possible  to  directly  measure  the  real  part,  ri(0), 
of  the  complex  reflection  coefficient  r(Q)  as  a  function  of 
wavevector  Q  and  to  mathematically  invert  it  to  obtain  the 


0  40  80  120 

x//o 

Fig.  1.  Diagram  of  model  (top)  for  which  the  computed  scattering  length 
density  p  (x)  is  the  thick  grey  curve.  Trends  converge:  p.^(x), 
pjx),  .  .  .,^p(x). 


scattering  length  density  (SLD)  depth  profile  p(x)  of  the 
film.  There  remains  much  to  do,  however,  in  defining  the 
reliability  of  the  results,  especially  in  accessing  the  un- 
avoidable reduction  of  spatial  resolution  induced  by  data 
truncated  at  a  maximum  wavevector,  Q^^x-  Wavelet 
analysis  provides  a  systematic  and  useful  approach  to  the 
problem  [2]. 

In  using  Wavelet  MultiResolution  Analysis  (WMRA), 
SLD  profiles  are  characterized  by  a  fmest  length  scale 
4  ~  1 A  and  coarser  scales,  ^  =  2-j£q,  with  j  negative.  p{x) 
is  viewed  as  a  coarse  "trend"  pj{x)  at  resolution  level  J, 
with  added  "detail"  Ajp(x)  giving  the  trend  at  the  next  finer 
scale,  Pj+i(x)  =  Pj(x)+Ajp(x).  Relative  to  a  base  scale  of 
resolution  ^j,  the  SLD  profile  thus  can  be  represented  by 
the  trend  plus  all  remaining  detail,  p(x)  =p^(xj+Xz\  p(x). 
The  experimental  p(x)^  associated  with  a  given  ^max 
determined  by  the  instrument,  is  a  blurred  representation  of 
the  veridical  SLD.  It  can  be  thought  of  as  a  coarse  image 
bracketed  by  neighboring  trends  for  this  2max- 

WMRA  provides  spatially  localized  orthonormal  bases 
for  this  description.  A  family  of  wavelets  called 
Daubechies-8  seems  well  suited  to  NR  analysis.  For 
illustration,  we  use  a  realistic  SLD  profile  obtained  by 
molecular  modeling  to  represent  a  hybrid  lipid  membrane 
on  a  thin  gold  film  (diagram  on  top  of  Fig.  1),  a  biomimetic 
system  typical  of  those  being  studied  in  many  laboratories. 
The  model  p(x),  seen  in  the  back  panes  of  Figs.  1  and  2, 
consists  of  the  gold  layer,  a  hydrogenated  alkanethiol  layer, 
and  a  deuterated  lipid  monolayer. 

Figures  I  and  2  depict  the  convergence  of  the 
WMRA  descriptions  of  p(x)  as  trend  and  detail.  The 
"overall"  shape  of  p(x)  effectively  is  determined  by  the 
trend  p.4ix)  and  detail  A_^p{x),  i.e.,  by  trend  p.^(x). 
However,  emergence  of  the  prominent  double  peaked 
structure  of  the  lipid  head  group  near  x/1q  =100  needs  detail 
A_^p(x).  Subsequent  detail  mainly  acts  to  sharpen  the  edges 
between  the  film's  components. 

Figure  3  shows  the  effective  contributions  of  the 
trend  and  the  successive  spatial  detail  to  the  reflection 
coefficient  r(Q),  each  calculated  exactly.  The  edge- 
sharpening  detail  seen  in  Fig.  3  is  not  revealed  in  the 
reflection  spectrum  below  Q^q~  0.6.  Figure  4  shows  the 
"smeared"  p(x)  obtained  by  direct  inversion  of  the  reflec- 
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Fig.  2.  Trend  plus  detail  converges:  Pi(x)+A  iP(x)+A  3p(x),  .  .  .,- 
p{x). 
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Fig.  3.  Effective  contributions  to  Q^r^(Q)  generated  by  the  trend  p4(x) 
(in  red)  and  the  details  A^p  (x)  for  /  =  -4,  -3,  -2,  -1  shown  in  Fig.  2.  The 
Q^r^(Q)  contributions  are  labeled  by  the  SLD's  that  produced  them. 


tion  from  p(x)  using  data  truncated  at  Q^zx-  0-2  A-^.  This 
is  seen  to  fall  "between"  the  trends  p.4ix)  and  p.^ix), 
stemming  from  low-pass  filters  with  roU-offs  at 
Q  w  0.2       and  Q  «  0.4  A-',  respectively.  This  is  expected 
from  the  fact  that  the  Fourier  transforms  of  trends  and 
detail  overlap  to  a  degree.  Thus,  a  "pure"  trend  cannot  be 
observed  in  the  truncated  data. 
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Fig.  4.  Inverted  r^(Q)  (black)  using  only  data  truncated  at  Q/o  =  0.2. 
Trends  p^{x)  (red)  and  p.-^i^)  (blue)  of  the  actual  p(x). 


Wavelet  analysis  thus  provides  a  systematic  method 
for  assessing  the  correctness  of  density  profiles  measured 
by  NR.  It  promises  to  add  reliability  in  unraveling  struc- 
tures of  importance  in  biomimetic  membranes. 
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Protein  Adsorption  to 
IMetal-Cheiating  Lipid  Monolayers 


he  adsorption  of  proteins  to 
surfaces,  especially  lipid  layers, 
plays  a  critical  role  in  biochemical  processes 
within  living  systems.  These  include  membrane 
signaling  and  recognition,  toxin  detection,  and 
mitigating  the  effects  of  toxin/biowarfare  agents. 
Specific  issues  to  be  addressed  are:  protein  orientation 
and  surface  density,  monolayer  versus  multilayer  adsorp- 
tion, specific  versus  non-specific  interactions,  and  lateral 
order  and  conformational  changes  upon  adsorption 
including  denaturization  of  protein  [1]. 

In  this  work  neutron  reflection  is  combined  with 
grazing  incidence  x-ray  diffraction  (GIXD)  to  study  the 
interactions  of  the  protein  myoglobin  with  Langmuir 
monolayers  of  synthetic  lipids.  Neutron  reflection  is  used 
to  determine  the  adsorbed  amount  and  to  obtain  informa- 
tion regarding  the  orientation  and  conformation  of  the 
adsorbed  myoglobin  [2].  GIXD  is  used  to  study  the 
response  of  the  lipid  layer  that  occurs  upon  protein 
binding. 

For  lipid  layers  we  have  used  the  synthetic  lipid 
distearyl  imino-diacetate  (DSIDA)  that  contains  receptors 
for  the  metal  ion  Cu2+.  The  use  of  metal  ion  coordination 
to  target  the  adsorption  of  proteins  to  lipid  membranes  has 
been  studied  extensively.  This  method  uses  specific 
coordination  interactions  between  Cu2+  and  naturally 
occurring  histidine  units  in  myoglobin. 

A  Langmuir  monolayer  of  100  %  DSIDA  was  spread 
onto  the  surface  of  buffered  D2O  or  H2O  (Fig.  1 )  and  then 
compressed  to  a  surface  pressure  of  »  40  mN/m  where  it 
forms  a  solid  condensed  phase  layer  (see  GIXD  data 
below).  A  dilute  solution  of  CuCh  was  injected  under  the 
condensed  DSIDA  monolayer.  The  Cu2+  ions  chelate  at  the 
imino-diacetate  site  as  indicated  in  Fig.  1 . 

Neutron  reflectivity  data  for  the  condensed  lipid  layer 
with  added  Cu2+  ions  are  shown  in  Fig.  2  for  both  the  HjO 
and  D2O  buffered  subphases.  Myoglobin  was  then 
injected  into  the  subphase  and  neutron  reflectivity  data 
were  taken  repeatedly  to  monitor  the  time  dependence  of 
protein  adsorption.  The  adsorbed  protein  layers  reached 
quasi-equilibrium  at «  14  h  after  adding  myoglobin  into  the 
solution.  Only  final  sets  of  reflectivity  data  with  adsorbed 
protein  (for  both  H2O  and  DjO  subphases)  are  shown  in 
Fig.  2.  The  very  large  change  in  the  reflectivity  after 
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Fig.  1.  Schematic  of  histidine  groups  on  myoglobin  adsorbing  to 
distearyl  imino-diacetate  (DSIDA)  molecules  at  chelated  Cu^^  ion  sites. 

protein  injection  indicates  a  significant  adsorption  of 
protein  under  the  lipid  layer. 

Solid  lines  in  Fig.  2  are  fits  to  the  reflectivity  data 
from  which  the  segment  concentration  profile  of  protein 
can  be  obtained.  Simultaneous  fits  to  both  sets  of  data 
show  that  the  dimension  of  the  quasi-equilibrium  adsorbed 
myoglobin  layer  is  36  A  +/-  2  A,  and  that  the  segment 
volume  fraction  of  protein  is  ~  50  %.  Unit  cell  dimensions 
of  myoglobin  from  its  crystal  structure  are  (a,  b,  c)  = 
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myoglobin 
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Fig.  2.  Neutron  reflectivity  for  100  %  DSIDA  on  buffered  DjO  and  HjO 
subphases  with  and  without  myoglobin  in  the  subphase.  Solid  lines  are 
fits  to  the  data  as  described  in  the  text. 
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Fig.  3.  Surface  diffraction  peal<  at  a  constant  surface  pressure  of  40  mN/m:  (a)  1 00  %  DSIDA  and  also  with 
Cu2^  (b)  1.9  h  (no  peai<)  and  4.4  h  (peak  returns)  after  myoglobin  injection. 


(64.84,  30.98,  34.92)  A.  Thus,  the  protein  in  all  likelihood 
is  bound  to  the  lipid  layer  in  its  native  state,  i.e.,  not 
denatured. 

If  no  CUCI2  is  added  to  the  subphase,  then  no  change 
in  the  neutron  reflectivity  signal  is  observed  upon  injection 
of  myoglobin.  Alternatively,  if  we  inject  a  strong  Cu 
scavenging  agent  (EDTA)  in  the  subphase  after  the  protein 
adsorption,  the  reflectivity  reverts  back  to  that  of  the  pure 
condensed  lipid  layer.  This  clearly  shows  that  we  are 
observing  very  specific  binding  of  the  protein  to  the 
chelated  Cu2+  sites  on  DSIDA. 

Important  insight  has  been  revealed  in  complemen- 
tary GIXD  data  from  the  DSIDA  system.  Figure  3a  shows 
the  surface  diffraction  peak  from  the  tails  of  a  condensed 
DSIDA  monolayer  at  an  in-plane  wavevector  Q^y=  1.50  A-'. 
This  corresponds  to  a  hexagonal  packed  lipid  layer  with  a 
nearest  neighbor  spacing  of  4.21  A.  Insertion  of  Cu^+ions 
into  the  layer  does  not  seem  to  affect  the  2-D  crystalline 
order  of  the  lipid  layer  except  for  a  slight  change  in  lattice 
spacing. 

When  myoglobin  is  added  to  the  subphase  at  a 
constant  surface  pressure  of  40  mN/m,  the  diffraction 
peak  arising  from  the  hexagonal  packing  of  the  DSIDA 
tails  (at  Q^^=  1.50  A-i)  is  drastically  reduced  at  short  times 
(Fig.  3b).  This  indicates  a  large  perturbation  in  the  2-D 
crystallinity  of  the  lipid  layer  upon  insertion  of  protein.  The 
surface  diffraction  peak,  however,  reappears  after  a  few 
hours  (Fig.  3b).  The  peak  is  slightly  shifted  to  a  higher 
value  of  Q^y,  and  the  distribution  of  intensity  normal  to  the 
surface  (not  shown)  indicates  a  substantially  increased  tilt 
of  the  lipid  tails.  By  contrast,  no  change  in  the  diffraction 
peak  is  observed  after  addition  of  myoglobin  when  the 
surface  layer  is  maintained  at  a  constant  area  (after  initial 
compression  to  40  mN/m). 


We  hypothesize  that  the 
protein  is  able  to  insert  into  the 
lipid  monolayer  when  it  is 
maintained  at  a  constant 
pressure  of  40  mN/m,  but  does 
not  insert  when  at  constant  area 
(surface  pressure  in  this  case 
rises  from  40  mN/m  to 
43  mN/m).  Thus  far,  all  the 
neutron  reflection  data  have 
been  obtained  at  constant 
surface  area.  In  future  work 
we  will  examine  adsorption  to 
lipid  monolayers  composed  of 
100  %  DSIDA  at  constant  surface  pressure.  We  expect  to 
see  very  different  arrangement  of  the  protein  adsorbed 
layers  for  conditions  at  which  the  protein  inserts  into  the 
lipid  monolayer  versus  when  it  does  not. 

In  summary,  a  combination  of  neutron  reflectivity 
and  GIXD  has  proven  to  be  a  very  powerful  way  to  study 
both  the  protein  adsorption  to  lipid  monolayers  as  well  as 
structural  changes  that  the  lipid  layer  itself  undergoes  as 
the  protein  is  adsorbed.  Because  of  different  contrast 
conditions  for  neutrons  and  x-rays,  neutron  reflectivity 
profiles  are  very  sensitive  to  the  adsorption  of  the  protein, 
while  GIXD  is  ideally  suited  for  examining  in  detail  any 
changes  in  the  2-D  ordered  structure  of  the  lipid  layer  itself 

This  successful  study  opens  up  a  rather  large  field  to 
examine  other  important  lipid-protein  interactions  such  as 
insertion  of  biological  toxins  into  a  lipid  layer.  Preliminary 
studies  of  neutron  reflectivity  and  GIXD  to  examine  the 
adsorption  of  cholera  toxin  to  monolayers  of  glycolipids  (a 
class  of  lipids  containing  sugar  molecules  which  are  found  on 
cell  surfaces)  show  extraordinary  promise  in  this  regard. 
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Hydration  State  of  Single 
cytochrome  c  Monolayers  on  Soft 
Interfaces  via  Neutron  Reflectivity 


pter  is  critical,  not  only  for  the 
correct  folding  of  proteins  but 
also  for  the  maintenance  of  this  folded  struc- 
ture. The  internal  molecular  motions  m  proteins, 
which  are  necessary  for  biological  activity,  are  very 
dependent  on  the  degree  of  flexibility  which  is 
determined  by  the  level  of  hydration.  The  nmnber  of 
water  molecules  hydrating  a  functioning  protein  is  thus 
an  issue  of  great  interest.  Given  an  appropriate  model 
system,  neutron  scattering,  particularly  sensitive  to 
hydrogen,  provides  a  way  to  measure  this  water  content. 

Previous  optical  spectroscopy  studies  have  shown 
that  yeast  cytochrome  c  (YCC)  covalently  bound  to  a  soft 
interface  (Fig.  1)  and  partially  hydrated  by  a  moist  helium 
atmosphere  can  be  fully  fimctional  with  respect  to  the 
oxidation-reduction  chemistry  of  its  iron  porphyrin  pros- 
thetic group  [1,  2,  5,  6].  This  system  thus  provides  an 
opportunity  to  examine  the  water  distribution  required  to 
maintain  the  structvire  and  function  of  the  protein  YCC. 

Unlike  x-rays,  neutrons  scatter  very  differently  from 
hydrogen  and  deuterium.  This  fact  allows  the  water 
distribution  in  the  monolayer  profile  structure  of  YCC  to  be 
obtained  by  comparison.  The  neutron  scattering  length 
density  profile  derived  from  neutron  reflectivity  for  a  YCC 
monolayer  hydrated  by  D^O  is  compared  to  the  profile  for 
identical  hydration  with  HjO. 

The  silicon  surface  layer  of  an  iron-silicon  (Fe/Si)  or 
iron-gold-silicon  (Fe/Au/Si)  multilayer  solid  substrate  can 
bind  a  self-assembled  monolayer  (SAM)  to  form  the  soft 


surface  used  to  tether  the  YCC  monolayer  (See  diagrams 
at  the  top  of  Fig.  3.)  Such  a  multilayer  substrate  has  two 
key  advantages  in  a  neutron  (or  x-ray)  reflectivity  experi- 
ment: a  multilayer  substrate  dramatically  enhances  this 
scattering  for  momentum  transfer  normal  to  the  substrate 
surface,  and  a  multilayer  substrate  also  provides  an 
important  reference  profile  structure  for  the  unique 
interferometric  phasing  of  the  reflectivity  data. 

Neutron  reflectivities  (Fig.  2)  were  collected  on  the 
NG-1  reflectometer  for  both  H2O  and  D2O  hydration  cases 
for  two  such  SAMA^CC  samples.  One  SAM  formed  a 
nonpolar  surface  (-CH3/-SH  =  6:1  mixed  endgroups),  and 
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Fig.  1 .  Molecular  dynamics  "snapshot"  of  YCC  tethered  to  the  soft 
surface  of  an  uncharged-polar  self-assembled  monolayer  (SAM). 


Fig.  2.  Top:  normalized  reflectivity  data  with  incident  neutron  spins 
parallel  to  the  iron  magnetization  for  the  bare  multilayer  substrate, 
substrate  plus  nonpolar  SAM  plus  YCC/DjO,  and  substrate  plus 
nonpolar  SAM  plus  YCC/HjO.  Bottom:  data  for  similar  arrangements 
with  the  uncharged  polar  SAM  on  a  Fe/Au/Si  substrate.  In  the  top  plot 
the  H2O  data  are  offset  by  150  units  and  the  DjO  data  by  300  units  on  the 
ordinate,  and  in  the  bottom  plot  the  HjO  data  are  offset  by  75  units  and 
the  D2O  data  by  150  units  on  the  ordinate.  (Note:  IstnOIX.) 
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the  other  SAM  formed  an 
uncharged  polar  surface 
(-0H/-SH  =  6:1  mixed 
endgroups).  These  data 
were  analyzed  using  a  new 
interferometric  phasing 
method  that  makes  use  of 
two  features:  the  neutron 
scattering  contrast  be- 
tween the  Si  and  Fe  layers 
in  a  single  reference 
multilayer  structure,  and  a 
constrained  refinement 
approach  using  the  finite 
extent  of  the  gradient  of 
the  profile  structures  for 
the  systems  [3].  The 
water  distribution  profiles 
for  the  two  SAIVLTCC 
monolayers  provided  by 
this  analysis  are  shown  in 
Figs.  3a  and  3b. 

For  hydration  with 
D2O,  these  profiles  show 
that  the  protein  monolayer 
is  3  A  to  4  A  closer  to  the 

substrate  surface  for  the  uncharged-polar  SAM  compared 
to  the  nonpolar  SAM.  This  finding  is  in  excellent  agree- 
ment with  simulations  of  these  systems  [3,  4],  and  arises 
because  residues  on  the  protein's  surface  interact  strongly 
with  the  polar  SAM's  hydroxyl  endgroups  via  hydrogen- 
bonding,  thus  drawing  the  YCC  closer  to  the  SAM  surface. 

Given  these  water  distribution  profiles,  the  number  of 
water  molecules  hydrating  the  YCC  monolayer  at  each 
SAM  can  be  calculated.  Allowing  for  proton  exchange  in 
the  cytochrome  c  molecule  itself  (  *  1 7  polypeptide  back- 
bone hydrogens  and  «  104  side  chain  hydrogens),  we 
obtained  values  of*  167  water  moleculesAfCC  at  the 
uncharged  polar  SAM  (exposed  to  He  at  81  %  relative 
humidity)  and  «  297  water  moleculesA' CC  at  the  nonpolar 
SAM  (exposed  to  He  at  88  %  relative  humidity)  with 
relative  errors  of  order  20  %  to  25  %.  These  findings 
allow  quantitative  comparison  to  molecular  models, 
opening  an  important  window  to  understanding  the  role  of 
water  in  protein  functioning. 


Fig.  3.  The  absolute  neutron  scattering  length  density  profiles  for  partial  hydration  with  DjO  and  H^O  and  their 
difference  profile  for  both  the  nonpolar  SAM  (a)  and  the  uncharged  polar  SAM  (b)  cases.  The  boundaries  for  the 
cytochrome  c  protein  region  of  the  profiles  used  for  calculation  of  the  amount  of  water  hydrating  the  protein  are 
z  =  1 0  A  and  z  =  60  A.  Schematics  of  the  composite  structures  are  shown  above  their  respective  scattering 
length  density  profiles  approximately  to  scale. 
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Direct  Measurement 
of  the  Reaction  Front  in 
Chemically  Amplified  Photoresists 


If 


he  semiconductor  industry  is  rapidly 
approaching  the  need  to  fabricate 
sub- 100  nm  structures  to  continue  perfor- 
mance increases  in  integrated  circuits.  Photoli- 
thography remains  the  enabling  technology  for  the 
fabrication  of  integrated  circuit  patterns.  Although 
industry  is  able  to  commercially  produce  130  nm 
features,  lithographic  materials  and  processes  are  not 
ftiUy  available  for  the  production  of  sub- 100  nm  struc- 
tures. New  materials  must  be  able  to  produce  structures 
with  dimensional  tolerances  of  2  nm  to  5  nm,  dimensions 
near  the  size  of  the  polymer  chain  molecules  in  the  imaging 
layer. 

Current  imaging  layers,  chemically  amplified  photore- 
sists, are  multi-component  materials  consisting  of  a 
polymer  resin  initially  insoluble  in  an  aqueous  base  devel- 
oper solution,  a  photoacid  generator  (PAG),  and  other 
additives.  A  schematic  diagram  of  the  fabrication  process 
and  an  example  structure  are  shown  in  Figs.  1  and  2.  The 
patterns  are  generated  in  the  resist  by  exposure  to  UV 
radiation  through  a  mask.  In  the  exposed  areas,  the  PAG 
decomposes,  forming  an  acid  species.  Upon  baking,  the 
acid  diffuses  and  catalyzes  a  deprotection  reaction  render- 
ing the  insoluble  resist  soluble  in  a  developer.  The  soluble 
regions  are  then  removed  with  the  aqueous  base  developer. 
Control  over  this  process  is  dominated  by  the  events  in  the 
transitional  region  between  exposed  and  unexposed  areas 
of  the  photoresist. 

The  initial  development  of  chemically  amplified 
photoresists  was  a  key  technological  breakthrough  further- 
ing the  continued  use  of  optical  photolithography  for  the 


Fig.  1 .  At  left,  SEM  image  of  a  litliograpliically  fabricated  structure  with 
a  nominal  critical  dimension  of  150  nm  and  300  nm  pitch.  At  right, 
schematic  illustration  of  patterning  through  a  mask  to  create  the 
structure. 


UV  radiation,  X<  248  nm 


1^     p  n  pif^jplr 


UUiUiUUUiii 


photoacid  generated 
base  soluble 


Reaction-diffusion  front 
profile  determines  both 
the  final  front  structure 
and  where  dissolution  is 
allowed 


§100% 


Q. 
CD 
Q 


0% 


Photoresist  film 


Silicon 


lithographic  interface 


Solubility  limit 


Fig.  2.  Schematic  diagram  of  the  reaction-diffusion  process  in 
chemically  amplified  photoresists.  The  transitional  region  between 
exposed  and  unexposed  regions  is  particularly  important  for  critical 
dimension  and  roughness  control. 

fabrication  of  sub  micrometer  features.  The  importance  of 
chemically  amplified  photoresist  concepts  is  illustrated  by 
the  industry-wide  use  of  these  materials  in  the  fabrication 
of  state-of-the-art  devices  today  and  into  the  foreseeable 
future.  With  the  imminent  need  for  sub- 100  nm  feature 
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Fig.  3.  General  deprotection  reaction  in  a  chemically  amplified 
photoresist.  Deuterated  molecular  parts  with  strong  contrast  to 
neutrons  are  shown  in  red. 
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Fig.  4.  (a)  schematic  diagram  of  the  processing  steps  including  the  application  of  each 
film  layer,  exposure  to  UV  radiation,  a  post-exposure  bake  (PEB)  and  development, 
(b)  and  (d)  show  the  NR  data  and  compositional  depth  profile  and  (c)  and  (e)  show  the 
XR  data  and  density  depth  profiles  from  each  step. 


sizes,  nanometer  level  control  over  the  spatial  evolution  of 
the  deprotection  reaction  front  and  the  subsequent  develop- 
ment steps  is  needed.  The  current  level  of  understanding 
of  the  acid-catalyzed  reaction-diffusion  process  is  not 
sufficiently  detailed  to  achieve  this  goal. 

To  meet  the  need  for  spatially  detailed  data,  we  have 
developed,  in  collaboration  with  the  IBM  T.  J.  Watson 
Research  Center  and  the  University  of  Texas  at  Austin,  an 
experimental  methodology  to  directly  measure  the  spatial 
evolution  of  the  deprotection  reaction  front  through  each 
processing  step.  By  using  model  photoresist  polymers 
with  a  deuterated  protection  group,  neutron  and  x-ray 
reflectivity  (NR,  XR)  measurements  are  able  to  follow  the 
compositional  and  density  profiles  of  the  reaction-diffusion 
process  in  an  idealized  transitional  region  in  a  bilayer  film 
stack.  Neutron  contrast  with  the  reaction  is  possible 
because  the  reaction  products  involving  the  deuterated 
protection  group  are  volatile,  as  illustrated  in  Fig.  3. 

The  NR  and  XR  resuhs  from  a  series  of  bilayer 
samples  are  shown  in  Fig.  4.  NR  and  XR  measurements 
were  taken  after  each  processing  step  from  the  spin- 
coating  of  each  film  layer  through  exposure  and  bake  to 
the  dissolution  of  the  upper  layer  and  any  deprotected 
polymer.  The  data  show  that  the  deprotection  reaction 
front  broadens  with  time  within  the  photoresist.  The  initial 


interfacial  width  was  approximately  4  nm  wide 
whereas  the  reaction  front  was  nearly  20  nm 
wide.  Upon  development,  the  data  show  that 
the  developed  film  remained  sharply  defined 
with  a  surface  roughness  of  approximately  4 
nm  as  well.  Although  the  reaction  front  was 
much  wider  than  the  final  interfacial  width,  the 
selectivity  of  the  dissolution  process  allows  for 
the  fabrication  of  well-defined  nanostructures 
with  diffusive  fransport  mechanisms. 

The  spatial  detail  afforded  by  the  devel- 
oped methodology  can  play  an  important  role 
in  not  only  the  qualitative  description  of  the 
reaction-diffusion  process,  but  also  with  the 
quantitative  determination  of  physical  param- 
eters and  changes  in  any  physical  or  chemical 
properties  over  the  nanometer  length  scales 
required  for  control  over  the  fabrication  of 
sub- 100  nm  structiires.  In  addition,  the 
strategy  developed  here  is  general  and  can  be 
adapted  for  the  study  of  candidate  photoresist 
materials  being  developed  for  future  applica- 
tions. 

"[This  work  [1]]  opens  a  window  of 
opportunity  to  construct  structure-property  relationships 
between  chemical  transport  mechanisms  and  ultimate  resist 
resolution.  It  may  also  lead  to  insights  into  the  ultimate, 
intrinsic  resolution  limits  and  critical  dimension  control  of 
polymer-based  imaging  materials."  [2] 
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Scaling  of  Membrane  Phase  Dynamics 


Membrane  phases  are  one  example 
of  the  variety  of  complex  fluid 
structures  formed  by  surfactant  molecules  in 
solution.  They  are  most  commonly  encountered 
in  biological  systems,  but  are  also  used  extensively 
in  a  variety  of  commercial  applications  from  personal 
care  products  to  foods  to  paints.  The  diagram  in  Fig.  1 
shows  the  regularly  stacked  lamellae  of  the  L„  phase  and 
the  convoluted  meandering  channels  of  the  "sponge"  or  L3 
phase.  The  topological  distinctions  between  these  two 
phases  result  in  very  different  macroscopic  properties. 
The  anisotropic  lamellar  phases  are  birefringent  and 
viscous,  with  complicated  flow  responses,  while  the 
isotropic  sponge  phases  are  clear  and  generally  flow  freely 
with  constant  viscosity.  However,  these  very  distinct 
phases  are  adjacent  in  the  phase  diagrams  of  many  sys- 
tems, and  very  dilute  sponge  phases  exhibit  a  transient 
birefringence  when  shaken,  which  suggests  a  shear- 
induced  transition  to  a  lamellar  state.  Due  to  the  experi- 
mental challenges  of  studying  both  the  structure  and 
rheology  of  such  dilute  phases,  the  exact  nature  of  the 
transient  birefringence  has  remained  somewhat  controver- 
sial. 

Both  sponge  and  lamellar  structures  have  scaling 
properties:  their  characteristic  lengths  (see  Fig.  1),  the 
channel  width,  ^,  and  the  smectic  periodicity,  d,  scale  with 
the  ratio  of  the  membrane  thickness,  5,  to  the  membrane 
volume  fraction,  (j).  The  dynamics  of  these  phases  are 
dominated  by  thermodynamic  fluctuations  of  the  mem- 
branes, which  also  scale  with  these  parameters.  Gates  and 
Milner  predict  the  scaling  relationship  and  shear  response 
to  be  expected  for  an  isotropic  diblock  copolymer  phase 
[1  ].  They  note  that  this  phase  is  topologically  similar  to 
sponge  phases,  and  therefore  expect  the  relationship  to 
extend  to  such  systems  as  well.  Their  analysis  indicates 
that  flow  should  dampen  fluctuations  normal  to  both  the 
flow  velocity  ( V)  and  the  velocity  gradient  (  W),  eliminat- 
ing channels  in  those  directions.  This  would  result  in  a 
quasi-first  order  phase  transition  to  a  lamellar  state  with  its 
membranes  parallel  to  the  {V,  W)  plane,  i.e.,  perpendicular 
to  Z,  (conventionally  the  "a"  orientation)  occurring  at  a 
critical  shear  rate:      ~  k^T/t]^^  ~  k^T(^/S//rj^^  where 

is  the  Boltzmann  constant,  T  is  the  temperature,  and  tj^ 
is  the  solvent  viscosity. 


From  the  equation,  it  is  clear  that  fransient  birefrin- 
gence, occurring  at  low  shear  rates,  requires  very  dilute 
samples,  i.e.,  having  small  ^.  Furthermore,  the  ^^/tj^ 
dependence  suggests  that  increasing  solvent  viscosity  to 
slow  membrane  fluctuations  would  relax  this  requirement. 
Taking  our  lead  from  previous  microemulsion  work  [2], 
we  found  that  dextrose  could  be  dissolved  at  volume 
fractions  up  to  40  %  in  the  heavy  brine  solvent  (0.2  mol  of 
NaCl  in  D2O  to  make  a  1  L  solution)  of  cetylpyridinium 
chloride/hexanol  membrane  systems  without  altering  their 
phase  diagrams.  The  maximum  dextrose  addition  changes 
^5  from  1.1  mPa-s  to  13.6  mPa-s.  Then  for  a  series  of 
sponge  samples  having  membrane  volume  fractions  from 


Fig.  1.  Representations  of  two  types  of  membrane  morphologies:  the 
sponge  (L3)  and  lamellar  (L  J  phases  with  their  corresponding 
characteristic  lengths. 


3  %  to  10  %,  the  factor  (/>^/%  can  be  varied  by  more  than 
two  orders  of  magnitude.  Plots  of  the  reduced  viscosity 
and  shear  stress  (Fig.  2b)  versus  the  rescaled  shear  rate 
parameter,  yrjj^^,  display  master  curve  behavior,  clearly 
showing  that  these  macroscopic  responses  measuring  the 
dynamics  of  this  system  scale  as  predicted.  The  corre- 
sponding scaled  micro-structural  response  parameters, 
explored  by  Couette  SANS  (Fig.  2c,  d)  likewise  exhibit 
excellent  scaling  behavior.  Moreover,  the  strong  correlation 
between  these  macroscopic  and  microscopic  responses 
leads  to  a  detailed  picture  of  the  mechanisms  involved. 

The  response  exhibits  three  distinct  regions  [3].  (See 
Figs.  2  and  3.)  At  low  shear  rates  (Region  I),  the  sponges 
show  their  familiar  constant  viscosity  Newtonian  behavior 
and  isotropic  scattering.  No  structural  changes  occur 
since  the  applied  shear  rate  is  much  less  than  the  dynamic 
response  rate  of  the  channels.  In  region  II,  when  the  shear 
rate  approaches  this  rate,  the  density  of  connections  in  the 
gradient  direction  (  W)  decreases,  leading  to  a  gradual 
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Fig.  2.  (a)  Rheological  data  for  a  range  of  membrane  volume  fractions  ((zij  and  dextrose  concentrations  (t// 
(b)  The  same  plotted  to  show  master  curve  behavior,  indicative  of  scaling,  (c)  The  corresponding  master 
curve  variation  of  the  scattering  anisotropy,  A^f^f^  =  (Iz -I  rv)/(/z  +  '  rv).  and  (d)  the  rescaled  correlation 
peak  position  and  intensity.  Regions  I,  II  and  III  are  discussed  in  the  text. 


development  of  smectic  order  in  that  direction,  and  thus 
causing  the  observed  decrease  in  viscosity.  At  shear  rates 
high  enough  to  totally  eliminate  the  channels  along  FFthe 
intensity  I^-  saturates,  indicating  the  appearance  of  a 
lamellar  phase  with  the  membranes  fully-aligned  parallel  to 
Z  In  view  of  the  theory,  this  is  an  unexpected  result:  the 
order  develops  in  the  so-called  "c" 
orientation,  indicating  that  it  is  fluctuations 
in  the  gradient  direction  (  W)  that  are 
being  suppressed  rather  than  in  the 
vorticity  (Z)  direction  of  the  predicted  "a" 
orientation.  Also  the  complete  lack  of  a 
rheological  stress  plateau  (Fig.  2a)  and  the 
gradual  change  of  the  peak  position 
(Fig.  2d)  both  preclude  the  possibility  of 
the  transition  being  first  order.  The 
aligned  lamellar  phase  remains  stable 


over  only  a  narrow  range  of 
shear  rates  in  region  II  before 
collapsing  at  the  onset  of 
region  III  to  a  structure  with 
correlations  on  length  scales 
»     which  is  only  observed 
for  samples  with  sufficiently 
high  solvent  viscosities  and  low 
membrane  concentrations. 
While  more  work  remains  to 
determine  the  nature  of  this 
final  structure,  we  note  that  it 
has  been  predicted  that  sup- 
pression of  fluctuations  by 
shear  should  eventually  destabi- 
lize a  bulk  lamellar  system  [4  ]. 

While  the  magnitude  of 
the  critical  shear  rates  are 
consistent  with  theory  and  we 
have  now  confirmed  the 
scaling  behavior  of  the  shear 
response,  the  nature  of  the 
transition  and  the  orientation  of  the  shear-induced  smectic 
state  are  not  as  predicted.  The  theory  does  not,  in  fact, 
extend  straightforwardly  to  the  topologically  similar  sponge 
phases  as  expected,  and  a  more  sophisticated  model  is 
needed  to  fully  explain  the  details  of  shear  response  in  this 
class  of  systems. 
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Fig.  3.  (a)  Schematic  of  the  beam  incidence  on  the  Couette  SANS  cell  in  radial  and  tangential 
geometries,  (b)  Scattering  patterns  in  radial  and  tangential  geometries  for  a  =  0.05,  ii/  =  0.4 
sponge  sample  at  a  range  of  shear  rates  over  the  three  rheological  regions. 
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Observation  of  Liquid-Glass-Liquid 
and  Glass-Glass  Transitions  in  the 
L64/D2O  Copolymer  Micellar  System 


'he  glassy  state  of  matter  is 
ubiquitous,  ranging  from  ordinary 
window  glass  to  polymers  to  proteins.  Yet,  the 
nature  of  the  transition  between  liquid  (L)  and 
glassy  (G)  states  remains  one  of  the  unsolved 
problems  in  the  physics  of  condensed  matter.  Mode- 
coupling  theory  (MCT),  originally  developed  for 
supercooled  simple  liquids,  makes  predictions  of  liquid- 
glass  (L-G)  and  glass-glass  (G-G)  transitions  in  a  colloidal 
suspension.  Here  we  present  small-angle  neutron  scattering 
(SANS)  and  photon  correlation  spectroscopy  (PCS) 
measurements  confirming  these  predictions. 

In  a  recent  series  of  papers  [1,2]  MCT  predicts  that 
in  colloidal  systems  with  a  sufficiently  short-range  effec- 
tive inter-particle  attractive  interaction  one  can  observe  two 
types  of  kinetic  glass  transitions  (KGT).  The  KGT  can  be 
the  result  of  a  "cage  effect,"  originating  trom  the  excluded 
volume  effect  of  the  particles  at  high  volume  fractions  and 
at  high  temperatures,  or  the  result  of  a  particle  bonding  to 
its  nearest  neighbors  at  low  k^T  and  at  all  volume  frac- 
tions. 

One  particularly  interesting  feature  in  this  novel  state 
of  matter  is  the  occurrence  of  two  distinct  amorphous 
glassy  states  called  respectively  the  "repulsive  glass"  (Gj^) 
and  the  "attractive  glass"  (G^),  as  the  result  of  the  two 
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Fig.  1.  The  lines  represent  phase  boundaries  (in  T*-  plane) 
determined  by  mode-coupling  theory  for  a  colloidal  system  with  a  short- 
range  attraction  (s  =  0.03).  The  phase  points  shown  are  for  the  L64/D2O 
micellar  solution.  (Liquid  state:  orange,  glassy  states:  blue,  red,  and 
green.)  The  results  confirm  the  existence  of  the  attractive  branch  in  the 
predicted  glass-transition  boundary,  and  the  G^-Gr  transition. 


distinct  mechanisms  for  the  structural  arrest  mentioned 
earlier.  With  PCS  and  SANS  measurements,  we  succeeded 
in  detecting  both  these  transitions  in  L64/D2O  micellar 
system  at  high  volume  fractions  [3]. 

L64  is  a  symmetric  tri-block  copolymer  of  composi- 
tion (PEO)]3(PPO)3o(PEO)i3.  Above  20  °C,  the  copoly- 
mers self-assemble  into  spherical  micelles  each  with  PPO 
segments  packed  into  a  compact  core  and  PEO  segments 
forming  a  corona  around  it.  The  short-range  attraction 
arises  from  the  overlap  of  PEO  chains  in  the  corona 
regions  when  two  micelles  approach  one  another.  The 
depth  of  the  attractive  well  increases  with  temperature 
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Fig.  2.  The  ISF  as  measured  by  photon  correlation  spectroscopy  at 
three  different  volume  fractions,  cp,  at  different  temperatures.  (Data 
follow  the  same  color  scheme  as  in  Fig.  1 .)  The  top  panel  shows  the  re- 
entrant L-G^-L  transition.  The  line  highlights  a  region  of  logarithmic 
time  dependence  preceding  the  plateau  in  the  liquid  state  just  before  the 
transition  [4].  Evidence  of  the  G^-Gp  transition  can  be  seen  in  the  middle 
panel.  By  comparing  the  long-time  limit  of  the  ISF's  with  the  predictions 
of  mode-coupling  theory  (MCT),  the  two  different  types  of  the  glasses 
can  be  identified  by  their  respective  Debye-Waller  factors.  The  bottom 
panel  is  for  a  volume  fraction  of  0.544,  which  is  the    end  point 
predicted  by  MCT.  f,^  and    become  identical,  confirming  the 
prediction. 
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because  water  becomes  an  increasingly  poor  solvent  as 
temperature  rises. 

If  the  interaction  is  modeled  by  a  square-well  poten- 
tial, the  crucial  control  parameters  of  the  phase  behavior 
are:  the  volume  fraction,  cp  =  7ipd'^/6,  the  effective  tem- 
perature,     =  k^T/u,  and  the  fractional  well- width 
parameter,  e  -  A/d.  p  is  the  number  density  of  the  spheres, 
d  its  hard  core  diameter,  A  its  attractive  well-width,  and  -u 
its  well-depth. 

The  calculated  glass-fransition  lines,  for  the  case  of 
£  =  0.03  [2],  are  shown  in  Fig.l  in  the  {T*,(p)  plane.  On 
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Fig.  3.  SANS  scaled  intensity  distributions  vs.  reduced  scattering  vector 
for  the  micellar  solutions  at  four  values  of  (p,  at  different  temperatures. 
(Data  follow/  the  same  color  scheme  as  in  Fig.  1 .)  The  top  panel  is  for 
(p=  0.460.  All  the  scaled  intensities  are  characterized  by  a  unique  length 
scale  and  collapse  into  a  single  master  curve  for  all  temperatures.  The 
second  panel  is  for  ^=0.535:  the  system  is  clearly  characterized  by  two 
length  scales.  While  the  narrower  peak  is  resolution  limited,  the  broader 
peak  is  similar  to  the  one  observed  in  the  first  panel.  If  we  identify  the 
broader  peak  to  represent  the  liquid  state,  then  the  narrower  peak  should 
represent  the  glassy  state.  The  system  thus  shows  a  re-entrant  L-G-L 
transition.  The  third  panel  is  for  (p=  0.542.  Again,  two  length  scales 
characterize  the  system.  The  narrower  peak  is  resolution-limited  and  the 
slightly  broader  peak  is  also  nearly  so.  The  results  can  be  interpreted  as 
showing  a  transition  between  two  amorphous  solid  states,  i.e.,  a  G-G 
transition.  The  last  panel  is  for  ^3=  0.544,  which  is  predicted  as  the  A3 
point,  at  temperatures  ranging  from  15  °C  to  50  °C.  The  two  different 
types  of  glasses  indeed  become  identical  at  this  volume  fraction. 


going  from  higher  to  lower  T*,  at  volume  fractions  less 
than  0.5365  L-G^  transitions  occur,  and  at  volume  frac- 
tions greater  than  this  value,  one  has  a  possibility  of 
observing  Gr-G^  transitions. 

Figure  2  shows  normalized  photon  correlation 
functions  (coherent  intermediate  scattering  frinctions,  ISF) 
measured  at  q  »  0.001  A-'  and  at  three  different  volume 
fractions,  as  a  function  of  temperature  [3].  For  the  (p 
where  the  L-G^  transition  is  possible,  one  can  see  that  in 
the  L  state,  the  long-time  limit  of  the  ISF,^,  is  zero,  while 
in  the  G^  state, /^'^  »  0.4  to  0.5,  which  is  the  value  of  the 
Debye-Waller  factor  of  the  glassy  state.  The  KGT  is  thus 
characterized  by  a  discontinuous  change  of  f^.  At 
<p  =  0.542,  we  observe  the  G^  (f^''  »  0.4)  to  G^  transition 
(/^^  «  .7).  At  (p=  0.544,  which  is  called  the  singularity 
end  point,  our  measured  ISF's  show  identical  values  of 
Z^*^"  «  0.4  in  the  two  glassy  states,  confirming  the  theory. 

Figure  3  shows  scaling  plots  of  SANS  intensity 
disfributions  at  four  volume  fractions,  and  at  different 
temperatures,  spanning  the  regions  of  the  phase  diagram 
where  the  L,  L-G-L,  G-G,  and  two  identical  glasses 
{A2  point)  are  located.  The  scale  is  set  by  a  unique  length 
parameter,  the  average  inter-particle  distance. 

Verification  of  mode-coupling  theory  predictions 
through  photon-correlation  spectroscopy  and  small-angle 
neutron  scattering  measurements  on  a  micellar  system 
provides  a  basis  for  improved  understanding  of  the  glass 
transition  in  general. 
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Early-Stage  Compositional 
Segregation  in  Polymer  Blend  Films 


etting  plays  an  important  role  in 
diverse  applications  ranging  from 
insecticide  sprays  to  detergency.  This  interest, 
in  part,  has  motivated  a  great  deal  of  fundamental 
research  involving  the  spreading  (wetting)  of  a 
liquid  on  a  solid  surface.  Polymer  blends  are  of  great 
commercial  importance  and  therefore  their  thermody- 
namic and  phase  separation  behaviors  have  been  widely 
studied. 

The  structure  and  properties  of  thin  blend  films  are 
governed  by  the  interplay  between  phase  separation  and 
surface  segregation  driven  by  polymer-surface  interac- 
tions. The  nature  of  the  segregation  depends  on  whether 
the  mixture  forms  a  stable  one-phase  fluid  or  phase- 
separates  [1,  2].  In  the  former  case,  the  surface-enrich- 
ment-layer (SEL)  grows  and  stabilizes  to  an  equilibrium 
value.  In  the  latter  case  the  initially  formed  SEL  breaks  into 
droplets  for  partial  wetting  and  grows  unsteadily  for 
complete  wetting  [1].  The  late  stage  evolution  of  the 
surface  enrichment  layers  thus  exhibits  "non-universal" 
growth  laws,  depending  on  the  phase  stability,  wetting 
characteristics  of  the  blend,  confinement,  as  well  as  on  the 
details  of  the  polymer-surface  interaction  potential  [2].  The 
universality  of  surface  segregation  at  the  early-stage, 
however,  has  not  been  explored,  partly  because  the 
compositional  changes  are  normally  rapid  [3,  4].  We 
therefore  designed  a  measurement  that  can  access  the 
short  time  and  small  length-scale  surface-enrichment 
regime. 

A  critical  blend  film  near  the  glass  transition  of  the 
polymers  is  utilized  so  that  the  chain  dynamics  are  greatly 
slowed  down  even  under  thermodynamically  strong 
segregation  conditions.  The  polymers  are  deuterated 
poly(methyl  methacrylate)  (dPMMA)  and  poly(styrene- 
ra^-acrylonitrile)  (SAN).  This  polymer  blend  has  a  lower 
critical  solution  temperature  (LCST)  below     as  revealed 
by  small-angle  neutron  scattering  measurements.  A  445  A 
film  of  dPMMA/SAN  (50/50  by  mass)  was  obtained  by 
spin-casting  fi"om  a  common  solution.  After  drying  at  a 
temperature  slightly  below  T^,  the  films  were  annealed  in  a 
vacuum  oven  at  (130  ±  0.1)  °C  for  various  times  ranging 
from  1  min  to  36  h.  After  drying  and  after  each  step  of 
annealing,  the  film  was  measured  using  neutron  reflectivity 


(NR).  An  atomic  force  microscopy  study  indicated  that  the 
film  roughness  was  below  6  A  throughout  the  measure- 
ment. 

Figure  1  shows  NR  spectra  after  several  annealing 
intervals.  The  spectra  were  analyzed  using  a  model-fitting 
scheme.  The  best  fits  are  shown  as  the  solid  curves,  and 
their  corresponding  compositional  profiles  are  shown  in 
Fig.  2.  The  inset  of  Fig.  2  shows  the  evolution  of  the 
dPMMA  volume  fraction  in  the  depletion  zone  (^^j)  and  at 
the  surface  {(f>^.  The  former  decreases  monotonically, 
suggesting  that  equilibrium  has  not  been  achieved,  while 
the  non-monotonic  variation  of  the  latter  is  due  to  the 
surface  relaxation  at  the  initial  time. 

The  surface  excess,  z*  =  J[<z^(x)  -  ^d\dx,  is  shown 
in  Fig.  3,  and  can  be  fit  to  a  power  law  f0.4i±o.oi  for  early 
times  (solid  line),  and  a  stretched  exponential  form  for  the 
entire  time  (dashed  curve).  The  latter  gives  z*(oo)  =  72  A 
at  the  saturation,  z*  starts  to  deviate  from  the  power-law 
behavior  at  about  9  h  (indicated  by  the  arrow  in  Fig.  3) 
corresponding  to  z*  =  40  A  and  4  =  57  A,  where  4  is  the 
surface  correlation  length  characterizing  the  decay  of  the 
surface  enrichment  composition.  These  crossover  scales 
in  the  kinetics  are  comparable  to  the  bulk  correlation  length 
length  4,  ~  60  A,  suggesting  a  transient  saturation.  After 
that,  the  SEL  will  grow  unsteadily,  entering  the  later  stage 
of  the  wetting  layer  growth  [4]. 

The  inset  of  Fig.  3  shows  that  ^  nearly  follows  a 
power  law  of  /0  25±o.oi^  except  for  initial  times.  Note  that 
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Fig.  1.  Selected  neutron  reflectivity  spectra  after  various  times  at  130  "C. 
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Fig.  2.  (z^tjPMMA  profiles  that  give  the  best  fit.  The  inset  shows  the 
evolution  of  surface  composition,  (circles)  and  depletion  layer 
composition,  (p^  (triangles). 
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Fig.  3.  z*  with  a  power-law  fit  (solid  line),  f  ''^,  and  a  stretched 
exponential  function  (dashed  curve).  The  arrow  denotes  a  z*  ^  40  A  and 
the  transition  to  slower  growth.  The  inset  shows  the  surface  correlation 
length. 


right  after  drying,  z*(0)  =  4  A,  indicating  a  segmental  scale 
of  saturation  from  processing  the  film  (spin-casting  and 
drying).  Thus,  the  initial  "anomaly"  of  (j)^  and  may 
reflect  the  relaxation  of  the  highly  non-equilibrium  compo- 
sitional profiles.  While  most  of  the  previous  work  focused 
on  the  kinetics  after  the  formation  of  the  initial  layer  on  the 
order  of  the  correlation  length,  the  transient  early-stage  of 
the  SEL  has  received  little  attention.  The  interpretation  of 
the  apparent  power  law  growth  will  thus  require  further 
investigation. 

In  summary,  during  the  early  stage  of  surface 
segregation,  enrichment  at  a  monomer  length-scale  occurs 
first,  followed  by  a  segregation  layer  that  grows  to  the 
order  of  the  bulk  correlation  length,  4^  which  coincides 
with  the  radius  of  gyration  of  polymers  deep  in  the  two- 
phase  regime.  Thus,  the  two  observed  distinct  surface 
saturation  states  correspond  to  the  characteristic  scales  of 
macromolecules.  However,  at  near  criticality,  previous 
studies  suggest  a  large  discrepancy  between  4  and 
pointing  to  the  need  to  consider  three  length  scales:  that  of 
the  monomer,  the  molecular  size,  and  the  size  of  composi- 
tional fluctuations.  The  extent  of  universality  in  the  early 
stage  growth  kinetics  therefore  needs  further  investigation. 
The  present  study  is  the  first  to  quantitatively  examine  this 
important  phenomenon. 
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Magnetic  Cluster  Sizes 
in  iVIagnetic  Recording  Disks 


■he  equivalent  of  Moore's  law  for 
magnetic  recording  is  that  the  areal 
density  of  magnetically  written  bits  increases  at 
60  %  to  100  %  per  year,  and,  of  equal  interest  to 
the  consumer,  the  cost  per  megabyte  decreases  at 
40  %  to  50  %  per  year.  To  achieve  this,  the  noise  in 
the  magnetic  media  must  be  continuously  reduced, 
which  is  accomplished  by  decreasing  the  size  of  the 
magnetic  clusters  —  regions  of  the  media  that  are  mag- 
netically coupled.  The  cluster  size  must  decrease,  since  the 
media  noise  primarily  results  from  the  finite  cluster  size; 
the  bit  transitions  are  narrower  when  the  clusters  are 
smaller.  Presently,  cluster  sizes  are  believed  to  be  about 
10  nm  to  15  nm,  but  accurate  knowledge  of  the  size 
distribution  and  even  the  average  size  is  lacking.  A  com- 
mon assumption  is  that  the  cluster  size  is  identical  to  the 
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Fig,  1.  SANS  intensity  for  magnetic  disl^s  So  and  Cy,  averaged  around  tlie  incident  beam  direction  whiichi 
was  perpendicular  to  surface  of  the  disks.  In  (a)  and  (b)  the  blue  lines  show  data  for  no  applied  field  (high 
noise  state),  while  red  lines  show  data  for  a  0.6  T  field.  Two  sample-to-detector  distances  were  used,  with 
Q  ranges  of  0.0025  A-i  to  0.025  A-^  and  0.025  A-^  to  0.3  A-^ .  In  (c)  and  (d)  the  blue  squares  are  the  magnetic 
SANS,  the  difference  of  the  zero  and  0.6  T  scattering  intensities  in  (a)  and  (b).  The  red  line  shows  the  best 
fit  to  a  model  described  in  the  text.  The  arrows  mark  the  approximate  Q  where  the  slope  of  the  data 
changes. 


crystalline  grain  size  of  the  media  [1],  but  this  has  not  been 
demonstrated. 

We  have  conducted  Small  Angle  Neutron  Scattering 
(SANS)  measurements  of  the  magnetic  cluster  size  in 
actual  recording  disks  and  quantified  the  size  distribution. 
The  measurements  were  performed  at  NG-7  30  m  SANS 
instrument  with  unpolarized  neutrons  with  wavelength  of 
7  A  and  a  relative  wavelength  spread  of  1 1  %.  Since  SANS 
is  caused  by  both  physical  and  magnetic  inhomogeneities, 
it  is  necessary  to  separate  these  contributions  and  to  isolate 
the  magnetic  cluster  SANS.  This  was  done  with  two 
separate  measurements.  Data  were  first  taken  for  disks  in 
an  initial,  high  noise  state  where  there  are  a  maximum 
number  of  interfaces  between  magnetic  clusters.  This  is 
accomplished  by  first  fully  magnetizing  the  disks  in  a 
circumferential  direction.  Then,  the  disks  were  scanned 

with  a  recording  head  field  in 
the  opposite  direction  and  with 
a  strength  approximately  equal 
to  the  remnant  coercivity  of 
the  media.  This  causes  ap- 
proximately half  of  the  mag- 
netic clusters  to  reverse 
direction,  leaving  a  maximally 
noisy  state.  In  this  state,  the 
SANS  contains  contributions 
from  the  physical  film  struc- 
ture (grains)  and  the  magnetic 
film  structure  (magnetic 
clusters).  Data  were  then 
obtained  with  a  magnetic  field 
of  0.6  T  applied  parallel  to  the 
disk  surface.  This  field  orients 
the  magnetic  moments  in  the 
clusters  predominantly  along 
the  field  and,  by  examining  the 
^-dependence  of  the  SANS, 
we  can  determine  the  SANS 
from  the  physical  film  struc- 
ture. Subtracting  the  zero  field 
and  0.6  T  spectra  leaves  the 
desired  magnetic  SANS  [2,  3]. 
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The  disks,  denoted  So,  To  and  Cy,  are  those  used  in 
three  product  generations.  The  magnetic  media  in  each 
two-sided  disk  was  a  Co-Pt-Cr  alloy  (25  nm  to  40  nm 
thick)  sandwiched  between  a  carbon  layer  and  one  or  more 
nonmagnetic  underlayers.  The  disks  were  sliced  into  eight 
pieces  and  stacked  to  create  1 6  magnetic  layers  for  the 
SANS  measurements.  The  relative  signal-to-noise  ratio 
(Sq/NJ  varied  from  0  db  (So)  to  3  db  (To)  to  6.5  db  (Cy). 

SANS  data  for  both  magnetic  fields  are  shown  in 
Fig.  la  and  lb.  Subtracting  these  two  spectra  leaves  the 
desired  magnetic  SANS,  which  are  shown  in  Fig.  Ic  and 
Id.  It  is  apparent  that  for  increasing  Sq/N^^  (So  to  Cy) 
there  is  an  increase  in  the  Q  where  the  slope  changes, 
which  establishes  that  the  disks  with  higher  Sq/N^  have 
smaller  average  cluster  size.  In  addition,  there  is  a  distribu- 
tion of  cluster  sizes,  since  there  are  no  oscillations  in  the 
data  and  the  change  in  slope  is  gradual. 

To  determine  the  cluster  size  distribution,  we  have 
modeled  the  cluster  shapes  as  cylinders,  since  TEM  shows 
the  physical  grains  adopt  a  columnar  morphology.  The 
magnetic  SANS  data  are  well  fitted  using  a  log-normal 
distribution  of  cluster  sizes.  These  distributions  are  shown 
in  Fig.  2,  while  the  fits  to  the  data  are  shown  in  Fig.  1. 
Such  magnetic  cluster  size  distributions  have  not  been 
reported  before  for  recording  media,  and  allow  us  to  draw 
several  important  conclusions. 

First,  these  results  show  that  the  media  Sq/N^  is 
inversely  correlated  with  the  magnetic  cluster  size.  This  is 
apparent  from  Fig.  2  when  one  recalls  that  So/N,„  increases 
from  So-To-Cy.  This  provides  firm  experimental 
verification  of  the  models  that  have  proposed  smaller 
clusters  give  larger  Sg/N^  [!]•  A  second  conclusion 
follows  from  a  comparison  of  the  magnetic  cluster  size  to 
the  physical  grain  size,  shown  in  the  inset  to  Fig.  2.  The 
average  cluster  size  is  slightly  larger  than  the  average  grain 
size,  and,  for  the  most  advanced  media  (Cy),  this  differ- 
ence is  <  10  %.  This  shows  that  the  intergranular  magnetic 
interactions  are  not  very  strong,  especially  for  the  ad- 
vanced media;  hence,  the  size  of  the  media  noise  sources 
scales  with  crystalline  grain  size.  This  result  validates 
previous  (implicit  or  explicit)  assumptions  [1,4]. 

The  cluster  size  distribution  is  quite  broad  with  a  log- 
normal  r  of  about  0.5.  Recent  models  of  the  magnetic 
recording  properties  of  media  have  begun  to  include 
microstructural  disorder  [1],  and  have  shown  that  the 
magnetic  cluster  size  distribution  has  significant  effects  on 
these  properties.  Thus,  the  exact  cluster  size  distribution  is 
important  for  these  models.  Our  results  provide  empirical 
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Fig.  2.  Green,  red  and  blue  curves  show  the  best-fit  distributions  for 
disks  Cy,  To  and  So,  respectively.  The  triangles  mark  the  average  cluster 
size.  The  log-normal  rare  0.56,  0.52  and  0.46  for  Cy,  To  and  So, 
respectively.  The  inset  shovi/s  the  average  magnetic  cluster  diameter 
versus  the  average  physical  grain  diameter,  obtained  from  x-ray 
diffraction  and  TEM. 

data  on  this  distribution,  which,  when  incorporated  into  the 
models,  will  improve  their  reliability. 

Our  results  have  implications  for  magnetic  recording 
technologies.  While  the  lack  of  coupling  is  desirable  from 
an  Sy/N^,  perspective,  the  broad  distribufion  is  not,  since 
large  clusters  will  contribute  to  noise  and  the  small  clusters 
will  be  thermally  unstable.  Reducing  this  distribution  will  be 
a  major  challenge  for  media  designers.  Finally,  the  method- 
ology we  have  developed  will  be  applicable  to  future  media 
where  the  average  cluster  sizes  are  <  10  nm. 
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Dynamical  Effects  of  Polar 
Nanoregions  in  Relaxor  Ferroelectrics 


I  elaxor  ferroelectrics  exhibit 
^exceptional  piezoelectric  properties. 
Single  crystals  of  some  of  these  materials  can 
achieve  ultrahigh  strains  under  applied  electric 
fields  that  are  fially  one  order  of  magnitude  larger 
than  those  attainable  in  conventional  lead-zirconate- 
titanate  (PZT)  ceramics,  which  are  widely  used  in 
solid-state  actuators  that  convert  electrical  energy  into 
mechanical  energy.  The  additional  properties  of  large 
electromechanical  coupling  and  low  dielectric  loss  suggest 
that  industrial  use  of  relaxors  may  revolutionize  the  areas 
of  medical  imaging,  naval  sonar,  and  other  acoustic 
applications.  Relaxors  are  further  characterized  by  a 
markedly  fi-equency-dependent  dielectric  susceptibility  that 
is  very  broad  in  temperature.  The  recent  neutron  scatter- 
ing studies  briefly  described  here  provide  new  insight  on 
the  behavior  of  these  systems. 

PZN 


> 


3 


0.2  0.3 

(2,  k,  0) 

Fig.  1.  Contour  plot  of  the  neutron  inelastic  scattering  intensity 
measured  in  PZN  in  its  cubic  phase  in  the  (200)  Brillouin  zone.  The 
color  scale  varies  logarithmically  with  intensity,  and  is  limited  in  range 
to  make  the  weaker  features  more  visible. 


A  key  feature  that  appears  to  be  common  to  all  of  the 
relaxor  compounds,  and  which  is  believed  to  play  a  funda- 
mental role  in  producing  the  enhanced  piezoelectricity,  is 
the  formation  at  high  temperature  of  tiny  regions  of  local 
and  randomly-oriented  electric  polarization  just  several  unit 
cells  in  size,  also  known  as  polar  nanoregions  (PNR).  The 
existence  of  these  PNR  was  inferred  from  the  observation 
that  the  optical  index  of  refraction  in  a  variety  of  disordered 
ferroelectric  systems  deviates  from  a  linear  temperature 
dependence  at  a  temperature     far  above  [2]. 

Pb(Zn,/3Nb2/3)03  (PZN)  and  Pb(Mgi/3Nb2/3)03 
(PMN)  are  two  of  the  most  studied  compounds  as  their 
piezoelectric  responses  rank  among  the  highest  [1].  These 
compounds  crystallize  in  the  complex  perovskite  structure 
Pb(B'B")03,  where  the  B'  and  B"  atoms  are  generally 
heterovalent.  Thus  the  relaxor  compounds  possess  not 
only  chemical  disorder,  but  usually  a  mixed-valence 
character  as  well. 

Our  studies  of  PZN  and  PMN  have  focused  on 
understanding  the  mechanism  underlying  the  formation  of 
the  PNR,  as  well  as  their  effect  on  the  lattice  dynamics. 
The  latter  are  clearly  manifest  in  neutron  scattering 
measurements.  The  long-wavelength  transverse  optic 
(TO)  phonons  are  heavily  damped  in  energy  below  in 
PZN,  which  gives  rise  to  the  anomalous  "waterfall"  feature 
shown  in  Fig.  1  [3].  This  damping  is  strongly  wavevector 
dependent,  and  becomes  significant  when  the  wavelength 
of  the  lattice  vibration  becomes  comparable  to  the  average 
size  of  the  PNR,  which  is  estimated  to  be  of  the  order  of 
30  A  to  40  A. 

Subsequent  experiments  were  performed  above 
=  620  K  in  PMN  to  look  for  a  soft  mode,  given  the 
obvious  similarities  to  the  classic  soft-mode  ferroelectric 
PbTi03.  Surprisingly,  the  phonon  damping  persists  well 
above  T^,  as  shown  in  Fig.  2  [4].  This  implies  that  the 
PNR  exist  as  dynamic  entities  above  T^,  but  are  still 
effective  at  inhibiting  the  propagation  of  long  wavelength 
phonon  modes.  Even  more  interesting  is  how  the  zone- 
center  TO  mode  softens  and  broadens  in  energy  as  the 
temperature  approaches     from  above,  and  that  only  the 
zone-center  TO  mode  becomes  overdamped  (no  peak  at 
non-zero  energy)  at  T^.  Modes  at  higher  wavevector  q 
become  overdamped  too,  but  only  below  T^.  This  behav- 
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ior  suggests  that  the  zone-center  TO  mode  condenses  at 
r^j,  thereby  giving  rise  to  the  PNR.  Consistent  with  this 
picture  is  the  observation  that  the  square  of  the  zone-center 
TO  mode  energy  varies  Hnearly  with  temperature,  as 
happens  in  normal  soft-mode  ferroelectrics.  This  is  shown 
in  Fig.  3,  which  summarizes  the  soft  TO  mode  dynamics 
in  PMN  between  1 0  K  and  1 1 00  K  [5] . 

An  anomalous  broadening  of  the  acoustic  modes, 
shown  in  the  upper  panel  of  Fig.  3,  also  occurs  at  T^,  as 
does  the  onset  of  diffuse  scattering  associated  with  the 
PNR.  The  difftise  scattering  structure  factors  can  be 
reconciled  with  those  of  the  TO  phonons  in  different 
Brillouin  zones  by  postulating  a  uniform  shift  of  the  PNR 
along  their  polarization  direction  [6].  Furthermore,  a  model 
that  couples  the  TA  and  TO  modes  is  able  to  describe  all  of 
the  observed  phonon  lineshapes  extremely  well  [7].  This 
then  suggests  a  novel  idea  in  which  PMN  exhibits  the 
condensation  of  a  coupled  soft  TO  mode  at  T^.  Such  a 
coupled-mode  would  carry  an  acoustic  component,  and 
thus  provide  a  natural  explanation  of  the  uniform  shift  of 
the  PNR.  If  true,  this  would  complete  an  elegant  descrip- 
tion for  the  dynamical  formation  of  the  PNR  in  PMN,  and 
most  likely  other  relaxors  as  well. 
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Emergent  Excitations  in  a 
Geometrically  Frustrated  Magnet 


A "Chinese"  checkerboard  having  an 
equal  number  of  black  and  white 
marbles  cannot  have  every  neighbor  being  the 
opposite  kind,  whereas  a  square  checkerboard 
can.  The  triangular  array  on  a  Chinese  checker- 
board offers  an  example  of  "geometrical  finstration". 
Frustration  (i.e.,  a  condition  in  which  some  relation 
among  the  components  cannot  be  satisfied  simulta- 
neously) is  an  important  concept  in  understanding  disor- 
dered states  of  matter.  Systems  in  which  a  large  diversity 
of  states  are  involved  are  common  in  biology,  chemistry 
and  physics  [1,2].  Notable  examples  are  glasses,  liquids 
and  proteins.  Magnetic  systems  offer  extreme  examples  of 
frustration  in  the  form  of  spin  lattices,  where  all  interac- 
tions between  spins  cannot  be  simultaneously  satisfied. 
Such  geometrical  frustration  can  lead  to  the  emergence  of 
qualitatively  new  states  of  matter,  having  "composite" 
degrees  of  freedom. 

To  explore  this  possibility,  we  examined  magnetic 
fluctuations  in  ZnCr204  [3].  The  B-site  of  this  spinel  lattice 
occupied  by  spin-3/2  Cr^+  leads  to  a  magnet  with  dominant 
nearest  neighbor  interactions  on  the  lattice  of  comer- 
sharing  tetrahedra  [4]  shown  in  Fig.  1 .  Because  the  spin 
interaction  energy  is  minimized  when  the  four  spins  on 
each  tetrahedron  add  to  zero,  interactions  do  not  call  for 
long-range  order,  but  simply  define  a  restricted  phase 
space  for  fluctuations.  Just  as  composite  fermions  can 
emerge  from  degenerate  Landau  levels  in  a  two-dimen- 
sional electron  gas,  the  near-degenerate  manifold  of  states 
in  a  frustrated  magnet  is  fertile  ground  for  emergent 
behavior  [5]. 

Neufron  scattering  provides  the  most  effective  tool 
to  study  possible  composite  spin  degrees  of  freedom  by 
directly  probing  the  form  factor  of  such  entities.  Fig.  2  (a,b) 
show  the  wave  vector  dependence  of  the  low  energy 
inelastic  neutron  scattering  cross  section  in  the  spin-liquid 
phase  of  ZnCr204.  The  data  exhibit  broad  maxima  at  the 
Brillouin  zone  boundaries,  signaling  the  emergence  of 
confined  nano-scale  spin  clusters.  Rather  than  Fourier- 
inverting  the  data,  we  consider  potential  spin  clusters  and 
test  the  corresponding  prediction  for  the  form  factor  against 
the  data. 

Individual  tetrahedra  would  be  prime  candidates  for 
such  clusters,  as  they  constitute  the  basic  motif  of  the 


pyrochlore  lattice.  However,  a  tetrahedron  is  too  small  to 
account  for  the  observed  features.  The  next  largest 
symmetric  structural  unit  is  the  hexagonal  loop  formed  by 
a  cluster  of  six  tetrahedra  (Fig.  3).  Two  spins  of  each 
tetrahedron  occupy  the  vertices  of  a  hexagon  while  the 
other  two  spins  belong  to  different  hexagons.  It  is  possible 
to  assign  all  spins  on  the  spinel  lattice  to  hexagons  simulta- 
neously, thus  producing  N/6  weakly  interacting  degrees  of 
freedom  (Fig.  1 ).  An  outstanding  fit  is  achieved  for  the 
antiferromagnetic  hexagonal  spin  loops,  as  displayed  in 
Fig.  2  (c,  d).  Thus,  rather  than  scattering  from  individual 
spins,  neutrons  scatter  from  antiferromagnetic  hexagonal 
spin  clusters.  In  effect,  ZnCr204  at  low  temperatures  is  not 
a  system  of  strongly  interacting  spins,  but  a  protectorate  of 
weakly  interacting  spin-loop  directors.  Since  the  six 
hexagon  spins  are  anti-parallel  with  each  other,  the  stag- 
gered magnetization  vector  for  a  single  hexagon,  which 
shall  be  called  the  spin  loop  director,  is  decoupled  from  the 
12  outer  spins,  and  hence  its  reorientation  embodies  the 


Fig.  1 .  The  lattice  of  corner-sharing  tetrahedra  formed  by  the 
octahedrally  coordinated  B  sites  in  a  spinel  structure  with  chemical 
formula  AB2O4.  A  periodic  assignment  of  all  spins  in  the  pyrochlore 
lattice  is  made  to  four  different  types  of  non-overlapping  hexagons, 
represented  by  the  colors  blue,  green,  red,  and  gold.  Every  spin  belongs 
to  just  one  hexagon  and  each  such  hexagon  carries  a  six-spin  director. 
The  resulting  tetragonal  structure  of  these  hexagons  has  a  unit  cell  of 
2a  X  2a  ^  3c  and  can  be  described  by  a  stacking  of  two  different  types  of 
three-layer  slabs  along  the  c-axis.  The  hexagon  coverage  on  consecutive 
slabs  is  in  fact  uncorrelated,  so  that  a  macroscopic  number  of  random 
slab-sequences  can  be  generated. 
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long-sought  local  zero-energy  mode  for  the  pyrochlore 
lattice. 

Composite  degrees  of  freedom  are  common  in 
strongly  interacting  many  body  systems.  Quarks  form 
hadrons;  hadrons  form  nuclei;  nuclei  plus  electrons  form 
atoms;  atoms  form  molecules  that  in  turn  are  the  basis  for 
complex  biological  functionality.  Planets,  stars,  galaxies 
and  galactic  clusters  are  examples  of  clustering  on  grander 
length  scales.  However,  to  our  knowledge,  the  emergence 
of  a  confined  spin  cluster  degree  of  freedom  has  not 
previously  been  documented  in  a  uniform  gapless  magnet. 
The  discovery  is  important  because  magnets  offer  an 
opportunity  not  afforded  by  the  aforementioned  systems, 
namely,  to  monitor  emergent  structure  in  complex  interact- 
ing systems  with  microscopic  probes  such  as  neufron 


scattering  and  NMR.  The  collapse  of  a  geometrically 
frustrated  magnet  into  a  director  protectorate  could  for 
example  be  a  useful  template  for  exploring  aspects  of 
protein  folding  [2J. 
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Fig.  2.  (a),  (b):  Color  images  of  inelastic  neutron  scattering  intensities 
from  single  crystals  of  ZnCr204  in  the  (hkO)  and  (hkk)  symmetry  planes 
obtained  at  7  =  15  Kforfify  =1  meV.  The  data  are  a  measure  of  the 
dynamic  form  factor  for  self-organized  nano-scale  spin  clusters  in  the 
material,  (c),  (d):  Color  images  of  the  square  of  the  form  factor 
calculated  for  antiferromagnetic  hexagon  spin  loops  averaged  over  the 
four  hexagon  orientations  in  the  spinel  lattice.  The  excellent  agreement 
between  model  and  data  identifies  the  spin  clusters  as  hexagonal  spin 
loops. 


Fig.  3.  Spin  cluster  surrounding  a  hexagon  (shown  in  gold)  in  the 
pyrochlore  lattice  of  Fig.  1. 
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Polaron  Dynamics 
and  the  Glass  Transition  in 
Colossal  Magnetoresistance  Materials 


The  recent  discovery  of  colossal 
magnetoresistance  (CMR)  effects  in 
the  manganese  oxide  class  of  materials  (such 
as  Laj.^A^MnOj  (A  =  Sr,  Ca,  Ba))  has  rekindled 
intense  interest  in  these  systems,  both  because  of 
their  fundamental  properties  as  well  as  their  techno- 
logical potential  for  read/write  heads  in  magnetic 
recording  media,  sensors,  and  spin-polarized  electronics. 
The  role  of  local  charge  correlations  and  their  competition 
with  magnetic  and  electronic  ground  states  is  one  of  the 
most  striking  new  features  of  the  transition  metal  oxides 
such  as  cuprate  superconductors  and  nickelates,  but  is 
most  elegantly  highlighted  for  these  CMR  manganites.  Our 
neutron  scattering  results  near  optimal  doping  show  that 
the  short-range  polaron  correlations  are  completely 
dynamic  at  high  T,  but  then  freeze  upon  cooling  to  a 
temperature  T„.  This  glass  transition  suggests  that  the 
paramagnetic/insulating  state  arises  from  an  inherent  orbital 
frustration  that  inhibits  the  formation  of  a  long  range 
orbital-  and  charge-ordered  state.  Upon  further  cooling 
into  the  ferromagnetic-metallic  state,  where  the  polarons 
melt,  the  diffuse  scattering  quickly  develops  into  a  propa- 
gating, transverse  optic  phonon. 

In  the  ground  state  these  CMR  materials  are  ideal 
isotropic  ferromagnets.  As  the  system  is  warmed  through 
the  Curie  temperature  of  this  half-metallic  system  the 
disorder  in  the  spins  begins  to  inhibit  the  hopping  of  the 
(magnetic)  conduction  electrons,  and  this  in  turn  provides 
enough  time  for  the  lattice  to  distort  and  form  polarons  that 
trap  the  electrons.  The  ferromagnetic-paramagnetic 
transition  is  thus  accompanied  by  a  metal-insulator  transi- 
tion; cooling  into  the  ferromagnetic  state,  or  alternatively 
inducing  a  magnetization  by  application  of  a  field,  converts 
the  insulator  back  into  a  metal,  producing  the  CMR  effect. 

To  investigate  the  dynamics  of  these  polarons,  we 
have  carried  out  inelastic  neutron  scattering  measurements 
on  a  single  crystal  of  optimally-doped  Lag  jCao  jMnOj 
{Tq  =  252  K)  [1].  Fig.  1  shows  the  energy  dependence  of 
the  polaron  scattering  at  the  peak  of  the  polaron  correla- 
tions [2].  The  data  reveal  that  there  is  a  purely  elastic 
component,  indicating  that  the  system  still  possesses  order, 
and  since  spatially  this  order  is  short  range  in  nature  we 
interpret  this  as  a  polaron  glass. 


The  data  in  Fig.  1  also  exhibit  dynamics,  in  the  form 
of  quasielastic  scattering  (peaking  at  zero  energy).  Figure 
2  shows  the  temperature  evolution  of  both  the  dynamic 
and  purely  elastic  components,  obtained  on  the  optimally- 
doped  CMR  bilayer  manganite  Laj  24Sri  76Mn207  [3].  On 
decreasing  the  temperature  from  460  K  we  see  that  the 
quasielastic  scattering  narrows  in  energy,  with  a  width  that 
decreases  approximately  linearly  until  T  «  7^  «  3 1 0  K 
while  the  integrated  intensity  of  the  scattering  is  approxi- 
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Fig.  1.  Polaron  dynamics  at  the  peal(  in  tlie  correlation  function  at  260  K, 
just  above    =  252  K,  in  Lag jCag  3Mn03.  The  data  reveal  a  pure  elastic 
signal  (indicating  an  ordered  state)  and  dynamic  quasielastic  scattering. 


mately  temperature  independent.  Below  «  3 10  K  the 
behavior  of  the  quasielastic  width  and  intensity  abruptly 
changes,  and  a  new,  purely  elastic  component,  develops  in 
the  spectrum.  The  elastic  scattering  develops  rapidly 
below  Jg,  and  appears  to  track  the  resistivity  of  the 
material.  This  increase  in  the  elastic  intensity  is  accompa- 
nied by  a  concomitant  decrease  of  the  quasielastic  inten- 
sity, while  the  energy  width  exhibits  a  change  in  slope, 
becoming  approximately  constant  for  Tq<  T <  T^. 

To  further  probe  the  origin  of  this  quasielastic 
scattering  we  compiled  the  map  of  the  dynamical  scatter- 
ing above  7*  shown  in  Fig.  3.  The  inelastic  response 
shows  that  the  dynamical  scattering  still  has  siarprisingly 
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well-defined  correlations  at  these  elevated  temperatures, 
demonstrating  that  the  polarons  are  not  isolated  Jahn-Teller 
Mn^'  deformations,  but  are  dynamically  correlated.  The 
dynamical  correlation  length  obtained  from  these  data  is 
5:  12  A,  comparable  in  size  to  the  static  correlation  length 
found  at  lower  T.  These  dynamic  correlations  are  still 
present  up  to  460  K,  the  highest  T  explored  so  far. 

From  our  measurements  we  can  identify  three 
regimes.  For  7"  >  7*  the  observ  ed  scattering  is  dominated 
by  correlations  that  are  completely  dynamic.  The 
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Fig.  2.  Temperature  dependence  of  the  elastic  (a)  and  quasielastic 
intensity  (b)  of  the  CMR  bilayer  manganite  (Tc  =  114  K).  (c)  Quasielastic 
width  of  the  Lorentzian  response.  Dashed  lines  are  guides  to  the  eye. 
Error  bars  correspond  to  standard  deviations  obtained  from  the  least 
squares  analysis. 


quasielastic  energy  width  from  these  correlations  varies 
linearly  with  T,  which  is  qualitatively  what  is  expected  for  a 
continuous  phase  transition.  The  second  regime  occurs 
for     <  r  <  7*,  and  is  dominated  by  frozen  static 
polaronic  correlations  as  revealed  by  the  rapid  development 
of  the  elastic  component  below  7^.  The  development  of 
the  elastic  component  is  accompanied  by  a  decrease  of  the 
quasielastic  intensity,  while  the  spatial  correlations  remain 
short  range  in  nature.  This  behavior  is  clear  evidence  for  a 
freezing  transition,  analogous  to  the  freezing  transition 
found  in  structural  glasses  or  geometrically  frustrated  spin- 
glasses.  Our  data  strongly  suggest  that  T*  represents  a 
phase  fransition  from  a  polaronic  liquid  to  a  polaronic 
glass.  This  polaron  glass  then  melts  at  Tq  as  the  ferromag- 
netic metallic  state  sets  in. 
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Fig.  3.  Dynamic  polaron  correlations  above  the  glass  transition. 
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Quantum  Impurities 
in  tiie  Two-Dimensional  Spin 
One-Half  Heisenberg  Antiferromagnet 


A phase  transition  is  the  transforma- 
tion of  a  system  from  one  state  to 
another,  e.g.  from  an  ordered  phase  to  a 
disordered  phase.  In  magnetic  iron  for  example, 
the  ordered  ferromagnetic  phase  gives  way  to  the 
disordered  paramagnetic  phase  when  fluctuations 
overwhelm  the  forces  tending  to  order  the  system.  In 
thermal  phase  transitions  the  amplitudes  of  these  fluctua- 
tions increase  with  temperature,  and  the  transition  occurs 
at  a  "critical"  temperature  Tq.  By  contrast,  quantum  phase 
transitions  occur  only  at  T  =  0  K  [1].  The  relevant 
fluctuations  are  thus  quantum  mechanical  in  nature  and  can 
overwhelm  the  tendency  to  order  when  some  other 
parameter  of  the  system  is  varied.  A  conceptual  example 
would  be  a  ferromagnetic  system  at  T  =  0  that  is  continu- 
ously diluted  with  some  non-magnetic  material  until  it 
undergoes  a  transition  to  the  paramagnetic  state. 

The  hope  for  insight  into  the  physics  of  high  tem- 
perature superconductors  has  generated  enormous  interest 
in  the  case  of  the  two-dimensional  (2D)  spin-Vz  square- 
lattice  Heisenberg  antiferromagnet  (5"  =  Vi  SLHAF)  because 
it  models  La2Cu04,  the  parent  compound  of  the  first  high- 
temperature  superconductor.  It  has  been  theorized  that  a 
quantum  phase  transition  might  occur  in  La2Cu04  by 
doping  with  static,  non-magnetic,  atoms  [2].  We  have 
investigated  the  properties  of  the  2D  5"  =  Yi  SLHAF  using 
complementary  neutron  scattering  and  numerical  simula- 
tions. Our  results  demonstrate  that  the  spin-diluted 
Heisenberg  antiferromagnet  La2Cuj.-(Zn,Mg)-04  is  a 
prototype  system  for  square-lattice  site  percolation  in  the 
extreme  quantum  limit  of  5"  =  Vi  [3].  In  the  pure  system  (z 
=  0),  magnetic  Cu2+  ions  form  a  2D  square  lattice  domi- 
nated by  an  antiferromagnetic  nearest-neighbor  exchange 
mediated  by  intervening  oxygen  ions,  which  leads  to  long- 
range  antiferromagnetic  order  at  low  temperature.  Re- 
placement of  the  Cu  spins  with  non-magnetic  Zn2+  and 
Mg2+  ions  randomly  dilutes  the  antiferromagnetic  sheets, 
weakening  the  antiferromagnetic  order.  With  sufficient  site 
dilution  the  system  undergoes  a  geometric  transition  from 
having  one  infinite  connected  cluster  of  Cu  spins  below  the 
percolation  threshold  to  only  finite  disconnected  clusters 
above  the  percolation  threshold  (around  41  %),  as  shown 
in  Fig.  1 .  Our  results  demonstrate  that  random  disorder 


and  quantum  fluctuations  both  play  a  strong  role  in 
determining  the  properties  of  this  system  near  the  percola- 
tion threshold. 

Previous  experiments  [4]  suggested  that 
La2CU].-(Zn,Mg),04  might  become  disordered  at  zero 
temperature  well  below  the  classical  percolation  threshold 
because  of  quantum  fluctuations.  This  question  remained 
unresolved  because  of  a  lack  of  suitable  samples,  but  high- 
quality  single  crystals  of  La2Cuj.^(Zn,Mg)^04  grown  at  the 
T.  H.  Geballe  Laboratory  for  Advanced  Materials  at 
Stanford  University  allowed  us  to  solve  this  quantum 
percolation  problem.  Figure  2  shows  the  dilution  depen- 
dence of  the  antiferromagnetic  ordering  (Neel)  temperature 
T^.  Although  a  purely  2D  system  cannot  have  long-range 
order  at  nonzero  temperature,  small  interlayer  couplings 
allow  long-range  3D  order  to  occur  when  the  2D  correla- 
tions become  large  enough.  Monte  Carlo  calculations  (also 
shown  in  Fig.  2)  indicate  that  this  occurs  at  correlation 
lengths  of  around  100  lattice  constants  for  both  the  pure 
and  diluted  systems.  Our  results  for  T-^{z)  agree  with 
previous  measurements  at  lower  concentrations  [4],  but 
we  find  that  at  higher  concentrations      falls  off  more 
slowly  and  does  not  reach  zero  until  the  percolation 
threshold. 


Fig.  1 .  A  schematic  of  finite  sections  of  the  infinite  lattice  at  dilution 
levels  well  below  (31  %),  just  below  (40.7  %),  and  above  (45  %)  the 
percolation  threshold.  Sites  on  the  infinite  cluster  are  shown  in  red, 
sites  on  finite  disconnected  clusters  are  blue,  and  diluents  are  in  white. 
The  inset  is  a  close-up  view  for  z  =  40.7  %,  showing  the  role  that 
magnetic  Cu  and  non-magnetic  Zn/Mg  ions  play  in  the  experimental 
system. 
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Fig.  2.  Neel  temperature  phase  diagram  of  La2Cu^.2(Zn,IVIg);04.  Single 
crystal  results  are  from  neutron  diffraction,  powder  results  are  from 
SQUID  Magnetometry.  Lines  are  Monte  Carlo  results  for  the  temperature 
at  which  spin-spin  correlations  reach  the  specified  length. 

To  study  the  2D  static  properties  of 
La^Cui..  (Zn,Mg).04,  we  measured  the  2D  instantaneous 
spin-spin  correlation  length  ^  using  a  2-axis  energy- 
integrating  technique  in  the  paramagnetic  phase  above  T-^. 
The  results  are  shown  in  Fig.  3  (colored  symbols)  for  a 
number  of  different  concentrations  and  across  a  broad 
range  of  temperatures.  Quantum  Monte  Carlo  results  for 
the  nearest-neighbor  Heisenberg  model  are  shown  as  black 
symbols.  There  are  no  adjustable  parameters  in  this 
comparison,  and  the  experimental  and  Monte  Carlo  results 
match  exceptionally  well  up  to  the  percolation  threshold. 
We  can  see  that  dilution  not  only  weakens  the  correlations 
at  a  given  temperature,  but  that  it  also  slows  down  the  rate 
at  which  they  grow  as  the  system  is  cooled.  At  low 
concentrations,  ^  scales  exponentially  in  inverse  tempera- 
ture. Near  the  percolation  threshold,  however,  <^  exhibits 
power  law  behavior  in  inverse  temperature,  and  above  the 
percolation  threshold  it  saturates  at  a  finite  length  at  low 
temperature. 

We  have  shown  that  La2Cui.,(Zn,Mg),04  is  an 
excellent  model  system  for  studying  percolation  in  a  2D 
S  =  Vi  Heisenberg  antiferromagnet.  A  long-range-ordered 
ground  state  persists  up  to  the  percolation  threshold  even 
though  quantum  fluctuations  increase  as  dilution  is  in- 
creased. Further  studies  of  the  dynamics  will  provide 
more  insight  into  the  role  of  disorder  in  this  quantum 
many-body  system. 
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Fig.  3.  Spin-spin  correlation  length  (in  units  of  the  lattice  constant)  vs. 
inverse  temperature  (scaled  by  the  nearest  neighbor  exchange  energy  J 
of  the  pure  system)  for  different  concentrations  of  non-magnetic  ions. 
Colored  symbols  are  experimental  results,  obtained  from  measurements 
ofthe  static  structure  factor.  Black  filled  symbols  are  quantum  Monte 
Carlo  results  for  z  =  (8, 20,  31, 35, 41,  and  46)  %  dilution.  Results  for 
z  =  0  %  are  from  [5,  6].  No  adjustable  parameters  are  used  in  the 
comparison. 
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Automobile  Sheet-Metal 
^^Springback":  Residual  Stress 
Measurements  and  Modeling 


A serious  impediment  to  the  use  of 
lighter-weight,  higher-strength 
materials  in  automobile  manufacturing  is  the 
relative  lack  of  understanding  about  how  these 
materials  respond  to  the  complex  forming  opera- 
tions that  go  into  shaping  a  blank  of  metal  into 
automobile  body  parts.  One  of  the  most  vexing  and 
costly  problems  is  "springback"  —  the  tendency  of  sheet 
metal  to  lose  some  of  its  shape  when  it  is  removed  from 
the  die.  Springback  is  very  pronounced  with  two  of  the 
likeliest  candidates  for  weight  reduction:  high-strength  steel 
and  aluminum  alloys,  than  it  is  with  standard  steel.  Unless 
it  is  well  managed  and  taken  into  account  when  the  dies  are 
designed,  it  leads  to  parts  that  are  ill-fitting  and  deviate 
excessively  from  design  intent. 

American  auto  manufacturers,  through  the 
Springback  Project  of  the  USCAR  consortium',  are 
engaged  in  a  major  effort  to  predict  springback  by  means 
of  sophisticated  finite  element  modeling  ("FEM").  How- 
ever, the  accuracy  of  predictions  of  large  strain  plasticity 
under  complex  load  histories,  such  as  those  applied  during 
stamping  processes,  is  uncertain  because  of  incomplete 
validation  of  the  FEM  programs.  Surprisingly,  calculated 
residual  stress,  one  of  the  key  mechanical  properties 
predicted  by  the  state-of-the-art  FEM  codes,  had  not  been 
compared  with  experimental  measurements.  The  present 
work  is  the  first  comprehensive  effort  to  determine  the 
residual  sfresses  of  interest. 

Diffraction  provides  a  powerflil  means  of  very 
accurately  measuring  microstructure,  strains  (from  which 
stresses  are  determined)  and  mechanical  behavior  in  a  way 
not  possible  with  other  techniques.  More  importantly, 
diffraction  facilities  available  to  the  NCNR  include  neufron 
diffraction,  laboratory  x-ray  diffraction  and  synchrotron 
x-ray  diffraction  (at  Argonne's  Advanced  Photon  Source). 
These  constitute  the  full  spectrum  of  diffraction  probes  of 
residual  stress  and  microstructure  for  surface,  sub-surface 
and  bulk  specimens.  The  test  specimens  employed  for  this 
study  were  two  deep-drawn  "Demeri"  cups:  one,  thin- 


'  USCAR  is  the  umbrella  organization  of  DaimlerChrysler,  Ford  and 
General  Motors,  which  was  formed  in  1992  to  further  strengthen  the 
technology  base  of  the  domestic  auto  industry  through  cooperative, 
pre-competitive  research. 
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Fig.  1.  The  deep-drawn  steel  cup.  The  aluminum  cup  was  similar, 
except  for  wall  thickness. 


walled  (a  1  mm)  6022-T4  aluminum;  the  second,  thicker- 
walled  (3.2  mm)  steel.  The  latter  is  shown  in  Fig.  1. 

The  objective  of  this  project  was  to  determine  the 
residual  stresses  in  the  "simple"  model  specimens  formed 
similarly  to  stamped  auto  parts.  Modelers  in  the  USCAR 
consortium  would  use  FEM  to  predict  the  stress  distribu- 
tions to  validate  their  codes.  At  this  time,  the  modeling  part 
of  the  project  is  still  in  progress. 

Two  distinct  experimental  studies  were  performed. 
The  first  utilized  synchrotron  x-rays  to  determine  the 
stress  distribution  in  a  ring  and  pieces  (Fig.  2)  cut  from  the 


Fig.  2.  Pieces  cut  from  the  aluminum  cup  and  examined  by  synchrotron 
x-rays. 
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synchrotron  results 
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Fig.  3.  Axial  and  hoop  stresses  for  the  0.9  mm  thici<  aluminum  cup  as 
determined  by  synchrotron  x-rays. 
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Fig.  4.  Hoop  stresses  as  a  function  of  depth  in  the  3.2  mm  thicl<  steel 
cup,  as  determined  by  neutron  diffraction. 


aluminum  cup  (which  except  for  wall  thickness,  was 
initially  like  the  steel  cup  shown  in  Fig.  1).  The  critical 
point  of  this  study  is  to  determine  the  stress  distribution  in 
the  ring  and,  ultimately,  whether  the  FEM  could  predict  it. 
Directly  related  to  this  was  how  the  measured  stress 
distribution  compared  with  the  simple  linear  depth  depen- 
dence used  in  analytical  calculations  to  predict  the  opening 
of  the  ring  when  cut. 

The  residual  stresses  determined  from  the  x-ray 
diffraction  measurements  are  shown  in  part  in  Fig.  3. 
Representative  neutron  diffraction  results  for  the  steel  cup 
are  shown  in  Fig.  4. 

The  more  complete  results  shown  for  the  aluminum 
ring,  the  first  such  measurements  on  deep-drawn  cups, 
satisfy  both  symmetry'  and  stress  balance  requirements. 
However,  the  stresses  vary  around  the  circumference  and 
in  the  axial  direction,  and  differ  strongly  from  ideal  bending 
stresses.  So  even  for  the  "simple"  model  system,  the 
plastic  deformation  process  and  the  resultant  stresses  are 
very  complex. 


In  summary,  these  results  provide  the  first  through- 
thickness  stress  distributions  by  which  springback  model 
predictions  of  residual  stress  can  be  tested.  Furthermore, 
synchrotron  radiation  and  neutrons  are  the  only  non- 
destructive methods  that  are  able  to  provide  the  necessary 
accuracy  and  spatial  resolution  needed  to  obtain  these 
results.  Finally,  successful  modeling  of  springback  requires 
successful  prediction  of  these  stress  distributions. 
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NCNR  Partnership 
Celebration 

This  year  marks  the  tenth  anniversary  of 
the  start  of  user  experiments  in  the  NCNR  guide 
hall.  The  NCNR  has  maintained  working  partner- 
ships with  a  number  of  government  agencies,  indus- 
trial laboratories,  and  academic  institutions  since  it  was 
commissioned  over  30  years  ago.  Several  key  partner- 
ships are  highlighted  in  this  section.  A  special  event 
celebrating  NCNR  partnerships  and  a  decade  of  success 
was  held  on  August  5,  2002,  with  a  number  of  distin- 
guished visitors  in  attendance.  The  inside  cover  of  this 
report  identifies  persons  addressing  the  audience  at  this 
celebration.  The  speakers  noted  the  key  role  of  partner- 
ships in  the  growth  and  success  of  the  NCNR  and  their 
impact  on  U.S.  neutron  science.  Many  of  their  comments 
referred  to  a  recent  report  from  the  Office  of  Science  and 
Technology  Policy  on  the  "Status  and  Needs  of  Major 
Neutron  Scattering  Facilities  and  Instruments  in  the  United 
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Fig.  1.  Numbers  of  NCNR  Research  Participants  overtime. 


States."  This  report  encourages  the  formation  of  partner- 
ships such  as  those  at  the  NCNR,  especially  between 
federal  agencies,  to  promote  effective  stewardship  of  large 
federal  research  resources.  The  speeches  and  a  guided 
tour  of  the  facility  for  the  visitors  provided  a  satisfying 
perspective  on  the  development  of  the  NCNR  into  one  of 
the  world's  leading  and  most  productive  neutron  scattering 
facilities. 

The  number  of  NCNR  research  participants  has  more 
than  doubled  over  the  last  decade  (Fig.  1)  and  has  contin- 
ued to  increase,  despite  an  extended  shutdown  for  the 
installation  of  the  advanced  cold  source  and  new  cooling 
towers.  As  detailed  in  the  recent  OSTP  report,  the  NCNR 
now  accommodates  more  than  half  of  all  neutron  users  in 
the  US. 

The  NCNR  User  Program 

During  the  past  year,  we  have  made  a  significant  shift 
toward  entirely  Internet-based  proposal  submission  and 
review.  Users  have  been  sending  proposals  to  us  on  our 
Web-page  form  for  more  than  six  years,  but  now  it  has 
become  the  exclusive  method  for  standard  submission  of 
beam  time  proposals.  Local  beam  time  requests  from 
NIST  staff  are  managed  through  the  same  system.  Peer 
reviews  of  proposals  are  also  communicated  to  us  elec- 
tronically through  a  Web  form,  resulting  in  faster  and  more 
efficient  decisions  on  instrument  time,  to  the  benefit  of 
users  and  their  research  programs. 

Two  calls  for  proposals  were  made  in  the  past  year, 
and  more  than  390  proposals  were  received  in  response. 
The  proposed  experiments  show  increasing  activity  in 
research  on  soft  condensed  matter  and  bioscience,  al- 
though more  mature  areas  such  as  magnetism  and  polymer 
research  have  held  their  own.  Proposals  are  now  nearly 
equally  divided  between  those  for  structural  investigations, 
mainly  SANS  and  reflectometry,  and  studies  of  dynamics, 
requesting  one  of  the  newer  high-resolution  specfrometers 
in  the  guide  hall. 

After  a  thorough  review  process  by  external  referees 
and  by  the  NCNR  Program  Advisory  Committee  (PAC), 
approved  proposals  were  allocated  beam  time.  The  PAC  is 
a  panel  of  accomplished  scientists  with  expertise  across  a 
broad  range  of  neutron  methods  and  scientific  disciplines. 
It  is  the  body  primarily  responsible  for  proposal  review  and 
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recommending  user  policies  for  the  NCNR,  working 
closely  with  the  Center's  Director  and  staff.  Its  current 
membership  includes  Sanat  Kumar  (Rensselaer  Polytechnic 
histitute,  chair),  Robert  M.  Briber  (University  of  Mary- 
land), Michael  K.  Crawford  (DuPont),  Yumi  Ijiri  (Oberlin 
College),  Dieter  K.  Schneider  (Brookhaven  National 
Laboratory),  Kenneth  Herwig  (Oak  Ridge  National  Labora- 
tory), Michael  Kent  (Sandia  National  Laboratory),  Robert 
Leheny  (Johns  Hopkins  University),  John  Tranquada 
(Brookhaven  National  Laboratory),  and  Sossina  Haile 
(California  Institute  of  Technology). 

The  Center  for  High  Resolution 
Neutron  Scattering 

One  of  the  most  important  partnerships  between  the 
NCNR  and  another  agency  is  with  the  National  Science 
Foundation  (NSF)  through  its  support  of  the  Center  for 
High  Resolution  Neutron  Scattering  (CHRNS).  Five 
instruments  in  the  NCNR  guide  hall,  and  one  in  the  thermal 
instrument  area,  are  operated  by  CHRNS,  comprising  a 
crucial  component  of  the  user  program.  The  instruments 
include  a  30  m  SANS  instrument,  the  SPINS  triple-axis 
spectrometer,  the  neutron  spin-echo  spectrometer,  the 
high-flux  backscattering  spectrometer,  the  disk-chopper 
time-of-flight  spectrometer,  and  the  ultra-small  angle 
scattering  diffractometer.  Most  of  the  instrument  time 
allocated  by  the  PAC  now  goes  to  experiments  carried  out 
on  CHRNS  instruments.  An  additional  SANS 
diffractometer  on  neutron  guide  NG-1,  which  up  to  now 
has  been  used  primarily  for  NIST  programmatic  research, 
is  being  upgraded  to  a  more  powerful  10  m  instrument 
with  a  new  detector,  and  will  be  made  available  to  CHRNS 
users.  In  aggregate,  the  CHRNS  instruments  provide 
structural  information  on  a  length  scale  from  1  nm  to 
~  1 0  fim,  and  dynamical  information  on  energy  scales 
from  «  30  neV  to  «  100  meV.  These  are  the  widest  ranges 
accessible  at  any  neutron  research  center  in  North 
America. 

Collaborations 

Direct  collaborations  on  specific  experiments  remain 
a  common  way  for  users  to  pursue  their  ideas  using 
NCNR  facilities,  accounting  for  approximately  half  of  the 
number  of  instrument-days.  The  thermal-neutron  triple- 
axis  spectrometers  are  mainly  scheduled  in  this  way.  Most 
of  the  time  reserved  for  NIST  on  these  and  all  other  NCNR 
instruments  is  devoted  to  experiments  that  are  collabora- 
tions with  non-NIST  users.  Collaborative  research 
involving  external  users  and  NIST  scientists  often  pro- 
duces results  that  could  be  not  obtained  otherwise. 


Fig.  2.  Sanat  Kumar,  Dieter  Schneider,  and  Rob  Briber  of  the  PAC 
consider  proposals  for  beam  time. 

Another  mode  of  access  to  the  NCNR  is  through 
more  formal  research  consortia.  Groups  of  researchers 
from  various  institutions  join  forces  to  build  and  operate  an 
instrument.  Typically,  a  substantial  fraction  of  the  time  on 
the  instrument  is  then  reserved  for  the  consortium  mem- 
bers, and  the  remaining  time  is  allocated  to  general  user 
proposals.  For  example,  a  group  including  ExxonMobil, 
the  University  of  Minnesota,  and  NIST  cooperates  on  the 
NG-7  30  m  SANS  instrument.  Similar  arrangements 
involving  other  consortia  apply  for  the  horizontal-sample 
reflectometer,  the  high-resolution  powder  diffractometer, 
the  filter-analyzer  spectrometer,  and  the  neutron  spin-echo 
spectrometer. 

Cold  Neutrons  for  Biology  and 
Technology 

An  important  new  partnership  with  five  universities, 
NIST  and  the  National  Institutes  of  Health  is  now  proceed- 
ing with  the  design  and  construction  of  a  new  reflectome- 
ter dedicated  to  bioscience  and  biotechnology.  The 
University  of  California  at  Irvine  leads  the  consortium 
whose  purpose  is  to  provide  a  new  national  resource  for 
investigations  in  structural  biology.  Research  leaders  from 
Penn  State,  Rice,  Carnegie  Mellon,  Duke,  and  the  Los 
Alamos  National  Laboratory  are  included  in  the  consor- 
tium. 

independent  Programs 

There  are  a  number  of  long  standing 
programs  located  at  the  NCNR  that  involve 
other  parts  of  NIST,  universities, 
industrial  laboratories,  or  other 
government  agencies. 


The  Nuclear 
Methods  Group 

(Analytical  Chemistry 
Division,  Chemical  Science  and 
Technology  Laboratory)  has  as  its 
principal  goals  the  development  and 
application  of  nuclear  analytical  techniques 
for  the  determination  of  elemental  composi- 
tions with  greater  accuracy,  higher  sensitivity 
and  better  selectivity.  A  high  level  of  competence 
has  been  developed  in  both  instrumental  and  radio- 
chemical neutron  activation  analysis  (INAA  and 
RNAA).  In  addition,  the  group  has  pioneered  the  use  of 
cold  neutron  beams  as  analytical  probes  with  both  prompt 
gamma  activation  analysis  (PGAA)  and  neutron  depth 
profiling  (NDP).  PGAA  measures  the  total  amount  of  a 
particular  analyte  present  throughout  a  sample  by  the 
analysis  of  the  prompt  gamma-rays  emitted  upon  neutron 
capture.  NDP,  on  the  other  hand,  determines  concentra- 
tions of  several  important  elements  (isotopes)  as  a  fianction 
of  depth  within  the  first  few  micrometers  of  a  surface  by 
energy  analysis  of  the  prompt  charged-particles  emitted 
during  neutron  bombardment.  These  techniques  (INAA, 
RNAA,  PGAA,  and  NDP)  provide  a  powerful  combination 
of  complementary  tools  to  address  a  wide  variety  of 
analytical  problems  of  great  importance  in  science  and 
technology,  and  are  used  to  help  certify  a  large  number  of 
NIST  Standard  Reference  Materials. 

A  large  part  of  the  group's  efforts  is  directed  towards 
the  exploitation  of  the  analytical  applications  of  the  guided 
cold-neutron  beams  available  at  the  NCNR.  The  Group's 
involvement  has  been  to  design  and  construct  state-of-the- 
art  cold  neutron  instruments  for  both  PGAA  and  NDP  and 
provide  facilities  and  measurements  for  outside  users, 
while  retaining  and  utilizing  our  existing  expertise  in  INAA 
and  RNAA. 

The  Neutron  Interactions  and  Dosimetry  Group 
(Physics  Laboratory)  provides  measurement  services, 
standards,  and  fundamental  research  in  support  of  NIST's 
mission  as  it  relates  to  neutron  technology  and  neutron 
physics.  The  national  and  industrial  interests  served 
include  scientific  instrument  calibration,  electric  power 
production,  radiation  protection,  defense  nuclear  energy 
systems,  radiation  therapy,  neutron  radiography,  and 
magnetic  resonance  imaging. 

The  Group's  research  may  be  represented  as  three 
major  activities.  The  first  is  Fundamental  Neutron  Physics 
including  magnetic  trapping  of  ultracold  neutrons,  opera- 
tion of  a  neutron  interferometry  and  optics  facility,  devel- 


opment of  neutron  spin  filters  based  on  laser  polarization  of 
3He,  measurement  of  the  beta  decay  lifetime  of  the  neu- 
tron, and  investigations  of  other  coupling  constants  and 
symmetries  of  the  weak  interaction.  This  project  involves 
a  large  number  of  collaborators  Irom  universities  and 
national  laboratories. 

The  second  is  Standard  Neutron  Fields  and  Applica- 
tions utilizing  both  thermal  and  fast  neutron  fields  for 
materials  dosimetry  in  nuclear  reactor  applications  and  for 
personnel  dosimetry  in  radiation  protection.  These  neutron 
fields  include  thermal  neutron  beams,  "white"  and  mono- 
chromatic cold  neutron  beams,  a  thermal-neutron-induced 
fission  neutron  field,  and  -^-Cf  fission  neutron  fields, 
both  moderated  and  unmoderated. 

The  third  is  Neutron  Cross  Section  Standards 
including  experimental  advancement  of  the  accuracy  of 
neutron  cross  section  standards,  as  well  as  evaluation, 
compilation  and  dissemination  of  these  standards. 

The  Polymers  Division  of  the  Materials  Science  and 
Engineering  Laboratory  has  several  program  elements  at 
the  NCNR  covering  the  technological  areas  of  electronics, 
multiphase  polymeric,  and  biological  materials.  The 
purpose  of  the  electronics  materials  program  is  to  help  the 
U.S.  microelectronics  industry  address  their  most  pressing 
materials  measurement  and  standards  issues.  The  struc- 
ture of  nanoporous  thin  films,  used  in  integrated  circuits  as 
low-dielectric-constant  insulators,  is  characterized  through 
a  combination  of  small  angle  neutron  scattering  (SANS), 
x-ray  reflectivity,  and  ion  scattering  measurements. 
Advances  such  as  vapor  contrast  matching  techniques 
provide  unique  information  about  critical  parameters  such 
as  the  porosity,  average  pore  size,  and  wall  mass  density. 
In  addition,  several  instruments  at  NCNR  provide  unique 
high-resolution  measurements  that  guide  the  development 
and  fundamental  understanding  of  next-generation  poly- 
meric photoresist  materials  for  sub- 1 00  nm  lithography. 
Measurements  include  the  segmental  dynamics  of  the 
polymers  that  govern  the  mobility  of  an  acid  catalyst, 
direct  measurement  of  the  deprotection  reaction  that  leads 
to  patterned  areas,  and  the  structural  characterization  of 
lithographically  prepared  structures. 

In  the  multiphase  polymeric  materials  program,  the 
objective  is  to  understand  underlying  principles  of  phase 
behavior,  local  structure,  and  flow  that  control  morphology 
and  structure  during  processing  of  multiphase  materials 
such  as  polymer  blends  and  nanocomposites  (polymers 
filled  with  clay  platelets,  carbon  nanotubes,  or  triblock 
copolymer  micelles).  SANS  and  reflectivity  measurements 
in  equilibrium,  in  transient  conditions,  and  under  external 
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fields,  provide  essential  information  for  general  understand- 
ing as  well  as  for  specific  application  of  polymer  blend/ 
alloy  systems.  For  example,  the  structure  induced  by 
shear  alignment  of  nanoparticle/polymer  solutions  has  been 
elucidated  with  SANS,  small  angle  light  scattering,  and 
rheological  measurements.  Customers  include  material 
producers  and  users,  ranging  fi^om  chemical,  rubber,  tire, 
and  automotive  companies,  to  small  molding  and  com- 
pounding companies.  The  research  focus  on  polymeric 
materials  includes  commodity,  engineering  and  specialty 
plastic  resins,  elastomers,  coatings,  adhesives,  films, 
foams,  and  fibers. 

The  biological  materials  initiative  is  committed  to 
demonstrating  the  use  of  inelastic  neutron  scattering  (INS) 
as  a  unique  measurement  tool  for  understanding  the 
dynamics  in  several  important  biological  systems.  For 
example,  INS  has  been  used  to  elucidate  changes  in  protein 
dynamics  that  occur  in  the  final  stages  of  protein  folding. 
Likewise,  INS  and  low  fi-equency  Raman  scattering  have 
been  developed  as  measurement  tools  to  identify  the 
dynamical  aspects  of  a  good  sugar  or  polyalcohol  for 
protein  stabilization  and  cryopreservation.  Effective 
preservatives  must  suppress  the  high  frequency  motions 
that  are  precursors  to  protein  unfolding  and  denaturation. 

The  ExxonMobil  Research  and  Engineering 
Company  is  a  member  of  the  Participating  Research  Team 
(PRT)  that  operates,  maintains,  and  conducts  research  at 
the  NG-7  30  m  SANS  instrument  and  the  NG-5  Neutron 
Spin  Echo  Spectrometer.  Their  mission  is  to  use  those 
instruments,  as  well  as  other  neutron  scattering  techniques 
available  to  them  at  NCNR,  in  activities  that  complement 
research  at  ExxonMobil's  main  laboratories  as  well  as  at  its 
affiliates'  laboratories  around  the  world.  The  aim  of  these 
activities  is  to  deepen  understanding  of  the  nature  of 
ExxonMobil's  products  and  processes,  so  as  to  improve 
customer  service  and  to  improve  the  return  on  sharehold- 
ers' investment.  Accordingly,  and  taking  full  advantage  of 
the  unique  properties  of  neutrons,  most  of  the  experiments 
use  SANS  or  other  neutron  techniques  to  study  the 
structure  and  dynamics  of  hydrocarbon  materials,  espe- 
cially in  the  fields  of  polymers,  complex  fluids,  and 
petroleum  mixtures.  ExxonMobil  regards  its  participation 
in  the  NCNR  and  collaborations  with  NIST  and  other  PRT 
members  not  only  as  an  excellent  investment  for  the 
company,  but  also  as  a  good  way  to  contribute  to  the 
scientific  health  of  the  nation. 

The  U.S.  Army  Research  Laboratory  is  the  current 
manager  of  a  30-year  NCNR-Army  partnership  sponsoring 
a  research  team.  This  group  performs  materials  research 


and  engineering  in  support  of  Army  needs  and  missions 
and  works  jointly  with  NIST/NCNR  in  developing  and 
improving  neutron  methods  and  instrumentation.  A  key 
contribution  of  Army  researchers  at  the  NCNR  is  a  highly 
successful  neutron  facility  to  map  residual  stress  and 
texture  within  large  components  such  as  jet  turbine  blades 
and  structural  welds.  This  resource  is  critical  to  develop- 
ment and  performance  of  industrial  and  transportation 
systems  and  military  hardware. 

The  Center  for  Food  Safety  and  Applied  Nutri- 
tion, U.S.  Food  and  Drug  Administration  (FDA),  directs 
and  maintains  a  neutron  activation  analysis  (NAA)  facility 
at  the  NCNR.  This  facility  provides  agency-wide  analytical 
support  for  special  investigations  and  applications  research, 
complementing  other  analytical  techniques  used  at  FDA 
with  instrumental,  neutron-capture  prompt-gamma,  and 
radiochemical  NAA  procedures,  radioisotope  x-ray  fluores- 
cence spectrometry  (RXRFS),  and  low-level  gamma-ray 
detection.  This  combination  of  analytical  techniques 
enables  diverse  multi-element  and  radiological  information 
to  be  obtained  for  foods  and  related  materials.  The  NAA 
facility  supports  agency  quality  assurance  programs  by 
developing  in-house  reference  materials,  by  characterizing 
food-related  reference  materials  with  NIST  and  other 
agencies,  and  by  verifying  analyses  for  FDA's  Total  Diet 
Study  Program.  Other  studies  include  the  development  of 
RXRFS  methods  for  screening  foodware  for  the  presence 
of  Pb,  Cd  and  other  potentially  toxic  elements,  use  of 
instrumental  NAA  to  investigate  bromate  residues  in  bread 
products,  and  use  of  prompt-gamma  NAA  to  investigate 
boron  nutrition  and  its  relation  to  bone  strength. 

The  Smithsonian  Center  for  Materials  Research 
and  Education  Nuclear  Laboratory  for  Archeological 
Research  (SCMRE)  has  chemically  analyzed  over  24  000 
artifacts  by  INAA  at  the  NCNR  over  the  last  25  years. 
SCMRE 's  research  programs  draw  extensively  upon  the 
collections  of  the  Smithsonian,  as  well  as  those  of  national 
and  international  institutions.  The  chemical  analyses 
provide  a  means  of  linking  these  diverse  collections 
together  to  study  continuity  and  change  involved  in  the 
production  of  ceramic  objects.  INAA  data  are  used  to 
determine  if  groups  of  ceramics  have  been  made  from 
the  same  or  different  raw  materials.  The  ceramics 
then  can  be  attributed  to  geographic  regions, 
specific  sources,  workshops  and  even  indi- 
vidual artists.  The  ability  to  combine 
chemical  data  for  semi-complete  or 
whole  vessels  with  that  derived  from 
the  analysis  of  materials  recov- 


ered  from  recent 
excavation  programs 
provides  many  new  opportu- 
nities for  study.  Museum  vessels 
that  are  found  to  be  similar  and 
attributable  to  some  geographic  location 
provide  a  basis  for  exploring  changing 
aspects  of  style,  iconography,  textual  history 
and  even  ideology.  In  an  integrative  manner,  the 
ESfAA  program  enhances  the  importance  of 
collection  materials  for  the  study  of  the  past. 
A  number  of  universities  have  also  established 
long-term  programs  at  the  NCNR.  The  University  of 
Maryland  is  heavily  involved  in  the  use  of  the  NCNR, 
and  maintains  several  researchers  at  the  facility.  Jolins 
Hopkins  University  participates  in  research  programs  in 
magnetism,  soft  condensed-matter  physics,  bioscience, 
and  in  instrument  development  at  the  NCNR.  The  Univer- 
sity of  Pennsylvania  is  working  to  help  develop  biological 
applications  of  neutron  scattering.  It  is  also  participating  in 
the  second  stage  construction  of  the  filter  analyzer  neutron 
spectrometer,  along  with  the  University  of  California  at 
Santa  Barbara,  DuPont,  Hughes,  and  Allied  Signal. 
The  University  of  Minnesota  participates  in  two  PRTs,  the 
NG-7  30  m  SANS  and  the  NG-7  reflectometer.  The  Univer- 
sity of  Massachusetts  also  participates  in  the  latter  PRT. 


Eighth  Annual  Summer  School 

With  support  from  the  National  Science  Foundation, 
NCNR  and  CHRNS  held  their  annual  Summer  School  on 
Neutron  Scattering  from  June  3-7,  2002.  The  subject 
matter  of  the  school  alternates  in  successive  years  between 
techniques  designed  to  investigate  structure,  and  those 
devoted  to  spectroscopy  and  dynamics.  The  course  this 
year  focused  on  structural  studies  with  SANS  and  neutron 
reflectometry  (NR).  Forty  graduate  students  and 
postdoctoral  fellows  from  institutions  across  the  US 
attended  the  Summer  School  this  year.  The  format  of  the 
meeting  in  recent  years  has  emphasized  hands-on  experi- 
ence with  CHRNS  instruments.  For  much  of  the  week, 
students  in  small  subgroups  circulated  among  various 
instruments  set  up  to  perform  illustrative  experiments,  and 
analyzed  data  using  NCNR-supplied  resources.  Among  the 
experiments  was  a  characterization  by  SANS  of  the 
structure  of  clay  platelets  dispersed  in  water.  The  neutron 
reflectometry  experiments  comprised  measurements  on 
thin  films  and  multilayer  materials.  On  the  final  day,  the 
subgroups  made  presentations  on  the  results  of  their 
studies.  According  to  student  evaluations,  the  school  was 
highly  effective  in  introducing  them  to  the  principles  and 
practice  of  structural  studies  with  neufron  instruments. 


Fig.  3.  Eighth  Annual  Summer  School  on  Neutron  Scattering  participants  concentrate  on  neutron 
reflectometry  results  under  the  guiding  hand  of  NCNR's  Sushil  Satija. 
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T'he  NIST  neutron  source  (NBSR)  operated  for  153 
full  power  (20  MW)  days  or  approximately  94  %  of 
the  scheduled  time  for  this  year.  A  typical  operating  year 
consists  of  seven  cycles.  A  cycle  has  38  days  of  continu- 
ous ftill  power  operation,  followed  by  1 1  days  of  shutdown 
maintenance,  refiieling,  and  startup  preparations.  This 
year,  the  NBSR  was  in  a  scheduled  extended  shutdown  for 
maintenance  and  facility  improvements  until  March  6, 
2002.  Those  improvements  included  a  new  cold  source, 

new  transformers  and  

switchgear  for  the 
building  electrical 
system,  a  new  plume- 
abatement  cooling 
tower,  and  new  shim 
arm  seal  assemblies. 
The  new  seal  assem- 
blies have  characteris- 
tics superior  to  the 
old  assemblies  and 
provide  greater 
assurance  of  depend- 
able shim  arm 
operation.  The 
replacement  of  the 
electrical  equipment 

and  the  cooling  tower  will  increase  the  reliability  of  the 
NBSR  over  the  next  25  years.  In  pursuit  of  that  same 
goal,  detailed  inspections  of  vital  systems  and  components 
were  completed  using  ultrasonic  probes,  visual  examina- 
tions, and  remote  examinations  with  miniature  cameras  and 
fiber-optic  imaging  equipment.  The  only  major  mainte- 
nance item  scheduled  for  completion  this  shutdown  was  a 


Seymour  H.  Weiss,  Chief, 
Operations  &  Engineering 


determination  of  the  status  of  the  Thermal  Column  Tank 
Cooling  System.  After  a  thorough  investigation,  the 
capability  of  the  Thermal  Column  was  preserved  by 
Reconfiguring  the  cooling  system  and  replacing  original 
components  within  the  biological  shield.  A  new  tank  has 
been  ordered  and  will  be  available  if  it  becomes  necessary 
to  replace  the  tank  in  the  future. 

A  number  of  important  personnel  issues  were 
addressed  in  FY2002,  including  the  effect  of  a  surge  of 

impending  retire- 
ments on  the  NBSR. 
All  licensed  opera- 
tors were  requalified 
this  year.  Three 
new  senior  operators 
were  licensed  by  the 
Nuclear  Regulatory 
Commission,  a  new 
deputy  chief  for 
Operations  was 
appointed,  and 
significant  progress 
was  made  on 
addressing  recruit- 
ment and  retention 
questions. 

Substantial  changes  to  facility  security  were  com- 
pleted this  year  to  fulfill  the  intent  of  the  President's 
directives  on  homeland  security.  A  vehicle  exclusion  area 
was  established  and  surveillance  equipment  was  improved. 
Other  changes  are  on-going  and  Operations  &  Engineering 
is  consulting  with  the  U.S.  Nuclear  Regulatory  Commis- 
sion on  security  matters  on  a  regular  basis. 


Tliomas  J.  IVlyers,  Deputy  Chief, 
and  Chief  of  Operations 
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DAVE  and  Other  Software 
Developments 

Scientists  at  the  NCNR,  with  support 
from  the  National  Science  Foundation  via  CHRNS, 
have  developed  a  new  software  tool  for  the  reduc- 
tion, visualization,  and  analysis  of  neutron  inelastic 
scattering  data.  DAVE,  short  for  the  Data  Analysis  and 
Visualization  Environment,  is  an  integrated  suite  of 
interactive  software  tools  with  a  visual  interface  for 
treating  and  analyzing  neutron  inelastic  scattering  data  sets. 
Using  a  powerful  graphical  interface,  users  can  reduce 
their  data  from  one  of  the  inelastic  spectrometers,  make 
cuts  through  the  scattering  ftinction,  and  fit  them  with  the 
lineshape  of  their  choice  from  a  library  of  model  ftinctions. 

The  goal  of  the  DAVE  software  package  is  to  allow 
users,  whose  neutron  scattering  experience  ranges  from 
casual  to  expert,  to  reduce,  visualize,  and  analyze  their 
inelastic  neutron  scattering  data  sets  with  a  minimum  of 
effort  by  using  a  set  of  visually  intuitive  tools.  The 
software  has  been  in  use  on  many  of  the  inelastic  neutron 
scattering  instruments  since  January  2002  and  was 
released  to  the  user  community  in  July.  The  source  code  is 
available  from  the  NCNR  website  to  run  on  PC,  Linux,  and 
Mac  platforms.  Project  and  download  information  are 
available  at  http://vvww.ncnr.nist.gov/dave. 

New,  easy  to  use  software  for  reducing  and  fitting 
neutron  reflectivity  data  is  under  active  development  at  the 
NCNR.  The  new  data  reduction  tools  employ  an  intuitive 


S  (Q,  omega)  Intensity 


0.2    0.4    0.6    0.8    1.0  1.2 
mean  0(A-'') 


Fig.  1.  An  example  of  DAVE  output:  Dispersion  curve  of  a  magnetic 
excitation  as  measured  using  the  Disl<  Chopper  Spectrometer. 


graphical  interface  developed  using  open  source  software 
to  read  and  process  a  variety  of  raw  data  files,  including 
specular  scans,  background  measurements,  rocking 
curves,  polarized  beam,  and  slit  scans.  Using  these  new 
tools,  the  user  can  perform  a  variety  of  data  correction  and 
manipulation  tasks  to  produce  reduced  data  sets  of  the 
sample  reflectivity  as  a  ftinction  of  wavevector.  The 
software  supports  a  wide  range  of  data  formats  from  both 
the  NG-7  and  NG- 1  reflectometers  and  single  channel  or 
multichannel,  position-sensitive  detectors. 
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Fig.  2.  Graphical  interface  image  of  neutron  reflectivity  data  obtained 
using  a  position  sensitive  detector  with  corresponding  slices  of  constant 
wavevector. 

A  new  graphical  interface  has  also  been  developed 
for  the  reflectivity- fitting  program  'mlayer'  [1].  This 
program  allows  users  to  model  and  fit  scattering  length 
density  profiles  to  measured  reflectivity  data.  The  new 
interface  for  this  program  provides  a  wide  set  of  graphical 
tools  for  adjusting  and  refining  model  parameters,  control- 
ling the  fit,  and  constraining  the  model  profile.  A  similar 
graphical  interface  has  also  been  developed  for  the  pro- 
gram that  treats  polarized  neutron  reflectivity  data.  Further 
information  on  these  tools  is  available  at  the  website: 
www.ncnr.nist.gov/programs/reflectometry/software. 

New  Monochromator  Drum  Shields 

This  year  the  NCNR  ordered  two  new  monochroma- 
tor drum  shields  for  the  new  BT-7  neutron  triple  axis 
spectrometer  and  for  the  BT-4  filter  analyzer  instrument. 
Fabrication  and  testing  of  one  of  these  drum  shields  has 
been  completed  and  delivery  of  both  is  planned  for  Dec. 
2002. 

Monochromator  drum  shields  are  a  key  component  in 
crystal  spectrometers,  forming  a  thick,  moving  shield  that 
surrounds  the  intensely  illuminated  monochromator  and 
allowing  the  diffracted  beam  to  be  directed  at  a  variety  of 


54 


Facility  Developments 


take-off  angles  towards  the  sample.  The  design  of  the 
new  drum  shields  is  based  on  a  concept  originally  devel- 
oped for  the  Advanced  Neutron  Source  project  at  Oak 
Ridge  National  Laboratory  which  uses  pairs  of  counter- 
opening  shield  wedges  to  allow  the  large  diameter  neutron 
beams  to  enter  the  drum.  For  the  NCNR  drums,  this 
design  concept  was  adapted  for  the  specific  dimensions  of 
the  NCNR  source  and  to  optimize  the  manufacturability  of 
the  shield.  This  double  wedge  design  uses  a  precision  cam 
track  to  position  the  shield  wedges  around  the  20  cm  beam 
opening  as  the  drum  is  rotated  to  various  take-off  angles. 
The  drum  shield  is  rotated  using  a  stepper  motor  driven 
gear  system.  The  take-off  angle  range  is  between  0°  and 
120°  and  the  drum  can  be  built  into  either  clockwise  or 
counter-clockwise  rotating  models.  The  NCNR  will  be  the 
first  neutron  laboratory  to  implement  this  double-wedge 
design  for  a  monochromator  shield. 

Sample  Environment  Equipment 

The  suite  of  available  sample  support  equipment  at 
the  NCNR  continues  to  be  upgraded  and  enhanced.  This 
year,  the  superconducting  magnet  system  with  dilution 
refrigeration  was  upgraded  from  the  seven  Tesla  coils 
available  last  year  to  the  full  design  field  coils  that  provide 
11.5  Tesla.  This  system  has  demonstrated  a  22  mK  base 
temperature  and  full  magnetic  field  strength. 

The  NCNR  also  received  a  new  nine  Tesla  horizontal 
field  superconducting  magnet  with  single  crystal  windows 
using  a  design  specifically  optimized  for  small  angle 
scattering  experiments.  This  system  uses  a  split  coil 
design  and  windows  in  either  direction  to  allow  the 
magnefic  field  to  be  directed  parallel  or  perpendicular  to  the 
neutron  beam  axis.  Extensive  heat  shielding  extends 
operation  of  the  magnet  to  the  higher  sample  temperatures 
needed  for  biology  or  polymer  physics  experiments 
without  significantly  altering  liquid  helium  consumption. 

Other  new  equipment  made  available  this  year 
includes  a  pulse  tube  closed  cycle  refrigerator  (CCR)  with 
an  operating  temperature  range  of  3  K  to  325  K,  two  low 
temperature  CCRs  (10  K  to  325  K),  one  high  temperature 
CCR  (15  K  to  800  K),  one  50  mm  helium  cryostat  (1.5  K 
to  300  K)  and  fifteen  additional  temperature  controllers. 

Facility  Improvements 

Taking  advantage  of  the  opportunity  provided  by  an 
extended  interruption  in  routine  operations,  the  NCNR  staff 
completed  several  improvements  to  the  facility  during  the 
outage  early  this  year.  The  helium  gas  blanketing  system, 
which  provides  a  continuous  flow  of  helium  gas  to  the  first 


Fig.  3.  George  Baltic  (center,  NCNR)  confirms  completion  of  a 
monochromator  drum  with  Paul  Brand  (left,  NCNR)  and  Rich  Porter, 
president  of  Ingersol  Contract  Manufacturing  Company. 

1.5  m  of  neutron  guide  closest  to  the  cold  source,  was 
replaced  with  a  new  gas  control  system.  Sensors  and 
controllers  on  the  guide  vacuum  systems  were  replaced 
and  upgraded  during  the  cold  source  installation.  A 
computerized  supervisory  and  control  system  for  the 
facility  was  expanded  to  include  monitoring  of  the  guide 
vacuum  pumping  system. 

Several  improvements  were  made  in  CI 00,  the 
location  of  the  thermal  neutron  instruments.  These  include 
the  installation  of  wiring  and  breaker  panels  that  are  part  of 
a  new  power  distribution  system  for  the  thermal  instru- 
ments. The  new  system  provides  each  neutron  station 
with  independent  panels  and  outlet  locations  at  the  floor 
level  adjacent  to  the  experimental  work  area.  A  new 
catwalk  system  was  also  designed  for  accessing  the  face 
of  the  biological  shield  above  the  thermal  instruments. 
This  will  be  a  significant  safety  and  performance  improve- 
ment for  maintenance  activities  in  this  area.  A  prototype 
for  a  perimeter  shield  was  also  accepted  this  year.  The 
perimeter  shields  surround  the  thermal  neutron  spectrom- 
eters and  separate  the  operating  area  of  the  spectrometers 
from  other  work  areas  on  the  floor. 

In  the  Guide  Hall,  the  concrete  floor  was  repainted. 
The  original  surface  was  badly  worn  and  needed  replace- 
ment. A  three-color  design  was  adopted  to  help 
differentiate  work  areas,  walkways,  and  fire 
protection  equipment. 
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Instruments  and  Contacts 


High  resolution  powder  diffractometer  (BT-1) 

B.  H.  Toby,  (301)  975-4297,  brian.toby@nist.gov 
C.  Jones,  (301)  975-4507,  camille@nist.gov 
Q.  Huang,  (301)  975-6164,  qing.huang@nist.gov 
J.  K.  Stalick,  (301)  975-6223,  judy.stalick@nist.gov 

DARTS,  Residual  stress  and  texture  diffractometer  (BT-8) 

H.  J.  Prask,  (301)  975-6226,  liank@nist.gov 
T.  Gnaeupel-Herold,  (301)  975-5380,  thomas.gnaeupel- 
herold@nist.gov 

30-m  SANS  instrument  (NG-7) 

C.  J.  Glinka,  (301)  975-6242,  cglinka@nist.gov 

B.  S.  Greenwald,  (301)  975-5795,  bsgreen@nist.gov 
L.  Porcar,  (301)  975-5049,  iionel.porcar@nist.gov 

30-m  SANS  instrument  (NG-3)  (NIST/NSF-CHRNS) 

B.  Hammouda,  (301)  975-3961,  hamnnouda@nist.gov 
S.  R.  Kline,  (301)  975-6243,  steven.kline@nist.gov 

D.  Ho,  (301)  975-6422,  derek.ho@nist.gov 
T.  K.  Misra,  (301)  975-6650,  tmisra@nist.gov 

8-m  SANS  instrument  (NG-1) 

D.  Ho,  (301)  975-6422,  derek.ho@nist.gov 

C.  J.  Glinka,  (301)  975-6242,  cglinka@nist.gov 

J.  G.  Barker,  (301)  975-6732,  john.barker@nist.gov 

USANS,  Perfect  crystal  SANS  (BT-5)  (NIST/NSF-CHRNS) 

J.  G.  Barker,  (301)  975-6732,  john.barker@nist.gov 
IVI.  H.  Kim,  (301)  975-6469,  kimmh@nist.gov 
C.  J.  Glinka,  (301)  975-6242,  cglinka@nist.gov 

Cold  neutron  reflectometer-vertical  sample-polarized  beam 
option  (NG-1) 

C.  F.  iVIajkrzak,  (301)  975-5251,  cmajkrzak@nist.gov 
J.  A.  Dura,  (301)  975-6251,  jdura@nist.gov 

Cold  neutron  reflectometer-horizontal  sample  (NG-7) 

S.  K.  Satija,  (301)  975-5250,  satija@nist.gov 

Y.-S.  Seo,  (301)  973-5391,  young-soo.seo@nist.gov 

Triple-axis  polarized-beam  spectrometer  (BT-2) 

J.  W.  Lynn,  (301)  975-6246,  jeff.lynn@nist.gov 

Triple-axis  fixed  incident  energy  spectrometer  (BT-7) 

J.  W.  Lynn,  (301)  975-6246,  jeff.lynn@nist.gov 

Triple-axis  spectrometer  (BT-9) 

R.  W.  Emn,  (301)  975-6245,  renwin@nist.gov 
P.  M.  Gehring,  (301)  975-3946,  pgehring@nist.gov 

SPINS,  Spin-polarized  triple-axis  spectrometer  (NG-5) 
(NIST/NSF-CHRNS) 

S.  Park,  (301)  975-8369,  sungil.park@nislgov 
S.-H.  Lee,  (301)  975-4257,  seung-hun.lee@nist.gov 
P.  M.  Gehring,  (301)  975-3946,  pgehring@nist.gov 


FANS,  Filter-analyzer  neutron  spectrometer  (BT-4) 

T.  J.  Udovic,  (301)  975-6241,  udovic@nist.gov 

C.  IVI.  Brown,  (301)  975-5134,  craig.brown@nist.gov 
J.  Leao,  (301)  975-8867,  jleao@nist.gov 

D.  A.  Neumann,  (301)  975-5252,  dan@nist.gov 

FCS,  Fermi-chopper  time-of-flight  spectrometer  (NG-6) 

C.  IVI.  Brown,  (301)  975-5134,  craig.brown@nist.gov 
T.  J.  Udovic,  (301)  975-6241,  udovic@nist.gov 

DCS,  Disk-chopper  time-of-flight  spectrometer  (NG-4) 
(NIST/NSF-CHRNS) 

J.  R.  D.  Copley,  (301)  975-5133,  jcopley@nist.gov 
I.  Peral,  (301)  975-6235,  inma@nist.gov 

HFBS,  High-flux  backscattering  spectrometer  (NG-2) 
(NIST/NSF-CHRNS) 

R.  M.  Dimeo,  (301)  975-8135,  robert.dimeo@nisLgov 
Z.  Chowdhuri,  (301)  975-4404,  zema.chowdhuri@nist.gov 

D.  A.  Neumann,  (301)  975-5252,  dan@nist.gov 

NSE,  Neutron  spin  echo  spectrometer  (NG-5) 
(NIST/NSF-CHRNS) 

N.  S.  Rosov,  (301)  975-5254,  nrosov@nist.gov 
D.  Bossev,  (301)  975-4662,  dbossev@nist.gov 

Prompt-gamma  neutron  activation  analysis  (NG-7) 

R.  M.  Lindstrom,  (301)  975-6281,  dick.lindstrom@nist.gov 
R.  L.  Paul,  (301)  975-6287,  rpaul@nist.gov 

Other  activation  analysis  facilities 

R.  R.  Greenberg,  (301)  975-6285,  rgreenberg@nist.gov 

Cold  neutron  depth  profiling  (NG-0) 

G.  Lamaze,  (301)  975-6202,  lamaze@nist.gov 

Instrument  development  station  (NG-0) 

D.  F.  R.  iVlildner,  (301)  975-6366,  mildner@nist.gov 

H.  H.  Chen,  (301)  975-3782,  chenmayer@nist.gov 

Neutron  interferometer  (NG-7) 

M.  Arif,  (301)  975-6303,  muhammad.arif@nist.gov 
D.  Jacobson,  (301)  975-6207,  jacobson@nist,gov 

Fundamental  neutron  physics  station  (NG-6) 

M.  S.  Dewey,  (301)  975-4843,  mdewey@nist.gov 

Theory  and  modelling 

N.  F.  Berk,  (301)  975-6224,  nfb@nist.gov 
T.  Yildirim,  (301)  975-6228,  taner@nist.gov 

Sample  environment 

D.  C.  Dender,  (301)  975-6225,  dender@nist.gov 
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NIST  Center  for 
Neutron  Research 


For  copies  of  or  information  on  this  report,  contact: 

Ronald  L.  Cappelletti 

(301)  975-6221 
ron.cappelletti@nist.gov 

For  additional  information  on  the  facility,  contact: 

J.  Michael  Rowe 

(301)  975-6210 
mil<e. rowe@nist.gov 

John  J.  Rush 

(301)  975-6231 
john.rush@nist.gov 

To  obtain  guidelines  for  preparing  proposals 
to  conduct  research  at  the  facility,  contact: 

William  A.  Kamitakahara 

(301)975-6878 
vi/illiam. kamitakahara@nist.gov 

Location  of  all  contacts: 

NIST  Center  for  Neutron  Research 
100  Bureau  Drive,  Mail  Stop  8562 
Gaithersburg,  MD  20899-8562 

Please  visit  our  vi^ebsite  for  further  information: 
wviM.ncnr.nist.gov 


